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I.  Introduction 

\ 

An  array  antenna  is  an  assembly  of  radiating  elements  in  one  of  many  possible 
geometrical  configurations  with  the  outputs  of  the  individual  elements  in  the  array 
combined  to  produce  radiation  of  desired  pattern  shape  and  gain.  The  advantages  of 
array  antennas  over  a  single  antenna  include  the  following:  1)  high  gain  and  high 
resolution  without  increasing  the  physical  size  of  each  element,  2)  electronic  control  of 
radiation  patterns,  such  as  scanning  of  main  beam  and  shaping  of  radiation  pattern,  and 
3)  graceful  degradation.  ^ 

^Many  techniques  have  been  developed  in  the  synthesis  of  array  patterns.  Most  of 
the  techniques,  however,  ignore  mutual  coupling  between  array  elements.  Mutual 
coupling  is  the  electromagnetic  interaction  between  array  elements.  In  addition,  unless 
the  antenna  elements  are  perfectly  matched  in  impedance  to  the  connected  transmission 
lines  and/or  perfect  isolation  power  dividers  are  used,  there  will  be  coupling  through  the 
feed  network.  The  two  fold  coupling  problem  (element-to-element  and  feed  network 
coupling)  then  becomes  difficult  for  analysis  and  very  difficult  for  synthesis.  This  is 
because  coupling  from  one  element  to  another  can  travel  through  the  feed  network  and  \ 
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reappear  in  other  elements  leading  to  further  antenna  coupling;  this  is  an  endless 
coupling  effect.  — - -  (  f-:  'f  - 

The  coupling  effects  are  minimum  in  the  synthesis  of  simple  patterns  such  as  sum 
and  difference  patterns  since  these  patterns  are  not  very  sensitive  to  the  element  current 
variations.  However,  in  cases  of  more  complex  patterns  such  as  low  side  lobe  patterns 
and  shaped  beams,  mutual  coupling  can  cause  significant  degradation  of  the  radiation 
pattern  from  the  desired  pattern. 

This  report  presents  new  methods  to  compensate  for  mutual  coupling  effects  of 
both  element  and  feed  coupling.  The  techniques  may  be  classified  as  synthesis  methods 
since  they  lead  to  realization  of  specific  current  distributions  across  the  array  in  the 
presence  of  practical  coupling  effects.  The  synthesis  techniques  described  in  this  report 
use  the  scattering  parameter  representation  of  array  antennas  and  were  initially 
formulated  by  Smith  and  Stutzman  [S-l].  Element  currents  can  be  calculated  in  the 
presence  of  mutual  coupling  (array  analysis);  also  a  network  that  compensates  for 
mutual  coupling  effects  cam  be  determined  (array  synthesis).  A  system  of  nonlinear 
equations  was  derived  by  Smith  and  solved  numerically  for  a  compensation  network. 
Two  approaches  were  used  to  compensate  for  mutual  coupling:  1 )  the  conventional  feed 
network  is  unaltered  and  amplitude  and  phase  shift  devices  are  placed  between  it  and  the 
array  at  each  port,  2)  the  parameters  of  the  feed  network  are  altered  slightly  in  the  design 
process  to  anticipate  compensation  for  coupling. 

This  report  marks  the  finish  of  a  four-year  project  supported  by  the  U.S.  Army 
Research  Office  through  the  University  of  Massachusetts.  Companion  experiments  were 
conducted  at  the  Physical  Science  Laboratory  of  New  Mexico  State  University  by 
Russell  P.  Jedlicka  and  his  colleagues.  The  problem  attacked  was  very  challenging, 
requiring  analysis,  numerical  computations,  and  experimental  verification.  Some  aspects 
of  the  experimental  results  are  presented  here  as  they  relate  to  comparison  with  theory. 
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However,  this  report  focuses  on  analytical  and  numerical  methods  of  array  synthesis. 
The  first  major  document  was  Smith  and  Stutzman  [S-1J.  Our  initial  focus  was  on 
microstrip  antenna  arrays  as  a  testbed  (hence,  the  project  title),  but  the  formulation  is 
completely  general  and  both  numerical  computations  and  measurements  involved  both 
microstrip  and  wire  antennas. 

Chapter  2  reviews  the  network  representations  of  array  antennas  and  the  general 
synthesis  technique  using  network  representations.  Also  in  Chapter  2  the  system  of 
nonlinear  equations  derived  by  Smith  will  be  reviewed.  Chapter  3  reviews  current 
literature  and  experiments  on  mutual  coupling  compensation.  In  Chapter  4  three  types 
of  compensation  networks  will  be  presented.  Chapters  5,  6  and  7  modify  the  nonlinear 
equations  for  each  type  of  compensation  network  so  that  the  equations  can  be  solved 
numerically.  In  Chapter  8  the  compensation  technique  is  applied  to  an  eight-element 
microstrip  array.  Chapter  9  presents  some  numerical  examples  using  a  moment  method 
computer  program  for  wire  antennas.  Chapter  10  gives  conclusions  and 
recommendations  for  future  study. 

Reader  note:  This  report  and  its  predecessor  [S-l]  are  lengthy  and  detailed.  Chapter 
10  of  this  report  was  written  to  aid  the  reader  in  extracting  the  contributions  and 
technical  results.  For  a  comprehensive  summary  of  the  findings  from  the  study  refer  to 
Chapter  10. 
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II.  Overview  of  Array  Modeling  and  Synthesis 


2.1  Scattering  Matrix  Representation  of  Array  Antennas 


An  array  antenna  in  its  most  general  form  consists  of  N-elements  connected  to  an 
N-way  feed  network  and  a  generator  as  shown  on  Fig.  2.1-1.  There  are  several  ways  to 
represent  these  multiport  networks  such  as  in  terms  of  impedance,  admittance  and 
hybrid  parameters.  The  mutual  coupling  in  a  network,  however,  is  most  easily 
represented  by  scattering  parameters,  or  S-parameters.  Unlike  other  parameters  which 
are  defined  in  terms  of  voltages  and  currents,  S-parameters  are  defined  in  terms  of 
incident  waves  {a,}  and  reflected  waves  {b,}  as  shown  on  Fig.  2.1-2.  S-parameters  are 
normalized  reflection  and  transmission  coefficients  under  the  condition  that  all  other 
ports  are  match  loaded  and  undriven  [G-l]: 


A  | 

ai  'h-o,  k*j 


(2.1-1) 


A  network  can  be  represented  with  S-parameters  in  matrix  form  as  follows: 
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Figure  2.1-1.  Network  representation  of  an  N-element  array  antenna. 
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Cb3  =  [SDCa] 


(2.1-2) 


where  [a]  and  [b]  are  vectors  of  incident  waves  and  reflected  waves  at  each  port,  and  [S] 
is  the  S-matrix  of  the  network.  These  wave  variables  are  related  to  voltages  and  currents 
by 

V,-^(a,  +  bi)  (2.1-3) 


(a,  -  to 


(2.1-4) 


where  Zo,  is  the  characteristic  impedance  of  the  port  i.  From  (2.1-3)  and  (2.1-4)  the 
average  power  flowing  into  port  i  is 


Pi-yRefyl*} 

-  Y  (a,a*  -  b,b* ) 


(2.1-5) 


It  is  often  convenient  to  use  normalized  voltages  and  currents,  where  the  voltage 
is  normalized  by  dividing  by  the  square  root  of  the  reference  impedance  Z,  and  the 
current  is  normalized  by  multiplying  by  the  square  root  of  the  reference  impedance. 
When  normalized  with  Z,  ■  Z«,  (2.1-3)  and  (2.1-4)  become 


'  ai  +  bj 

(2.1-6) 

®i  ”  bj 

(2.1-7) 

The  power  expression  (2.1-5),  however,  remains  the  same.  On  the  other  hand,  when  the 
characteristic  impedance  of  the  port  i  and  reference  impedance  are  different  (Z,  #  ZJ  a 
step  change  in  the  characteristic  impedance  at  port  i  must  be  assumed  before  the 
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voltages  and  currents  can  be  normalized  by  Z*  Throughout  this  report  the  normalized 
definition  of  currents  and  voltages  are  assumed. 

When  two  or  more  networks  represented  by  S-matrices  are  interconnected,  a  single 
S-matrix  that  describes  the  overall  network  can  be  formulated.  In  general,  q  ports  of  an 
M-port  network  A  and  an  N-port  network  B  may  be  connected  as  shown  on  Fig.  2.1-3. 
The  combined  network  is  described  by  the  matrix  equation 


'[bA]' 

’[SA] 

[0]' 

V]' 

[bB] 

o1 

LJ 

_ 1 

[SB]_ 

_[aB]_ 

The  order  of  the  scattering  matrix  in  the  above  equation  can  be  reduced  using  the 
procedure  discussed  by  Smith  (S-l,  Sec. 3. 1.2].  The  resulting  scattering  matrix  for  the 
overall  network  is 

[SA/B]  -  [SPP  +  SpJT  -  Scc)  lSCP]  (2.1-9) 


where  the  matrix  [SA/B]  is  a  square  matrix  of  the  order  N  +  M  —  2q  .  The  superscript 
A/B  represents  S-parameters  for  the  combined  network  that  consists  of  the  network  A 
and  the  network  B.  The  S-parameters  on  the  right  hand  side  of  (2.1-9)  are  submatrices 
of  original  S-matrix  in  (2.1-8)  grouped  in  external  (unconnected)  and  internal 
(interconnected)  ports.  These  matrices  are  defined  by  Gupta  [G-l]  as 


'[bP]' 

[SppJCSpc] 

'[aP]' 

.Od. 

.[Scp][Scc]_ 

.[ad. 

(2.1-10) 


where  [aP]  and  [bp]  are  the  wave  variables  of  the  external  ports,  and  [a<J  and  [bj  are 
wave  variables  of  the  internal  ports.  [Spp]  and  [Scc]  are  the  S-parameters  of  external  and 
internal  ports,  respectively,  and  [SpJ  and  [SCT]  are  the  S-parameters  from  the  external 


11.  Overview  of  Array  Modeling  and  Synthesis 


8 


Figure  2.1-3.  Arbitrary  connection  of  M-port  network  A  and  N-port  network  B. 
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ports  to  the  internal  ports  and  vice  versa.  The  matrix  [Spp]  is  a  square  matrix  of  order 
M  +  N  —  2q  ,  and  [Scc]  is  a  square  matrix  of  order  2q.  The  connection  matrix  [1“]  is 
defined  as  r„  —  1  if  port  i  is  connected  to  port  j  and  r„  —  0  otherwise.  For  example,  the 
submatrices  for  the  networks  A  and  B  in  Fig.  2.1-3  can  be  written  as 


[SPP] 


*q-H.q+l  aq+l,q+2 


>M.q+l 

0 


>M,q+2 

0 


...  S 


q+J.M 


>M.M 


n  B  oB 

aq+I,q+l  aq+1,q+2 


>N.q+l  aN,q+2  • 


.  s 


q+l.N 


.  s 


N.N 


(2.1-11) 


[Sec] 


M.i 


3q.i 

0 


>q.2 


0  0 


0 


'q.q 


>u 


0 


cB 


0 


0 


SB 

•• 


A  cB  cB  eB 

.  U  aq,l  aq,2  •••  aq,q 


(2.1-12) 


[Spd 


*q+l,t  °q+U 


>M.l 

0 


>q+l,q 


>M.q 


Vu  Jq+u 


>N,1 


1 

;B 

>N,2 


q+'.q 


>N,q 


(2.1-13) 
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where  row/column  i  for  i*  l...q  in  T  represent  ports  I  through  q  of  the  network  A  and 
row/column  i  for  i*q+1...2q  represent  ports  1  through  q  of  the  network  B.  For 
example,  rM+l  —  1  since  port  1  of  network  A  and  port  1  of  network  B  are  connected  as 
shown  in  Fig.  2.1-3. 

When  an  external  port  is  excited  by  a  generator,  (2.1-8)  and  (2.1-9)  can  be  extended 
to  allow  calculation  of  wave  variables  impressed  by  the  generator  at  the  remaining  port. 
Figure  2.1-4  shows  a  generator  connected  to  port  i  of  an  N-port  network.  In  general,  a 
generator  is  expressed  in  S-parameters  as 


-  M 
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(2.1-16) 


SGa°  +  c° 


where  S°  is  reflection  coefficient  of  the  generator  and  c°  is  the  wave  excited  by  the 
generator.  The  wave  variable  a°  is  the  incident  wave  into  the  generator  from  a  network 
connected  to  the  generator.  Equation  (2.1-16)  is  a  statement  of  the  fact  that  the  signal 
that  comes  out  of  the  generator,  bG,  is  a  superposition  of  waves  produced  by  the 
generator  and  waves  reflected  from  the  generator.  Ideally,  the  reflection  coefficient  of 
a  generator  is  small  (SG  =  0).  In  this  case  the  reflected  wave  from  the  generator  is  equal 
to  the  wave  induced  by  the  generator  (b°  ■*  c°).  The  multiport  combined  with  the 
generator  (as  in  Fig.  2.1-4)  can  be  expressed  as 


bf 
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(2.1-17) 


If  subscript  P  is  used  to  denote  external  ports  and  subscript  C  to  denote  internal  ports, 
(2.1-17)  is  of  the  form 


(.[be] 

The  above  equation  can  be  reduced  to  two  types  of  equations,  one  that  relates  the 
external  ports  and  one  that  relates  the  internal  ports.  The  relationship  at  the  external 


[sPP]  cspcjir capDi  r[or 

j  +  (2.1-18) 

[SCP]  [Sec]  M  Ccc] 
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ports  (for  example,  ports  1  through  i-1  and  i+1  through  N  in  Fig.  2.1-4)  can  be 
described  by 


[b A/G]  -  [SA/°][aA/G]  -l-  [cA/G]  (2. 1-19) 

where  [SA/°]  is  combined  S-parameters  for  the  network  and  the  generator  which  is  a 
square  matrix  of  the  order  N  —  1  obtained  using  (2.1-9).  The  wave  variables  [aA/c]  and 
[b^0]  are  equivalent  to  [aP]  and  [bP]  in  (2.1-18),  respectively.  The  waves  induced  by  the 
generator  at  external  ports  of  the  combined  networks  are  obtained  from  [G-l] 

[cA/G]  -  [Spc][r  -  See]'1  [cj  (2. 1-20) 

where  T  defines  the  connection  between  the  generator  and  the  feed  network.  The  second 
relationship  derived  from  (2.1-19)  describes  the  incident  waves  and  the  reflected  waves 
at  the  internal  ports  of  the  combined  network  (for  example,  port  i  and  generator  port 
of  the  network  in  Fig.  2.1-4)  [G-l]: 

[a<J  “  [T  -  Scc]  '[cj  (2.1-21) 


and 


[bd-CnW  (2.1-22) 

Suppose  the  network  of  Fig.  2.1-4  has  4  ports  and  port  3  is  excited  (N=4,  i=  3). 
The  waves  induced  by  the  generator  at  external  ports  of  the  combined  network  (ports 
1,  2,  and  4)  are  obtained  using  (2.1-20)  as  follows: 
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[cA/G]  -  [sPC][r  -  scc]_l[ccj 


(2.1-23) 


The  incident  waves  at  the  internal  ports  of  the  combined  network  (ports  3  and  5)  are 
obtained  using  (2.1-21) 


DkJ  * 


-  [r  -  SccT'ccj 


(2.1-24) 


The  reflected  waves  at  the  connected  ports  are  obtained  from  the  incident  waves  using 


(2.1-22) 


[be]  - 


-CHM 
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The  above  techniques  can  be  used  to  analyze  array  antennas.  In  particular,  the 
currents  at  each  antenna  element  in  the  presence  of  mutual  coupling  can  be  calculated 
in  five  simple  steps. 

1.  Find  the  S-parameters  of  the  antenna  array,  SA,  the  feed  network,  Sp,  and  the 
generator,  Sc.  This  may  be  done  by  measurements,  by  using  a  mathematical  model  for 
the  networks,  or  by  using  a  CAD  procedure  such  as  MCAP  (see  Section  3.2). 

2.  Calculate  the  combined  S-matrix  for  the  feed  network  and  a  generator,  SF/G  ,  using 
(2.1-9)  and  the  waves  induced  by  the  generator  at  the  output  ports  of  the  feed  network, 
cF'G,  using  (2.1-20). 

3.  Calculate  S-parameters  for  the  array/feed/ generator  network,  SA/F/G  using  (2.1-9). 

4.  Calculate  incident  and  reflected  waves  at  antenna  ports  for  the  array/feed/generator 
network  using  (2.1-21)  and  (2.1-22). 

5.  Calculate  the  element  currents  with  (2.1-7). 
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2.2  The  Mutual  Coupling  Problem 


Mutual  coupling  is  the  electromagnetic  interaction  between  two  antenna  elements. 
Mutual  coupling  can  significantly  alter  the  performance  of  an  array  antenna.  For 
instance,  it  can  change  the  active  impedance  of  the  antenna  elements  (the  input 
impedance  of  the  elements  when  the  array  is  fully  excited),  alter  the  element  excitation 
currents,  and  lead  to  deterioration  of  the  active  element  patterns.  All  of  these  changes 
lead  to  degradation  of  the  radiation  pattern  from  that  obtained  with  no  mutual  coupling. 
The  major  effects  on  the  radiation  pattern  are  changes  in  sidelobe  level,  and  shilling  and 
filling  of  nulls.  In  addition,  changes  in  gain  and  beamwidth  of  the  main  lobe,  and  other 
undesirable  effects  can  occur  due  to  mutual  coupling.  In  general,  these  effects  on  the 
radiation  pattern  increase  as  the  magnitude  of  coupling  increases. 

A  comprehensive  summary  of  mutual  coupling  studies  is  given  by  Smith  [S-l].  In 
this  section  the  mechanism  of  mutual  coupling  in  an  array  antenna  and  the  effect  of 
mutual  coupling  on  radiation  pattern  will  be  discussed.  In  particular,  we  will  consider 
mutual  coupling  on  microstrip  array  antennas.  The  basic  theory,  however,  is  applicable 
to  any  array  antenna  configuration. 

2.2.1  The  mechanism  of  mutual  coupling  in  microstrip  arrays 

There  are  two  types  of  mutual  coupling  that  exist  in  an  antenna  array:  the  mutual 
coupling  among  antenna  elements  and  the  coupling  through  the  feed  network.  Figure 
2.2-1  shows  the  coupling  mechanisms  in  a  four  element  array  antenna.  The  dominant 
mutual  coupling  among  the  antenna  elements  occurs  through  space  waves.  This  type  of 
coupling  is  strongly  dependent  upon  the  geometry  of  the  array.  The  interelement 
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spacing  is  the  most  important  parameter  in  space  wave  coupling.  The  mutual  coupling 
in  an  eight  element  microstrip  array  consisting  of  0.312>i  spaced  quarter- wave  patches 
which  is  used  in  one  of  the  experiments  performed  for  this  project  (described  in  Chapter 
8)  is  shown  in  Fig.  2.2-2.  The  magnitude  and  phase  of  the  scattering  coefficients  S„ 
between  element  1  and  j  are  plotted  versus  element  number  j.  Both  magnitude  and  phase 
of  the  mutual  coupling  show  strong  correlation  to  the  spacing  between  antenna 
elements.  This  result  agrees  with  the  mutual  coupling  studies  conducted  by  Jedlicka  and 
Carver  (J-l.J-2].  In  addition  to  the  element  spacing,  the  shape  of  antenna  elements  and 
their  orientation  are  important  in  the  space  coupling  [A-7^-1]. 

Secondary  mutual  coupling  in  a  microstrip  array  occurs  through  surface  waves  on 
the  array  substrate  [J-l].  In  this  coupling  the  substrate  acts  as  a  dielectric  waveguide 
permitting  signals  to  propagate  between  antenna  elements.  The  structure  of  the  array 
board  such  as  the  thickness  of  substrate,  the  permittivity  of  dielectric  material  and  the 
length  of  array  are  important  variables  in  surface  coupling.  However,  Jedlicka  and 
Carver  [J-1J-2J  have  shown  that  a  change  in  substrate  thickness  has  little  effect  in 
mutual  coupling.  It  implies  that  the  degree  of  surface  coupling  is  very  small  compared 
to  space  coupling.  Surface  wave  coupling  also  occurs  in  a  microstrip  feed  network  as  a 
secondary  mutual  coupling  effect. 

The  third  mutual  coupling  component  among  antenna  elements  is  return  loss 
represented  by  S„.  This  is  reflection  coefficient  of  an  element  when  all  other  elements  are 
match  loaded.  The  return  loss  does  not  alter  the  element  currents  when  there  is  no 
element  to  element  coupling  or  when  each  element  is  driven  by  an  independent 
generator.  The  most  significant  effect  in  these  cases  are  loss  of  power  efficiency. 
However,  when  mutual  coupling  exist  among  elements  and  in  the  feed  network  the 
return  loss  significantly  changes  element  current  distribution. 
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re  2.2-2.  Measured  magnitude  and  phase  of  mutual  coupling  as  a  function  of 
element  number  for  0.312/10  spaced  quarter  wave  microstrip  patches 
described  in  Section  8.1. 


The  primary  coupling  effect  in  the  feed  network  occurs  from  the  cross  talk  between 
two  ports  due  to  nonideal  power  dividers  used  in  the  network.  In  general,  the  cross  talk 
in  a  power  divider  cannot  be  completely  eliminated  without  use  of  active  devices.  It  is 
also  well  known  that  passive  multiport  networks  other  than  two  port  networks  cannot 
be  matched  on  all  ports.  This  implies  that  any  feed  network  for  array  antennas,  even  for 
a  simple  two  element  array,  cannot  have  output  ports  matched  to  the  antenna  elements. 

When  mutual  coupling  in  the  antenna  and  the  coupling  in  the  feed  network  are 
combined  very  complex  feedback  loops  are  created.  The  overall  effect  of  mutual 
coupling  on  the  element  excitation  currents  is  nonlinear.  This  will  be  illustrated  later  in 
Section  5.4. 

2.2.2  The  effect  of  mutual  coupling  on  the  radiation  pattern 

The  degradation  of  a  radiation  pattern  due  to  mutual  coupling  in  an  array  antenna 
includes  increased  side  lobes  and  shifting  and  filling  of  nulls.  These  effects  are  especially 
severe  when  precise  control  of  excitation  currents  is  necessary  as  in  the  case  of  a  low  side 
lobe  pattern.  The  effects  on  the  radiation  pattern  will  be  illustrated  here  using  numerical 
simulations. 

Consider  an  eight  element  array  with  0.312  wavelength  interelement  spacing.  We 
will  assume  the  measured  mutual  coupling  in  Fig.  2.2-2.  An  average  element  pattern  of 
005*^(0.90)  is  assumed.  Two  types  of  excitations  are  considered:  a  sum  pattern  and  a  35 
dB  Dolph-Chebyshev  pattern. 

The  synthesis  of  the  sum  pattern  requires  equal  magnitude  and  equal  phase 
excitation  on  all  antenna  elements: 
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(2.2-1) 


This  is  the  simplest  form  of  array  excitation.  An  8-way  corporate  feed  network  consists 
of  seven  3  dB  power  dividers  as  shown  in  Fig.  2.2-3  is  assumed.  The  S-parameters  of 
power  divider  feed  networks  are  derived  in  Section  6.2.  The  element  currents  in  the 
presence  of  mutual  coupling  are  calculated  using  the  Fortran  program  MCAP  described 
in  Section  4.2.  Using  MCAP  the  currents  are  obtained  as 


■ 


1 


f 


Ii  -  1.145 19/22.6°  I5  -  0.90666/1.9° 

I,  -  0.62643/ 11.5°  I*  —  1. 14473/29.9° 

-  -  (2.2-2) 

I3  -  0.96579/28.0°  I7  -  0.84026/26.6° 

I4~  0.8295 1/7JT  Ig  -  1.16862/11.1° 

The  plot  of  normalized  radiation  patterns  for  the  case  with  mutual  coupling  and  for  the 
case  without  mutual  coupling  are  shown  in  Fig.  2.2-4.  The  figure  shows  that  the  major 
effect  on  the  radiation  pattern  is  filling  of  nulls.  However,  other  effects  such  as  change 
in  side  lobe  levels  are  very  small. 

The  35  dB  Dolph-Chebyshev  array  requires  the  following  current  distribution  [S-2]: 

Is  =  1.00000/0.0° 

16  =  0.78425/0.0° 

(2.2-3) 

17  *  0.46362/0.0° 

1,-0.19154/0.0° 


1,-0.19154/0.0° 

12  =  0.46362/0.0° 

13  =  0.78425/0.0° 
I4=  1.00000/0.0° 
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The  feed  network  consists  of  an  8-way  corporate  feed  network  consists  of  seven  2-way 
variable  power  dividers  as  shown  in  Fig.  2.2-3.  The  power  split  ratio  K,  in  the  figure  are 
set  to  yield  desired  currents  when  there  is  no  coupling  in  the  array.  A  detailed  discussion 
on  N-way  power  divider  network  is  given  in  Section  6.2.  The  currents  in  the  presence 
of  mutual  coupling  obtained  using  MCAP  are 

I,  -  0.32639/-36.70  I5  -  Q.63545/-75.30 

I,  -  0.50824/- 1 1.3*  I*  -  0.993 16/-7 1 .4° 

(2.2-4) 

13  -  0.73737/-52.3*  I,  ■  0.279 18/- 11 3.9° 

14  -  1 .09558/-62.4*  I8  -  0. 1 1 985/-71.78 

The  calculated  radiation  pattern  for  both  sets  of  currents  (2.2-3)  and  (2.2-4)  are  shown 
in  Fig.  2.2-5.  The  figure  clearly  shows  that  the  side  side  lobes  increase  to  approximatly 
-25dB  for  the  array  with  mutual  coupling.  A  change  in  the  beamwidth  is  also  apparent. 
Compared  to  the  sum  pattern  this  low  side  lobe  pattern  requires  more  precise  control 
of  excitation  currents.  It  can  be  concluded  from  these  results  that  an  array  that  requires 
precise  control  of  excitation  currents  is  affected  more  by  mutual  coupling. 
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Figure  2.2-5.  Normalized  radiation  patterns  of  35dB  Dolph-Chebyshev  array  with 
main  beam  steered  to  -10  degrees,  a)  without  and  b)  with  mutual 
coupling  using  currents  from  (2.2-3)  and  (2.2-4)  and  a  cosOJ(O.90) 
element  pattern. 
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2  J  Array  Antenna  Pattern  Analysis  Approaches  Including  Mutual  Coupling 

Most  array  analysis  is  based  on  the  assumption  that  elements  are  connected  to 
independent  generators  with  either  forced  or  free  excitation  [O-IJ.  However,  array 
antennas  are  usually  excited  by  a  feed/generator  network,  in  which  case  the  above 
assumption  is  not  applicable.  We  begin  this  section  with  a  discussion  of  independent 
generator  excitation  methods,  and  follow  it  with  a  more  general  treatment  that  includes 
feed  network  effects. 

2.3.1  Independent  generator  excitations 

In  the  independent  generator  excitation  case  elements  of  the  array  are  assumed  to 
be  identical  and  to  be  driven  by  independent  generators.  It  is  only  necessary  to  consider 
the  coupling  among  antenna  elements. 

Forced  Excitation 

In  forced  excitation,  the  antenna  elements  are  fed  by  constant-voltage  sources  or 
constant-current  sources.  In  practice  it  is  difficult  to  create  such  sources  to  drive  the 
elements,  but  they  are  convenient  from  a  modeling  standpoint.  In  the  case  of  voltage 
sources,  as  shown  in  Fig.  2.3-1  the  generator  imposes  constant  voltages  across  the 
antenna  elements.  The  element  currents  are  then  found  from  the  voltage  and  the 
impedance  matrices  using 

DD-Czr’Cv]  (2.3-1) 
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where  Z„  represents  self  impedance  of  antenna  element  i  and  Z„  for  i  *  j  represents  the 
mutual  impedance  between  element  i  and  element  j.  The  radiation  pattern  for  this  array 
is  given  by  [K-I.W-l] 


N 


n»l 


(2.3-2) 


where 


—  xn  sin  6  cos  <f>  +  yn  sin  6  sin  0  +  z„  cos  0 

for  an  array  of  identical  elements  with  arbitrary  locations  (x,,,  y„,  zj.  The  coordinate 
system  used  here  is  defined  in  Fig.  2.3-2.  The  currents  {ij  are  obtained  from  (2.3-1)  and 
may  contain  a  linear  phase  that  steers  the  beam  to  location  (0„,  <£„).  The  isolated  element 
pattern,  g,(0,  <f>),  is  obtained  in  the  short  circuited  array  environment  where  only  one 
element  is  fed  by  a  generator  and  all  other  elements  are  shorted.  This  is  illustrated  in 
Fig.  2.3-3. 


Free  Excitation 

In  the  free  excitation  mode  each  antenna  element  is  driven  by  constant  incident 
power  sources  as  shown  in  Fig.  2.3-4.  This  is  a  more  realistic  method  of  feeding  array 
antennas  than  forced  excitation  method  [0-1].  Because  of  the  inclusion  of  the  generator 
impedances,  the  change  in  input  impedance  of  the  elements  alters  both  the  voltage  and 
current  in  each  element.  In  this  case  the  system  can  be  described  better  with  scattering 
parameters  as 

[b]  *  [S][a]  (2.3-3) 
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Figure  2.3-2.  Geometry  of  an  arbitrary  array  antenna 
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Figure  2.3-3.  Array  with  all  but  one  elements  shorted  (forced  excitation). 
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where  a's  and  b's  represent  incident  and  reflected  waves  for  the  array  as  discussed  in 
Section  2.1.  The  elements  of  the  scattering  matrix  [S]  represent  coupling  coefficients 
between  pairs  of  elements  of  the  array  with  all  others  match  loaded.  The  currents  at  each 
element  can  be  calculated  using  (2.1*7)  as 

[i']*[a]-[b]  (2.3-4) 

The  radiation  pattern  for  free  excitation  is  [A-1J 

N 

(2.3-5) 

n— 1 

where  g,  is  the  isolated  element  pattern  obtained  as  in  the  free  excitation  case.  In  (2.3-5) 
all  mutual  coupling  effects  are  contained  in  the  element  currents  (i'J. 

Alternatively,  the  radiation  pattern  for  the  free  excitation  can  be  obtained  from 
[A-l.W-U-3] 

N 

F(0,^)  =  X!g^(0,0)rnejJ"  (2.3-6) 

n—1 

where  g^,  is  the  active  element  pattern  for  the  n-th  element  which  is  obtained  when  only 
the  n-th  element  is  driven  and  all  other  elements  are  terminated  with  their  generator 
impedance  (see  Fig.  2.3-5).  I',  is  the  current  that  would  be  flowing  without  coupling 
from  other  generators  which  is  proportional  to  generator  voltage  Vn.  In  other  words,  the 
effect  of  mutual  coupling  is  contained  only  in  the  active  element  pattern  and  not  in  the 
element  currents  [A-l.J-3].  The  active  element  pattern  of  each  element  can  be  calculated 
from  (P-3, P-4] 
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(2.3-7) 


N 

8„(«,  *)  -  8,(0. 

m— 1 

where  ig*  is  the  current  in  the  m-th  element  induced  by  n-th  generator  when  all  other 
elements  are  match  loaded. 

The  relationship  in  (2.3-6)  is  one  of  superposition  and  relies  on  linearity  of  the 
array  environment,  which  is  the  usual  case.  The  mutual  coupling  here  arises  only  from 
scattering  from  the  array  elements.  This  is  different  from  the  fully  excited  array 
discussed  next. 

2.3.2  Full  coupling  array  representation 

Most  arrays  are  fed  by  a  feed/generator  network  which  cannot  be  represented  by 
independent  generators.  This  general  situation  is  shown  in  Fig.  2.3-6.  The  basic 
difference  between  the  independently  driven  array  and  the  feed  network  driven  array  is 
the  existence  of  coupling  through  the  feed  network  in  the  latter  case  as  well  as  coupling 
between  elements.  To  account  for  the  element  and  feed  coupling,  the  scattering 
parameter  analysis  discussed  in  Section  2. 1  must  be  used.  Specifically,  the  element  and 
feed  S-matrices  must  be  combined  using  (2.1-9),  and  solved  for  the  incident  and  reflected 
waves  at  the  elements  using  (2.1-21)  and  (2.1-22),  then  solved  for  the  element  currents 
using  (2.1-7).  The  radiation  pattern  can  be  obtained  using  (2.3-6)  but  with  currents  I„: 

N 

F(M) -£*.*(*,  (2-3*8) 

n»i 
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Here  g'^,(d,  <p)  is  the  pattern  of  each  element  in  the  array  environment  when  the 
array  is  operating  in  a  fully  excited  state.  These  element  patterns  cannot  be  measured 
directly,  rendering  (2.3-8)  somewhat  useless.  It  is,  however,  useful  as  a  step  in  the 
process  of  accounting  for  all  coupling  effects  in  a  fully  excited  array.  The  next  section 
reveals  this. 

2.3.3  Approximate  approaches  to  full  coupling  representation 

The  first  level  of  approximation  we  introduce  is  to  approximate  the  full  array  active 
element  patterns  in  (2.3-8)  by  the  active  element  pattern  of  (2.3-6)  obtained  by  exciting 
elements  one  at  a  time.  The  current  distribution  on  an  element  in  the  array  is  affected 
by  neighboring  elements.  This  includes  element  type,  orientation,  and  position  as  well 
as  the  neighboring  element  current  amplitude  and  phase.  All  these  factors  affect  the 
current  distribution.  However,  the  current  of  an  individual  element  is  affected  primarily 
by  its  own  boundary  conditions  and  secondarily  by  mutual  coupling  effects.  Therefore, 
we  assume  that  the  element  pattern  is  very  well  approximated  by  the  case  where  it  is 
excited  and  all  other  elements  are  match  loaded  (e.g.  the  active  element  pattern).  Thus, 

N 

(2-3-9) 

n«l 

The  next  level  of  approximation  in  array  analysis  is  to  factor  the  array  pattern  into 
a  product  of  an  array  factor  and  an  "element  pattern"  which  we  call  the  array  element 
pattern  g*(0,  <£).  The  array  element  pattern  is  the  pattern  of  a  typical  antenna  element 
obtained  by  measuring  its  pattern  with  all  other  elements  match  loaded.  The  element 
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usually  selected  is  a  central  element  in  the  array.  Then  the  array  pattern  of  (2.3-9)  is 
approximated  using 

N 

(2.3-10) 

n«»l 

where  the  summation  is  the  classical  array  factor.  The  current  I„  of  (2.3- 10)  is  obtained 
using  the  S-parameter  array  analysis  technique  discussed  in  Section  2.1.  Employment 
of  the  principle  of  pattern  multiplication  requires  that  all  elements  be  alike  (in  type  and 
orientation),  that  coupling  is  not  excessive  (elements  not  extremely  close  together),  and 
that  array  edge  effects  are  not  large. 

Further  study  on  derivation  of  (2.3-8)  and  approximations  used  to  obtain  (2.3-9) 
and  (2.3-10)  are  necessary.  In  particular  numerical  experiments  using  computer 
programs  such  as  ESP  (described  in  Sec.  9.2)  must  be  conducted  to  verify  these 
equations. 

2.4  The  Synthesis  Problem 

As  seen  in  Section  2.2,  mutual  coupling  alters  the  radiation  pattern  of  an  array 
antenna.  To  compensate  for  mutual  coupling  a  network  can  be  placed  between  the 
elements  and  the  feed  network.  The  network  is  designed  such  that  correct  currents  are 
obtained  at  the  input  to  the  antenna  elements  in  the  presence  of  mutual  coupling.  A 
more  useful  realization  of  a  compensation  effect  is  to  not  have  a  separate  compensation 
network,  but  to  include  mutual  coupling  compensation  into  the  design  phase  of  the 
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conventional  feed  network  which  can  be  handled  within  the  framework  of  the  general 
approach  introduced  here. 

In  the  most  general  form  an  array  has  N-elements,  an  M-way  feed  network 
connected  to  a  generator,  and  an  N  +  M  port  compensation  network  as  shown  on  Fig. 
2.4-1.  This  array  can  be  fully  described  using  S-parameters.  The  antenna  elements  and 
the  compensation  network  are  represented  by  matrices  [SA3  and  [Sc3,  respectively.  The 
combined  network  of  a  feed  network  and  a  generator  is  described  by  a  matrix  CSF/G]  and 
a  vector  CcF'°]  obtained  using  (2.1-9)  and  (2.1-20),  respectively. 

In  the  synthesis  of  a  compensation  network,  the  matrix  CSC3  represents  the 
unknowns,  and  the  desired  (design)  element  currents  are  denoted  as  CIAD  •  The  S-matrix 
of  the  elements,  CSA],  and  the  feed/ generator  network,  fSF/G],  are  assumed  to  be 
known.  The  relationship  between  the  incident  and  reflected  waves  of  the  compensation 
network  is  given  by 
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(2.4-1) 


where  the  notation  used  here  is  consistent  with  that  in  Fig.  2.4-1.  In  addition,  it  is  clear 
from  the  figure  that 


(2.4-2) 
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Figure  2.4-1.  General  form  of  an  array  antenna  with  a  compensation  network. 
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(2.4-3) 


for  1  ^  i  £  N.  The  wave  variables  [aA]  and  [b*]  are  related  to  the  (specified)  desired 
currents  [IA]  by 

[1A]  -  Oa]  -  [bA]  (2.4-4) 

or 

[aA]-[IA]  + [bA]  (2.4-5) 

The  wave  variables  are  also  related  to  the  antenna  S-matrix  [SA]  by 


[bA]-[SA][aA]  (2.4-6) 

Substituting  (2.4-5)  into  the  above  equation  and  solving  for  the  vector  [b*]  proceeds  as 
follows 


[bA] 


[SA][IA  -  bA] 


[bA]  -  [SA][IA]  -  [SA][bA] 
[I  -  SA][bA]  -  [SA][IA] 


[bA]  -  [I  -  SAr'[SA][IA]  (2.4-7) 

where  [I]  is  the  square  identity  matrix.  It  follows  from  (2.4-5)  and  (2.4-7)  that  wave 
variables  [aA]  can  be  represented  as 

[aA]  -  [IA]  +  [I  -  SA]"![SA][IA]  (2.4-8) 
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Equations  (2.4-7)  and  (2.4-8)  or  (2.4-5)  relate  the  wave  variables  to  the  known  S-matrix 
[SA]  and  the  desired  currents  [IA]  . 

At  the  junction  between  the  feed/ generator  network  and  the  compensation  network 
(see  Fig.  2.4-1), 

aS.N-bf/G  (2-4-9) 

b^-N  “  af/G  (2.4-10) 

for  1  ^  i  <,  M.  Using  (2.1-19),  the  wave  variables  at  the  feed  network  are  related  by 

[bF'°]  -  [sF'°][aF'°]  +  [cF'°]  (2.4-11) 

It  should  be  noted  that  the  incident  wave  vector  from  the  compensation  network  into 
the  feed/generator  network  Ca"G]  is  not  a  constant  vector.  Rather,  it  depends  on  the 
actual  compensation  network  CSC]  .  This  vector  is  not  necessary  in  the  final  solution. 
However,  these  variables  play  an  important  role  as  the  intermediate  unknowns  in  the 
solution  process.  To  emphasize  the  fact  that  the  incident  waves  need  to  be  treated  as 
unknowns,  the  vector  entries  are  changed  as  follows: 

U4  -  af/G  for  l£i£M  (2.4-12) 

Substituting  (2.4-2),  (2.4-3)  and  (2.4-9)  through  (2.4-12)  into  (2.4-1)  gives 
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(2.4-13) 


From  this  the  resulting  system  of  nonlinear  equations  is 


0 — 

i-l  i— I 


i-i 


for  1  <,  k  <.  N 


(2.4-14) 


which  follow  from  the  top  portion  of  (2.4-13)  and 
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N  M 

_Uk-M  +  + 

i-l  i-l 


JIG 


forN+l^k^N  +  M 


(2.4-15) 

which  follow  from  the  lower  portion  of  (2.4-13).  The  system  of  equations  (2.4-14)  and 
(2.4-15)  are  in  general  form.  The  equations  are  solved  for  S£  which  are  the  scattering 
parameters  of  the  compensation  network.  The  variables  a,,  b„  SfJG,  and  cF'°  are  known, 
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and  U,  are  intermediate  unknowns.  The  nonlinear  property  in  (2.4-14)  and  (2.4-15) 
arises  from  the  the  terms  S^U,  where  both  {S£}  and  {U,}  depend  on  unknown  variables 
{X,}.  When  generator  value  G  is  treated  as  an  unknown  (for  example,  power  divider 
compensation)  the  product  S^N+,c,F/c  also  contributes  to  the  nonlinear  property  of  the 
equations.  These  equations  cannot  be  solved  for  a  compensation  network  until  the 
relationship  between  physical  parameters  and  S-parameters  of  the  compensation 
network  is  defined.  In  other  words,  we  must  define  functions  of  the  form 

Sj-^X.X, . XN)  (2.4-16) 

where  {X,}  are  complex  variables  representing  network  parameters.  Three  different  types 
of  compensation  networks  considered  in  this  report  are  discussed  in  Chapter  4.  In 
Chapters  5  through  7  we  proceed  with  solving  for  required  compensation  configurations 
for  three  networks. 

The  system  of  nonlinear  equations  are  then  solved  for  compensation  network 
parameters  using  the  Damped  Newton's  method  (D-4J  which  is  an  iterati .  •'  numerical 
technique.  A  review  of  the  Newton's  method  and  the  Damped-Newton's  method  is  given 
by  Smith  [S-l,  Secs.  5.2-5.4).  The  numerical  technique  is  implemented  in  programs 
SANE  and  SANE-PODCON  to  solve  for  attenuator/phase  shifter  network  and  power 
divider  network  which  are  described  in  Appendix  A  and  Appendix  B,  respectively. 


2.5  Calculation  of  Element  Currents  at  the  Patch  Edges 

Microstrip  arrays  are  usually  fabricated  with  a  section  of  transmission  line  attached 
between  the  antenna  patches  and  the  board  input,  which  will  be  called  the  feed  lines  as 
indicated  in  Fig.  2.5-1.  These  feed  lines  should  not  be  confused  with  the  feed  network 
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Definition  of  antenna  patches  and  their  feed  lines. 


that  is  fabricated  on  a  separate  substrate.  The  purpose  of  the  feed  lines  are  to  excite  the 
antenna  patches  with  correct  mode  and  polarization. 

The  S-parameters  of  the  array  board  represent  the  combined  network  of  antenna 
patches  and  their  feed  lines.  The  currents  calculated  from  the  antenna  and  feed  network 
matrices  using  (2.1-7),  (2.1-21)  and  (2.1-22)  are  the  values  at  the  feed  line  -  feed  network 
junction.  In  the  absence  of  coupling  these  currents  are  equivalent  to  the  array  excitation 
coefficients.  However,  in  the  presence  of  mutual  coupling  each  antenna  element  in  an 
array  has  a  different  active  impedance.  As  a  result,  the  relative  currents  at  the  patch 
edges  are  different  from  that  of  the  feed  line  input.  To  determine  the  correct  currents 
at  patch  edges,  the  feed  line  currents  need  to  be  transformed. 

Two  types  of  current  transformations  are  necessary  for  analysis  and  for  synthesis. 
In  the  analysis  case,  the  calculated  currents  at  the  board  input  need  to  be  transformed 
to  the  currents  at  patch  edges.  In  the  synthesis  case,  on  the  other  hand,  the  desired 
currents  at  patch  edges  need  to  be  transformed  to  the  board  input  currents  to  compute 
the  compensation  network. 

In  either  case,  it  is  necessary  to  assume  that  S-parameters  of  feed  lines  [SFL]  are 
known.  Since  the  feed  lines  cannot  be  separated  from  the  antenna  patches  to  measure 
S-parameters,  an  accurate  mathematical  model  for  the  feed  line  is  required.  In  this 
report  the  S-parameters  obtained  from  MCAP  are  assumed  to  be  accurate  enough  since 
the  structure  of  feed  line  is  simple.  It  is  also  necessary  to  assume  that  the  currents  at 
patch  edges  actually  correspond  to  the  excitation  coefficients  of  the  antenna  patches. 
More  rigorous  analysis  of  the  patch  -  feed  line  interaction  can  be  accomplished  using  the 
method  of  moment  as  reported  by  Pozar  [P-9]. 

In  the  analysis  case  the  wave  variables,  the  incident  waves  a  and  reflected  waves  b 
are  known  for  the  specified  S-matrices  of  an  antenna  array,  the  feed  network  and 
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generator  by  (2.1-21)  and  (2.1-22).  The  wave  variables  at  the  patch  edge  can  be  obtained 
from  the  general  S-parameter  equation: 


CM  -  [SFL]Ca3 


(2.5-1) 


or 


binput 

boutput 


CFL  -FL 
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(2.5-2) 


where  a,^  and  blnpu(  are  wave  variables  at  board  input  which  are  known,  and  a^^,  and 
bOTIp«t  are  wave  variables  at  the  feed  line  near  the  patch  edge.  The  matrix  [SFL]  represents 
S-parameters  of  a  feed  line  which  may  include  loss  in  the  line.  It  is  clear  from  Fig.  2.5-2 
that  the  incident  waves  to  the  patches  are  the  reflected  waves  from  the  feed  line  output 
and  vice  versa,  i.e.,  b*^,  =  a^  and  a^^  =*  b^.  Using  these  relationships  in  (2.5-2) 
and  solving  for  the  wave  variables  at  patches, 
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Then,  the  current  at  the  patch  edge  is 


(2.5-4) 
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(binput  ^11  ainput) 


In  the  synthesis  case  a  similar  transformation  cannot  be  applied  directly  because 
only  the  desired  currents  and  not  wave  variables  are  specified  at  the  patch  edges.  It  is 
clear  from  (2.1-7)  that  there  is  an  infinite  number  of  pairs  of  a„  and  b„  which  produce  the 
specified  value  of  currents  I„.  To  calculate  wave  variables  at  the  patch  edges,  the 
scattering  matrix  for  antenna  patches  without  feed  lines  is  necessary.  This  can  be 
obtained  by  subtracting  the  feed  line  portion  of  S-parameters  from  the  S-matrix  of 
antenna  patches  with  feed  lines. 

Consider  the  problem  of  removing  a  transmission  line  section  of  known 
S-parameters  from  the  port  1  of  a  general  N-port  as  indicated  in  Fig.  2.5-3.  We  will 
designate  the  original  N-port  scattering  matrix  by  [Sm]  and  the  reduced  N-port  by 
[S(1>].  The  feed  line  section  has  an  S-matrix  of  the  form 
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Using  the  interconnection  matrix  equation  (2.1-9), 

[s(0)]  —  [sPp  +  sPC(r  —  scc)  1  sCP] 


(2.5-7) 


where 
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Sp,  represents  the  scattering  matrix  for  the  network  with  external  ports  1,2,. ..,N  and  Scc 
represents  the  connected  ports.  The  connection  of  ports  2°-  and  1'  is  expressed  by 

cn-r  M  (25-9) 

Evaluation  of  the  right  hand  side  of  (2.5-7)  proceeds  as  follows: 


(r  -  sccr' 
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Equating  (2.5-10)  to  the  original  S-parameter  matrix,  we  can  obtain  four  relations. 
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Solving  the  above  four  equations  for  the  S-parameters  S®  gives 
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Generalizing  (2.5-15)  to  (2.5-18)  for  i-th  feed  line, 
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To  obtain  the  S-matrix  of  antenna  patches  without  feed  lines,  the  above  process 
must  be  applied  sequentially  to  all  feed  lines;  see  Fig.  2.5-4.  Thus,  the  S-matrix  for  the 
array  shown  at  the  bottom  of  Fig.  2.5-4  is  [S®0]  with  entries  given  by  (2.5-19)  to  (2.5-22) 
and  evaluated  sequentially  for  i»  1,...,N.  The  wave  variables  can  now  be  calculated  for 
the  antenna  patches  using  (2.4-7)  and  (2.4-8),  and  desired  currents  referred  to  the  board 
input  are  obtained  from  a  process  similar  to  the  analysis  case. 

Subroutines  were  developed  in  the  computer  programs  SANE  (Appendix  A)  and 
SANE-PODCON  (Appendix  B)  implementing  the  analysis  and  synthesis  steps  discussed 
in  this  section.  In  particular,  equations  (2.5-19)  through  (2.5-22)  which  calculate 
S-matrix  of  microstrip  patches  given  S-matrix  of  an  array  board  with  feed  lines  are 
implemented  in  subroutine  PREPRC  in  both  SANE  and  SANE-PODCON. 
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III.  Existing  Methods  to  Compensate  for  Mutual 
Coupling 


Many  techniques  have  been  developed  to  compensate  or  partially  compensate  for 
mutual  coupling  effects  in  array  antennas.  Most  of  these  techniques,  however,  consider 
only  the  element  coupling  and  neglect  the  second  order  effect  of  feed  coupling.  In 
addition,  many  techniques  are  applicable  only  for  a  specific  array  geometry  or  for 
specific  radiating  elements.  This  chapter  presents  a  review  of  some  existing 
compensation  techniques. 


3.1.  Compensation  by  an  Added  Network 


Perhaps  the  most  popular  method  to  compensate  for  mutual  coupling  is  to  place 
a  network  between  the  array  element  network  and  feed  network.  The  simplest  of  such 
networks  are  circulators  and  isolaters  [H-1 1].  However,  these  networks  have  power  loss 
in  reflection  [H-4]  and  are  non-reciprocal  causing  difference  on  transmit  and  receive. 
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They  are  also  expensive  components.  A  more  practical  method  to  compensate  for 
mutual  coupling  is  to  place  a  matching  network  between  the  array  and  the  feed  network. 

Hannan,  et  al.  [H-4]  developed  a  theory  that  accomplishes  impedance  matching 
using  networks  that  connect  feed  lines.  The  theory  was  derived  from  the  proof  by 
Hannan  (H-3)  that  for  an  infinite  array  it  is  possible  to  match  all  elements  over  all  scan 
angles  for  a  single  frequency  and  polarization. 

Andersen  and  Rasmussen  [A-2]  presented  an  example  for  synthesis  of  a  connecting 
network  applied  to  a  finite  array.  Their  technique  uses  a  network  that  has  the  negative 
of  the  impedance  matrix  of  the  array  elements.  The  network  is  connected  in  series  with 
the  elements  which  reduces  the  impedance  matrix  of  the  overall  network  to  zero.  For 
the  connecting  network  to  be  lossless  and  passive  it  requires  that  the  element  impedance 
matrix  be  purely  reactive.  This  is  accomplished  by  changing  the  array  geometry.  A  two 
element  experiment  showed  a  verification  of  the  technique.  However,  for  a  three 
element  array  complete  compensation  was  not  achieved  because  reactive  mutual 
impedance  between  two  outer  elements  and  adjacent  elements  cannot  be  obtained 
simultaneously. 

Davies  [D-3]  presented  an  application  of  a  matching  network  at  the  input  ports  of 
a  Butler  matrix  network.  He  showed  that  for  the  case  of  small  circularly  symmetric 
planar  arrays  such  as  those  used  as  feeds  to  reflector  antennas,  mutual  coupling 
compensation  can  be  accomplished  by  a  set  of  fixed  phase  shifts.  This  is  achieved  by 
breaking  down  the  array  excitation  into  series  of  orthogonal  circularly  symmetric  phase 
modes,  and  impedance  matching  phase  mode  ports  with  phase  shifters.  He  claimed  that 
these  networks  are  independent  of  changes  in  the  excitation  which  are  necessary  to  scan 
the  main  beam. 

The  compensation  techniques  discussed  above  are  all  based  on  the  knowledge  of 
coupling  coefficients  between  antenna  elements.  Another  approach  to  mutual  coupling 
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compensation  is  to  use  the  radiation  pattern  of  the  array.  Davis  [D-9J  showed  that 
element  currents  in  an  N-element  array  can  be  approximately  determined  from  the 
magnitude  and  phase  of  the  radiation  pattern  at  N  points.  This  calculation  assumes  that 
active  element  patterns  are  similar  for  all  elements.  Using  the  calculated  currents, 
attenuators  and  phase  shifters  are  added  between  the  elements  and  the  feed  network  to 
correct  for  the  error  in  the  currents.  As  it  will  be  shown  later  in  Section  5.4  this  process 
is  nonlinear.  Thus,  the  correction  to  the  element  currents  must  be  done  iteratively. 
Davis  applied  the  pattern  analysis  technique  to  synthesize  a  35  dB  Dolph-Chebyshev 
pattern  using  an  8-element  array.  Using  the  technique  the  side  lobe  level  were  reduced 
from  -26  dB  to  -31  dB  after  several  iterations. 


3.2.  Compensation  by  Feed  Network  Modification 

The  methods  treated  in  Section  3.1  compensated  for  mutual  coupling  by  adding  a 
network  between  the  elements  and  feed  network.  In  this  section,  we  review  methods 
which  allows  synthesis  of  feed  networks  that  include  mutual  coupling  compensation. 
Kang  and  Pozar  [K-2]  developed  two  techniques  to  determine  terminal  voltages  required 
to  produce  the  desired  element  currents  using  moment  method  impedance  matrices  of 
array  antennas.  Then,  the  calculated  voltages  and  desired  currents  can  be  used  to  design 
a  feed  network.  The  first  technique  compensates  for  mutual  coupling  by  point  matching 
a  desired  pattern  in  a  given  plane.  A  more  complicated  second  technique,  on  the  other 
hand,  corrects  the  pattern  over  all  space.  They  show  a  numerical  example  of  synthesis 
of  a  30  dB  Dolph-Chebyshev  endfire  pattern  using  an  array  of  eight  parallel  dipoles. 
Application  of  these  techniques,  however,  may  be  difficult  for  some  radiating  elements 
since  the  techniques  require  a  moment  method  impedance  matrix  of  the  array.  Also  the 
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techniques  neglect  the  feed  network  coupling  effect  which  may  be  significant  in  a  low 
side  lobe  pattern. 

Another  technique  presented  by  Elliott  and  Stem  (E- 1]  compensated  for  mutual 
coupling  using  microstrip  dipoles  which  are  fed  with  electromagnetically  coupled  feed 
lines.  In  this  technique  the  spatial  relation  between  the  dipoles  and  the  feed  lines,  the 
length  of  dipoles,  and  location  of  dipoles  are  design  variables.  The  self  impedance  and 
mutual  impedance  of  the  dipoles  are  functions  of  the  design  variables  which  were 
determined  by  measurements.  Using  these  functions,  a  multivariable  equation  can  be 
solved  for  the  design  variables  given  a  desired  radiation  pattern,  and  then  realized  by  the 
geometry  of  the  array.  Experimental  verification  of  the  theory  was  achieved  with  a  32 
element  dipole  array  [S-3]. 

3.3.  Other  Mutual  Coupling  Compensation  Techniques 


There  are  many  other  compensation  techniques  that  are  specific  to  array 
architecture.  For  example,  mutual  coupling  compensation  can  be  accomplished  by 
software  for  a  digital  beam  forming  (DBF)  array  (S-6.H-5].  In  a  DBF  array,  the  received 
signals  are  detected  and  digitized  at  the  element  level.  The  digitized  signals  are  then 
processed  using  digital  computer  to  form  a  desired  radiation  pattern.  Steyskal  [S-6] 
showed  that  when  the  amplitude  and  phase  of  signals  at  each  element  in  an  array  as  well 
as  coupling  coefficients  between  elements  are  available  such  as  in  a  DBF  array,  coupling 
efTects  can  be  cancelled  by  taking  the  inverse  of  the  coupling  coefficient  matrix  and  then 
multiplying  the  inverted  matrix  with  the  signals  at  each  element.  Herd  [H-5]  used  this 
technique  for  an  eight  element  array  to  reduce  the  side  lobe  level  from  -20dB  to  -30dB. 
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IV.  Feed  Configurations 


4.1  Introduction 

To  compensate  for  mutual  coupling  the  system  of  nonlinear  equations  derived  in 
Section  2.4  must  be  solved.  These  equations,  however,  cannot  be  solved  until  the  form 
of  the  scattering  matrix  CSCD  is  specified  for  the  physical  configuration  of  the  network. 
In  this  study  we  considered  three  types  of  feed  network  architectures  that  involve  the 
following  basic  building  blocks:  attenuator/phase  shifter  network,  variable  power  divider 
network,  and  transmission  line  network.  This  chapter  develops  the  mathematical 
representations  necessary  to  solve  for  compensation.  In  general,  these  networks  will  be 
considered  as  ideal  in  the  sense  that  ohmic  loss  in  transmission  lines  and  phase  shifters 
are  ignored,  and  mismatch  in  the  phase  shifters  and  attenuators  are  also  ignored. 
However,  for  the  attenuator/phase  shifter  network,  the  non-ideal  parameters  of  the 
device  will  also  be  discussed. 
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4.2  Using  MCAP  to  Study  Feed  Components 

The  Microwave  Circuit  Analysis  Program  (MCAP)  which  was  originally  written 
at  Indian  Institute  of  Technology,  Kanpur,  is  a  general  purpose  program  for  microwave 
network  analysis  [G-l].  The  program  has  been  modified  at  New  Mexico  State  University 
Physical  Science  Laboratory  to  enhance  its  capability.  The  modifications  are  discussed 
in  detail  by  Henry  and  Jedlicka  [H-l].  Many  such  microwave  CAD  codes  exist,  but  this 
one  is  available  at  no  cost  and  is  documented  in  [G-l).  Access  to  the  source  code 
permitted  the  modifications. 

MCAP  evaluates  the  scattering  matrix  of  an  overall  network  given  the 
S-parameters  of  individual  components  in  the  network.  This  is  basically  a  computer 
implementation  of  the  interconnected  network  equation  in  (2.1-9).  The  program  is 
equipped  with  mathematical  models  of  various  stripline  and  microstrip  components. 
These  models  evaluate  the  S-parameters  of  the  components  from  their  physical 
dimensions  and  their  characteristic  impedance.  The  standard  components  available  in 
MCAP  include  transmission  line  section,  step  in  the  line  width  (impedance  step),  bend 
in  a  line  at  an  arbitrary  angle,  T-junction,  coupled  lines,  series  impedance  and 
admittance,  and  shunt  impedance  and  admittance.  In  addition,  user  supplied  S-matrices 
for  networks  with  up  to  10  ports  can  be  specified  as  components.  After  the  evaluation 
of  the  overall  S-matrix,  the  program  calculates  the  wave  variables  and  corresponding 
currents  at  all  internal  ports  when  an  external  port  is  connected  to  a  generator  with  the 
value  G  =  1/0^  and  all  other  external  ports  are  match  loaded  which  is  an  implementation 
of  (2. 1-21)  and  (2.1-22). 

MCAP  was  used  extensively  in  our  study  to  analyze  antenna  element  currents  from 
measured  S-parameters  of  the  antenna  elements  and  feed  networks.  For  example,  an 
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eight  element  array  can  be  treated  as  a  combination  of  elements,  a  feed  network  and  a 
generator.  The  internal  ports  for  this  network  are  the  element  input  ports  and  the  feed 
network  output  ports  which  are  connected  to  each  other.  The  element  currents  are 
calculated  at  the  junctions  between  the  element  ports  and  feed  network  output  ports 
using  MCAP.  The  program  was  also  used  to  investigate  and  verify  the  S-parameters  of 
compensation  networks. 


4.3  The  Attenuator/Phase  Shifter  Compensation  Network 

The  basic  configuration  of  the  network  consists  of  an  attenuator  in  series  with  a 
phase  shifter  as  shown  on  Fig.  4.3-1.  The  attenuator  reduces  the  magnitude  and  the 
phase  shifter  alters  the  phase  of  the  signal.  The  series  combination  of  these  two  devices 
provides  the  required  complex  control  of  the  signal  flow  to  compensate  for  mutual 
coupling.  The  scattering  matrix  of  an  ideal  attenuator/phase  shifter  is  given  as  [S- 1} 

[SA/P]  - 

where  A  is  the  attenuation  and  0  is  the  phase  delay.  This  equation  assumes  that  both 
the  attenuator  and  phase  shifter  are  matched,  so  that  no  reflection  occurs  (S„  =  0).  It 
also  assumes  that  the  attenuator  provides  only  attenuation  and  no  phase  delay,  and  that 
the  phase  shifter  provides  only  phase  shift.  In  the  subsequent  analysis,  the  ideal 
performance  via  (4.3-1)  is  assumed.  However,  representation  of  non-ideal  performance 
is  discussed  next  for  completeness. 

In  practice,  neither  the  attenuators  nor  the  phase  shifters  are  ideal  devices.  The 
attenuator  introduces  some  phase  shift  and  mismatch.  Similarly,  the  phase  shifter  has 
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Attenuator 


Phase  Shifter 


Figure  4.3-1.  Graphical  representation  of  attenuator/phase  shifter  compensation 
network.. 
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some  attenuation  and  mismatch.  The  S-matrix  of  attenuator/phase  shifter  which 
includes  these  non-ideal  parameters  was  derived  by  Smith  [S-l]  as  follows 


[S££  ideal] 


rtB2e-i2^  "I 

p  rpra-i  J 

-ABe~**‘ +  *** 

rPra- 1 


-ABe'i(^^i>) 

ry\-i 


r pA2e-i2^'  ' 

rpr.-i 


(4.3-2) 


where  A  is  the  attenuation,  <f>t  is  the  phase  shift  and  T,  is  the  reflection  coefficient  of  the 
attenuator.  Similarly,  B  is  the  attenuation,  <f>9  is  the  phase  shift  and  Tp  is  the  reflection 
coefficient  of  the  phase  shifter.  When  mismatches  of  these  devices  are  small  enough 
such  that  they  can  be  neglected  (T,  =  r„  =  0),  the  S-matrix  reduces  to 


C  ^no  mismatch 


0 


ABe-^* +  *** 
0 


(4.3-3) 


From  this  we  see  that  phase  shifter  loss  can  be  lumped  with  the  attenuator  and 
attenuator  phase  shift  can  be  lumped  with  the  phase  shifter.  Interconnecting 
transmission  loss  and  phase  shift  can  be  handled  similarly. 

In  the  experiment  associated  with  this  project,  the  S-matrix  of  the  form  (4.3-3)  was 
utilized  to  realize  required  attenuation  and  phase  shift  by  the  combination  of  an 
attenuator  and  a  section  of  coaxial  cable.  It  is  usually  difficult  to  obtain  an  attenuator 
with  the  exact  value  required  for  the  compensation.  One  can  use  precision  variable 
attenuators  to  obtain  exact  values.  However,  they  are  not  practical  in  a  permanent 
configuration.  In  the  compensation  experiment  associated  with  this  project  a  fixed 
attenuator  was  chosen  to  approximate  the  required  attenuation.  Then  a  section  of 
coaxial  cable  was  cut  to  realize  the  required  phase  shift  and  remaining  attenuation. 
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4.4  Variable  Power  Divider  Network 


The  variable  power  divider  network  compensates  for  the  mutual  coupling  in  a  very 
different  way  from  the  other  two  methods  discussed  in  this  report.  Unlike  the 
attenuator/phase  shifter  (Sec.  4.3)  or  transmission  line  (Sec.  4.5)  networks,  the  power 
divider  network  is  itself  a  feed  network.  In  this  compensation  technique,  the  feed 
network  and  the  compensation  network  are  synthesized  simultaneously.  Only  a 
knowledge  of  the  S-matrix  of  the  antenna  elements  is  required  to  solve  for  the 
compensation  network. 

A  general  representation  of  a  2-way  variable  power  divider  network  is  shown  in  Fig. 
4.4-1.  The  network  consists  of  two  parts,  a  variable  power  divider  and  a  section  of 
transmission  line.  The  amplitude  control  of  the  output  signal  is  achieved  through  power 
division  instead  of  attenuation.  A  separate  control  of  phase  is  necessary  since  the  power 
dividers  are  phase  balanced.  This  network  has  a  definite  advantage  over  the 
attenuator/phase  shifter  network  from  an  efficiency  standpoint  since  power  dividers  are 
ideally  lossless. 

4.4.1  S-matrix  representation  of  power  divider  network 


In  order  to  solve  for  the  power  divider  compensation  network  parameters,  it  is 
necessary  to  know  the  S-parameters  of  the  network  in  terms  of  K,  the  fraction  of  input 
power  transmitted  to  output  port  2.  We  will  assume  that  all  input  power  is  split  between 
the  two  output  ports  as  follows: 
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Figure  4.4-1.  Graphical  representation  of  the  2-way  power  divider  network. 
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where  P,  and  P2  are  the  output  powers  from  port  1  and  2,  respectively,  and  Ptn  is  the 
input  power  to  port  3.  The  power  flowing  into  port  i  is  related  to  the  wave  variables  as 

P i  “  aA*  -  bjb*  (4.4-2) 

Then  (4.4-1)  can  be  expressed  as 

bib,  —  a, a,  =*  (1  —  K)(a3a3  —  b3b3) 

.  .  .  .  (4.4-3) 

bjbj  a2a2  *■  K(a3a3  b3b3) 


In  general,  a  three  port  device  such  as  a  2-way  power  divider  has  an  S-matrix  of  the  form 


S|1  S,2  S13 
S2,  S22  S23 
S31  S32  S33 


(4.4-4) 


where  [b]  *  [Sro][a].  We  assume  that  port  3,  the  input  port,  is  matched  so  that 
Sj,  =»  0.  When  port  1  and  port  2  are  match  loaded,  a,  =  0  and  a,  =  0,  (4.4-4)  becomes 
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and  (4.4-3)  becomes 

b,b,  *  (1  —  K)a3a3 
b2b2  =  Ka3a3 

since  bj  *  0  from  (4.4-5).  Using  (4.4-6)  equation  (4.4-5)  can  be  expanded  as 
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(4.4-7) 
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With  these  results  and  the  application  of  the  reciprocal  property  of  power  dividers, 
(4.4-4)  becomes 


S„  S,2  Jl-K 

Sl2  S22  JK  e^33 
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(4.4-8) 


The  remaining  four  S-parameters  can  be  solved  for  using  the  unitary  property  of  lossless 
networks  which  is 


for  j »  k 
for  j  #  k 


For  j*  1  and  k-3, 

Sn^l-Ke^13  +  SuViTe^33  *  0 


Solving  for  Su  gives 
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Next,  for  j=*k*  1  in  (4.4-9), 
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Substituting  (4.4-11)  into  this  gives 


=  K 

S„S„  4  SjjS]]  ■»  K 

KSnSjj  4-  (1  —  K)SUS|,  *  K2 
S,iS„  =  K2 

S„  -  Ke**“  (4.4-13) 

Using  this  in  (4.4-1 1), 

S,2  -  -  v/K(l^-  K)  (4.4-14) 

«• 

Similarly,  for  j=*k=2, 

S12S12  4-  S22S22  4  K  =»  1  (4.4-15) 

Using  (4.4-14)  in  this  and  solving  for  Sa, 

S22S22  =*  1  —  K  —  St2S|2 

-  1  -  K  -  K(1  -  K) 
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Finally,  for  j  =  2  and  k*  3  in  (4.4-9)  and  using  (4.4-7)  and  (4.4-14), 

S21S31  +  S22S32  ™  0 

—  VfT(l  —  K)t*+" " ' +*>  +  (1  -  K)jK  ~+»>  =  0 

4>\\  ~  ^23  “  <^22  ~  023 

0u“  022  (4.4-17) 


Combining  the  forgoing  results  of  (4.4-13),  (4.4-14),  (4.4-16)  and  (4.4-17),  the  scattering 
matrix  of  the  2-way  power  divider  of  (4.4-8)  becomes 
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(4.4-18) 


Adding  a  section  of  transmission  line  of  length  t  to  port  2  of  the  power  divider  for  phase 
control  as  indicated  in  Fig.  4.4-1  yields  the  complete  S-matrix  for  the  variable  power 
divider  network 


Ke^11 

—  ^K(l  —  K)  e^*11  +  ~ 

J  l  -  K  e^13 


IV.  Feed  Configurations 


71 


-  VK(1  -  K)  t}f*ii  +  *n-*a-M  J\-K  fi*" 
(1  -  K)^"  ~2p/)  ~pe) 

0 


(4.4-19) 


4.4.2  Physical  realization  of  the  variable  power  divider 


Many  types  of  power  dividers  using  microstrip  technology  have  been  discussed  in 
the  literature  and  have  been  implemented.  These  power  dividers  can  be  categorized  into 
two  groups.  In  1960  Wilkinson  {W-4j  implemented  an  N-way  power  divider  using  strip 
line.  Unfortunately,  this  device  provides  only  equal  power  split  among  the  output  ports. 
For  mutual  coupling  compensation  variable  power  dividers  must  be  used.  Hanna 
[H-6.H-7]  has  designed  an  N-way  variable  power  divider  by  using  a  Wilkinson  power 
divider  and  attenuators.  This  design  has  inherent  power  loss  in  the  network  due  to 
attenuators  and,  thus,  is  not  practical  as  a  power  divider  compensation  network.  A 
general  N-way  variable  power  divider  does  not  seem  to  be  practical. 

Variable  2-way  power  dividers,  on  the  other  hand,  are  physically  realizable.  Branch 
guide  couplers  and  in-line  power  dividers  discussed  by  House  [H-2],  and  hybrid  ring  and 
modified  hybrid  ring  power  dividers  [A-8.A-9]  are  some  of  the  common  designs  for 
variable  power  dividers.  In  general,  these  power  dividers  realize  the  unequal  output 
power  split  by  the  ratio  of  characteristic  impedances  of  transmission  lines  in  the  device. 
These  devices  are  adequate  when  the  power  ratio  required  is  relatively  small.  In  the  case 
of  large  power  split  ratio  the  narrow  line  width  required  for  the  high  characteristic 
impedance  or  the  wide  line  width  required  for  the  low  characteristic  impedance  would 
be  impractical  to  implement. 
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The  quadrature  hybrid  power  divider  discussed  by  Bailey  [B-2]  uses  the  ratio  of  the 
length  of  transmission  lines  for  control  of  power  division.  See  Fig.  4.4-2.  In  this  device 
an  unequal  output  power  ratio  is  obtained  by  exciting  the  input  p-irts  of  the  quadrature 
hybrid  with  signals  of  equal  amplitude  but  different  phases.  The  phase  adjustment  is 
made  by  varying  input  line  length  .  Then,  the  signal  at  the  output  ports  are  the 
superposition  of  two  input  signals.  Any  desired  output  signal  amplitude  ratio  can  be 
achieved  by  selecting  the  proper  phase  difference  as  we  now  discuss. 

Bailey  has  shown  that  the  transmission  line  length  <?,  is  related  to  the  output  power 
ratio  by 


(4.4-20) 


The  output  powers  at  each  port  are  related  to  the  fraction  K  by  (4.4-1).  Thus,  the  length 
of  the  line  required  is  expressed  in  terms  of  K  as 


arctan 


1  —  K 
K 


(4.4-21) 


When  a  2/4  quadrature  hybrid  is  used  (/„  =  2/4),  the  equation  simplifies  to 


(4.4-22) 


The  value  for  /c  is  discussed  in  the  next  subsection. 


4.4.3  Verification  of  power  divider  network  by  MCAP 


Before  we  can  proceed  to  synthesis  using  power  divider  networks,  it  is  necessary  to 
confirm  that  the  network  has  acceptable  bandwidth.  Also,  the  three  phase  variables  in 
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(4.4-18),  4>n,  0,j  and  which  depend  on  the  physical  configuration  of  the  network 
must  be  obtained.  We  investigated  these  points  using  MCAP. 

Figures  4.4-3  and  4.4-4  plot  the  magnitude  and  phase  response  of  quadrature 
hybrid  power  dividers  as  a  function  of  frequency  calculated  with  MCAP.  These  power 
dividers  are  designed  for  2733  MHz  center  frequency  using  1/16  inch  substrate  with 
relative  dielectric  constant  et  =  2.5  .  Four  cases  are  considered:  1)  K  =  0.4,  2)  K  =  0.25, 
3)  K^O.l  and  4)  K*0.01.  The  magnitude  of  the  computed  output  power  ratio  plotted 
in  Fig.  4.4-3  shows  that  the  output  remains  close  to  the  design  values  of  K  over  about 
10  %  bandwidth.  One  the  other  hand,  the  phase  plot  in  Fig.  4.4-4  shows  that  for  the 
extreme  case  of  K  =  0.01  the  phase  difference  between  the  two  output  ports  deviates 
from  zero  for  operation  off  of  the  design  frequency.  Thus,  some  caution  must  be  used 
when  a  large  power  ratio  is  required  for  the  compensation  network. 

Next,  the  phase  variables  must  be  determined.  In  the  case  of  the  quadrature  hybrid 
power  divider,  the  length  in  Fig.  4.4-2  is  the  only  variable  undefined.  It  would  be 
reasonable  to  assume  that  the  phase  variables  depend  on  the  line  length.  Indeed,  this 
is  confirmed  by  MCAP  as  we  now  show. 

From  (4.4-18)  the  S-matrix  for  the  variable  power  divider  has  the  following  phase 
components 


011  n  —  011  +  013  ~  023  013 

n  ~  011  +  013  ""  023  011  023 

0 13  023  0 


(4.4-23) 


For  the  case  *  0,  the  phase  of  S-parameters  are  calculated  by  MCAP  as 
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Power  Ratio 


Frequency,  MHz 


Figure  4.4-3.  Magnitude  response  of  the  2-wav  variable  power  divider  in  Fig.  4.4-2 
for  K  =  0.4,  K  =  0.25,  K  =  0.1  and  K-0.01  designed  for  273?  MHz 
as  a  function  of  frequency.  The  substrate  is  1/16  inch  with  c,  -  2.5  . 
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Figure  4.4-4.  Phase  response  of  the  2-way  variable  power  divider  in  Fig.  4.4-2  for 
K-0.4,  K*0.25,  K  =  0,1  and  K*0.0l  designed  for  2733  MHz  as  a 
function  of  frequency.  The  substrate  is  1/16  inch  with  t,  -  2.5  . 
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Equating  (4.4-24)  with  (4.4-23)  it  is  seen  that 


0n  =0 

013  =  n/ 2  for  0SC  =»  0 

023  =  W  2 


For 


- «  0 
0  —  n 
0  0 


0 

0 

0 


and 


0ii  “  -  n 

0,3-0  for  /tfc  =  -2- 

023  ”0 

Comparing  between  (4.4-25)  and  (4.4-27)  it  is  seen  by  inspection  that 

0n  ™  — 

0 13-f-^e 

023  -f-^c 
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Using  the  above  result  in  (4.4-19),  the  S-matrix  for  the  variable  power  divider  network 
can  be  rewritten  as 


Ke“jWt 


-  VK(1  -  K) 


-  VK(1  -  K)  e~mf'  **>  Jl-K  ej( f 
(l-KJe-^*^  VKe^f-^^)) 


(4.4-29) 


We  see  from  this  result  that  /c  acts  to  control  the  phase  at  both  output  ports  1  and  2. 


4.5  Transmission  Line  Network 

The  third  device  considered  for  the  mutual  coupling  compensation  network  was  a 
section  of  transmission  line.  In  this  technique  the  transmission  line  is  placed  between 
the  feed  and  the  array  networks  as  in  attenuator/phase  shifter  compensation.  The  form 
of  the  S-matrix,  however,  is  quite  different  in  the  two  cases.  Unlike  the  attenuator/phase 
shifter  case,  a  section  of  transmission  line  cannot  be  matched  to  the  array  network  and 
the  feed  network  unless  its  characteristic  impedance  is  the  same  as  that  of  the  networks. 
Thus,  for  this  technique  both  the  transmission  coefficient  and  the  reflection  coefficient 
must  be  considered,  even  for  an  ideal  line. 

Figure  4.5-1  shows  a  typical  micro  strip  line  with  characteristic  impedance  Z  and 
length  (.  Smith  [S-1J  has  shown  that  this  general  line  has  the  S-matrix 
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Zo  Z  Zo 


Figure  4.5-1.  Transmission  line  characteristic  impedance  network. 
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SO 


(4.5-1) 


where 


r  “  reflection  coefficient 

jsin/?tf(Z2-Zg)  (4.5.2) 

2ZZq  cos  fit  +  j(Z2  +  Z^)  sin  fit 


and 


T  *  transmission  coefficient 
_ 2ZZg _ 

2ZZq  cos  fit  +  j(Z2  +  Z$)  sin  fit 


(4.5-3) 


Zb  in  the  above  equations  is  the  normalization  impedance  of  the  feed  network  and  the 
antenna  array.  In  microstrip  networks  Zo  =  50£2  is  usually  used. 

There  are  many  empirical  equations  that  relate  the  characteristic  impedance  of  a 
microstrip  line  to  the  physical  dimensions  of  the  line.  The  expressions  given  by  Wheeler 
[W-5]  and  Schneider  [S-l  1]  are 


Zo« 


60 


V®rT 

120* 

V£re 


,  W_\ 

H  W  +  4h  / 

[ -j-  +  1.393  +  0.667  ln^-^-  +  1.444)] 


for  W/h  <;  1 
for  W/h  >  1 


(4.5-4) 


where  W  is  the  microstrip  line  width  and  h  is  the  substrate  height  as  illustrated  on  Fig. 
4.5-2.  The  effective  relative  dielectric  constant  is  given  by  (W-5J 
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£ 


re 


£r  +  1 
2 


+ 


Cf~  1 
2 


lOh 

W 


(4.5-5) 


These  two  expressions  are  reported  to  have  a  relative  error  of  less  than  two  percent 
(G-1J. 

There  are  many  other  two  port  networks  which  can  be  used  in  place  of  the 
transmission  line  as  a  compensation  network.  We  have  considered  four  other  passive 
networks:  series  impedance  network,  shunt  admittance  network,  shunt  connected  open 
ended  stub  network  and  shunt  connected  short  circuited  stub  network.  These  networks 
are  shown  in  Fig.  4.5-3.  The  expressions  for  the  transmission  coefficients  and  the 
reflection  coefficients  of  these  networks  are  given  by  Gupta  [G-IJ.  For  the  series 
impedance  network, 


r51 


Z  +  2Zo 


(4.5-6) 


and 


tSi  2Zq 
“  Z  +  2Zo 


(4.5-7) 


Where  Z  is  the  impedance  of  a  component  connected  in  series.  The  shunt  admittance 
network  has  very  similar  expressions  as  the  above  two: 


fSA 


-  Y 

"  Y  +  2Y0 


(4.5-8) 


J 


J 


J 


ySA  ^ _ 


2Y0 

Y  +  2Y0 


(4.5-9) 
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where  Y  is  the  admittance  of  a  component  connected  in  shunt  with  the  ports.  As  in  the 
case  of  first  two  networks  the  shunt  connected  open  ended  stub  and  the  shunt  connected 
short  circuited  stub  have  similar  expressions  for  the  reflection  coefficients  and  the 
transmission  coefficients.  For  the  former  network  the  coefficients  are  given  by 


pOS 


Zp  tan  /?7 
-  Zq  tan  0S  +  \2Z 


(4.5-10) 


and 


jOS 


12Z 

-Zq  tan  /?/  -I-  j2Z 


(4.5-11) 


where  the  stub  has  the  characteristic  impedance  7  and  the  length  7  .  For  the  latter 
network. 


r58 

ySS 


Zq 

Zq  +  jZ  tan 

j2Z  tan  fit 
Zq  +  j2Z  tan  PS 


(4.5-12) 


(4.5-13) 


where  Z  and  7  are  defined  as  before. 
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V.  Compensation  with  an  Attenuator/Phase  Shifter 
Network 


5.1  Introduction 

In  this  chapter  the  system  of  nonlinear  equations  derived  in  Section  2.4  will  be 
solved  for  a  specific  compensation  network  using  the  S-matrices  obtained  in  the  previous 
chapter.  The  first  type  of  compensation  network  considered  uses  the  attenuator/phase 
shifter  network.  This  technique  was  discussed  extensively  by  Smith  [S-l].  He  has 
successfully  solved  the  equations  using  computer  program  Synthesis  of  Array  with 
Network  and  Element  Coupling  (SANE). 

In  the  first  stage  of  a  companion  experiment  at  New  Mexico  State  University  the 
attenuator/phase  shifter  network  was  used  to  verily  the  mutual  coupling  compensation 
technique.  This  network  was  selected  because  it  is  simple  to  implement  and  the 
S-parameters  of  the  compensation  network  can  be  easily  measured.  The  attenuator 
phase  shifter  network  is  also  the  simplest  to  deal  with  theoretically. 
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We  begin  this  chapter  with  the  derivation  of  nonlinear  equations  specific  to  the 
attenuator/phase  shifter  compensation  network,  followed  by  a  discussion  of  the  physical 
constraints  imposed  on  compensation  network  variables.  Then,  the  relationship  between 
the  generator  phase  and  power  efficiency  of  the  attenuator/phase  shifter  network  is 
shown.  Finally,  the  nonlinear  property  of  attenuator/phase  shifter  compensation 
network  is  demonstrated. 

5.2  Nonlinear  Equations  for  an  Attenuator/Phase  Shifter  Network 


This  section  develops  a  reduced  set  of  nonlinear  equations  for  the  attenuator/phase 
shifter  network.  The  scattering  matrix  for  a  single  attenuator/phase  shifter  combination 
was  derived  in  Section  4.3.  The  S-matrix  of  an  N  element  array,  [SA],  and  the  S-matrix 
of  N-way  feed  network,  [SF],  are  assumed  to  be  known.  Thus,  in  Fig.  2.4-1  and  (2.4-14) 
and  (2.4-15)  we  are  assuming  M  =*  N.  Figure  5.2-1  shows  the  structure  of  the 
compensation  network  for  the  N  element  array.  The  complex  variable  X,  represents  the 
combined  effect  of  an  attenuator  and  a  phase  shifter  between  port  i  and  port  N  +  i. 
When  ideal  devices  are  assumed,  the  magnitude  and  the  phase  of  X,  represent  an 
attenuator  and  a  phase  shifter,  respectively.  The  S-matrix  of  this  network  is  given  as 


[SC] 


[0]  [XN] 
[XN]  [0] 


where 


(5.2-1) 
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V.  Compensation  with 


[XN] 


Xj  0  0 

0  X2  0 

0  0  x3 


0 

0 

0 


0  0 


(5.2-2) 


Substituting  [S0}  into  (2.4-15)  the  equations  reduce  to 

0  =  -U,  +  X,bf 
0  =  -U2  +  X2b£ 

(5.2-3) 

0--UN  +  XNbj5 

Remember  that  {U,}  are  intermediate  unknowns  representing  reflected  waves  from  the 
compensation  network  to  the  feed/generator  network.  It  is  clear  from  the  above 
equations  that  the  variables  {U,}  are  related  to  the  unknowns  {X,}  by 

Uj  =  Xjbf  for  j  =  1  ...  N  (5.2-4) 

Substituting  (5.2-1),  (5.2-2)  and  (5.2-4)  into  (2.4-14)  and  noting  that  S£n+,  =  X,  for  k=  i 
and  zero  for  k  #  i  the  nonlinear  equations  become 

N 

0  -  -aA  4-  X,cf/G  +  ^SyGX,XjbjA  for  i  =  1  ...  N  (5.2-5) 
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In  the  essence,  N  +  M  »  2N  nonlinear  equations  have  been  reduced  to  N  equations  due 
to  specific  structure  of  the  attenuator/phase  shifter  scattering  matrix  [S°].  This  result 
can  now  be  solved  using  numerical  techniques. 

Since  (5.2-5)  imposes  no  constraint  on  the  range  of  X„  it  is  possible  that  a  solution 
represents  a  physically  impractical  network  such  as  an  attenuator  with  a  value  greater 
than  unity.  Such  a  network  represents  gain  that  requires  the  use  of  active  components, 
which  increases  the  complexity  of  compensation  network. 

The  use  of  active  components  can  be  avoided  by  reducing  the  magnitude  of 
attenuators.  It  is  tempting  to  correct  this  situation  by  normalizing  the  variables. 
However,  this  cannot  be  done  because  of  nonlinear  nature  of  the  system  of  equations 
[S-l,  Sec.  6.2].  Smith  has  solved  this  problem  in  SANE  by  varying  the  magnitude  of 
generator  G  until  the  largest  magnitude  of  (XJ  becomes  close  to  unity  within  a  specified 
margin.  This  is  accomplished  by  adding  a  constraint  to  the  nonlinear  equations  which 
can  be  written  as 

1.0-e<max{|Xi|}  £  1.0  (5.2-6) 

where  e  <  <  1.  This  ensures  that  the  variables  {X,}  represent  attenuation.  The  left  hand 
side  of  the  inequality  eliminates  one  attenuator  from  the  compensation  network  by  an 
appropriately  small  c  .  In  other  words,  at  most  N  -  1  attenuators  and  N  phase  shifters 
are  required  to  compensate  for  mutual  coupling  in  an  N  element  array  and  N-way  feed 
network. 
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5.3  Relation  Between  Compensation  Parameters  and  the  Generator  Phase 


The  nonlinear  equations  contain  one  more  variable  which  can  be  varied  to  optimize 
the  compensation  network  for  power  efficiency.  The  variable  is  the  absolute  phase  of 
desired  element  currents.  Equivalently,  the  phase  of  generator  G  can  be  varied  to  give 
the  same  effect.  Unlike  relative  phase  differences  among  element  currents  which  steer 
main  beam,  the  absolute  phase  of  the  currents  have  no  effect  on  the  radiation  pattern. 
However,  the  solution  values  for  {X,}  to  (5.2-5)  change  with  the  absolute  phase,  as  we 
now  illustrate  with  an  example. 

The  change  in  the  solution  vector  with  respect  to  the  absolute  current  phase  can 
be  illustrated  with  a  two  element  array  antenna.  Consider  a  two  element  array  with 
d  »  0.4/l0  element  spacing.  Also  assume  that  the  elements  have  -10  dB  mutual  coupling 
with  -/?d  phase  shift,  -17  dB  mismatch  with  -90°  phase  shift,  and  are  fed  by  an  ideal  tee 
as  shown  in  Fig.  5.3-1.  The  S-matrices  of  the  elements  and  the  feed  network  are 


[sA]- 

[sFj- 


0.00000 -jO.  14142 
-0.25583  -jO.  18587 

-0.50000  0.50000 

0.50000  -0.50000 
0.70711  0.70711 


-0.25583 -jO.  18587 
0.00000  -jO.  14142 

0.70711  ‘ 

0.70711 

0.00000 


(5.3-1) 


(5.3-2) 


Suppose  the  desired  currents  for  the  array  are 

I,  «  1 .0/0.0 
I, ,  l.Q/  -9Q,o 


(5.3-3) 
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This  excitation  gives  a  sum  pattern  with  main  beam  steered  to  51.3°  off  broadside. 

The  required  compensation  network  variables  were  solved  using  SANE  for  the 
generator  phase  between  0°  and  180°  for  every  5°.  Figure  5.3-2  shows  the  magnitude 
and  phase  of  variables  {XJ  with  respect  to  the  generator  phase.  The  plots  clearly  show 
the  variations  of  {XJ  with  generator  phase.  Figure  5.3-3  shows  the  generator  magnitude 
required  to  excite  the  elements  at  the  desired  currents  while  keeping  the  peak  I X,  I  at 
unity;  this  generator  magnitude  with  the  generator  phase  of  Fig.  5.3-2  produces  the 
required  element  currents  (5.3-3).  The  generator  magnitude  can  be  interpreted  as  being 
a  measure  of  the  efficiency  of  compensation  network,  the  lowest  generator  magnitude 
being  most  desireable.  In  the  array  considered,  the  compensation  network  is  most 
efficient  at  /Gad 20°  because  it  requires  the  minimum  generator  magnitude.  On  the 
other  hand,  the  compensation  network  is  least  efficient  at  /G=s45#.  Comparison  of  Figs. 
5.3-2  and  5.3-3  shows  that  the  lowest  generator  value  and  the  lowest  attenuator  setting 
occur  together  (at  the  same  /G).  This  property  is  seen  only  for  certain  2-element  arrays 
and  cannot  be  generalized  for  N-element  arrays. 

The  calculated  difference  in  gain  from  the  most  efficient  to  the  least  efficient 
compensation  networks  is  1.9  dB.  The  difference  is  small  in  this  case.  However,  for  an 
array  with  large  number  of  elements  the  difference  can  be  significant.  For  example,  an 
eight  element  35  dB  Dolph-Chebyshev  array  used  in  the  experiment  (described  in 
Chapter  8)  has  4.5  dB  difference  in  gain  from  the  least  efficient  compensation  network 
to  the  most  efficient. 
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.3-2.  The  magnitude  and  the  phase  of  compensation  network  value  with 
respect  to  the  absolute  phase  of  generator  for  two  clement  array 
example  of  Fig.  5.3-1. 


5.4  The  Nonlinear  Property  of  Attenuator/Phase  Shifter  Compensation 

In  attenuator/phase  shifter  compensation,  the  value  of  compensation  variables 
required  for  a  given  antenna  array  are  unpredictable  and  are  often  surprising.  This 
phenomenon  was  shown  with  a  simple  numerical  example  by  Smith  [S-l,  Sec.  6. 2]  and  is 
demonstrated  here  with  a  similar  example.  Consider  the  array  used  in  the  previous 
section  (shown  in  Fig.  5.3-1)  which  is  represented  by  the  S-matrices  (5.3-1)  and  (5.3-2). 
Analysis  of  the  array/feed  combination  using  MCAP  gives  the  following  element 
currents: 

I,  =0.91767  /  14.61° 

(5.4-1) 

I2  =  0.91767  /  14.61° 

These  excitations  have  equal  magnitude  and  equal  phase  because  both  the  feed  and  the 
element  coupling  are  symmetric.  Now  suppose  the  desired  currents  for  the  array  are 

1,  =«  1.00000  /  0.0° 

(5.4-2) 

I2~  1.00000  /  -  90.0° 

which  are  the  same  as  in  the  previous  section.  This  excitation  is  not  symmetric  and  if 
realized  would  give  a  main  beam  at  0O  =  51.318°  from  broadside. 

Intuitively,  one  would  think  that  the  compensation  network  requires  only  —90° 
phase  shift  for  element  2  since  the  magnitudes  of  the  uncompensated  currents  are 
identical  for  both  elements.  According  to  the  plots  in  Fig.  5.3-2  calculated  using  SANE, 
however,  an  attenuator  and  a  phase  shifter  are  required  for  the  element  1  whereas  only 
a  phase  shifter  is  required  for  the  element  2.  In  addition,  the  difference  in  the  required 
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phase  shifts  between  element  1  and  element  2  is  not  90°  .  For  example,  the  required 
compensation  variables  for  /G  =  0°  obtained  using  SANE  are 


X,- 0.76154  /  9.80° 
X2-  1.00000  /  258.32° 


(5.4-3) 


Clearly,  an  attenuator  is  required  only  for  element  1.  It  is  also  apparent  that  the 
difference  in  phase  shifts  is  not  90.0° . 

The  above  example  demonstrates  how  the  unexpected  can  happen  even  for  simple 
arrays.  Since  this  property  is  not  directly  apparent,  we  present  a  treatment  of  how  the 
compensation  network  affects  the  currents  in  the  antenna  element  in  a  step  by  step 
fashion  for  the  various  coupling  mechanism. 

In  general,  the  S-matrices  of  two  element  array,  feed  and  attenuator/phase  shifter 
compensation  networks  (see  Fig.  5.4-1)  can  be  represented  as 


[sc]  = 


SA 

311 

C/1 

z> 
_ 1 

SA 

321 

CA 

SF 

SF 

a12 

SF 

a13 

SF 

SF 

322 

CF 

•*23 

SF 
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SF 

032 

SF 

s33 

0 

0  X, 

°1 

0 

0 

0  X2 

X, 

0 

0 

0 

0 

X2 

0 

0 

(5.4-4) 


(5.4-5) 


(5.4-6) 
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Figure  5.4-1.  Current  components  in  a  two  element  array  with  matched  array  and 
feed  network. 
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5.4.1  Matched  Case 


The  nonlinear  effect  can  be  demonstrated  without  including  mismatches.  In  fact, 
even  for  the  two  element  case  tracing  individual  signals  becomes  very  involved  with 
reflections  present,  so  we  neglect  them  for  this  discussion.  The  compensation  network 
is  assumed  to  be  matched  to  the  antenna  and  to  the  feed  at  the  respective  ports.  The 
generator  is  also  assumed  to  be  matched  to  the  feed.  Thus,  Sfi  =  S&  =  0  and 
Sft  <■  S&  —  S£,  =  0.  It  is  physically  impossible  for  multiport  networks  other  than  two  port 
networks  to  have  all  ports  matched.  However,  this  constraint  is  neglected  in  the  feed 
network  assumptions  to  simplify  the  analysis.  Under  these  assumptions,  there  will  be  no 
reflected  waves  generated  at  the  junctions  between  antenna  elements  and  feed,  and 
between  feed  and  generator. 

Using  the  above  S-parameters,  the  current  in  antenna  element  1  can  be  represented 
as  an  infinite  sum  of  circulating  currents  that  are  shown  in  Fig.  5.4-1: 

i,  -  x,s3f,g  -  s£x2s3f2g  +  x1s2F,x2sf2x1s3F,G  -  s£x2sf2x,s£x2s3F2G  +  ••• 

oo 

-  £(X,S,f,XjS^)"  x,s,f,g  -  (s{iX,Sf2X,)"  S?,X,SF  G  <5  ) 

n«0 


Assume  that  the  antenna  and  the  feed  are  reciprocal  networks. 


(5.4-8) 


Then,  (5.4-7)  can  be  rewritten  as 
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oo 

I,  -  ^(X.SfjXjS*,)" X,S?,G  -  (X,sfjX,S*)"  sfjXjSfjG 

n-0 

oo 

-  (x,s3F,G  -  s*XjSfg) 

n— 0 

oo 

=  (X,sf,  -  XjSfjSfjjG^fX.XjSfjSf;) 


Showing  the  first  term  of  the  summation  explicitly, 


I,  -  (X,SjF,  -  XjSfjSjFj)G 


1  + 


OO 

£(X,X2SFSA)” 

n-l 


OO 

x,sf,g  -  XjS^G  +  (x,sF  -  X2s*  sf:)g^](x,xisf  s*)' 

n—1 


The  current  at  element  2  can  also  be  expressed  in  similar  way. 


OO 

XjSl’jG  -  X,SjA,SjF,G  +  (X2SFj  -  X,SjAiS,F,)G^(X,X2srjSl'J)' 


n— 1 


These  current  expressions  consist  of  three  terms: 


1  "  I  ideal  I  antenna  coupling  "b  ^ full  coupling 


where 


^1  ideal  **  XiS31G 


^1  antenna  coupling  —  ^2^12^32^ 


(5.4-9) 


(5.4-10) 


(5.4-11) 

(5.4-12) 

(5.4-13) 

(5.4-14) 


V.  Compensation  with  an  Attenustor/Phsse  Shifter  Network 


100 


(5.4-15) 


I,  -  (X.Sf,  -  XjSfjSfjjG^fX.XjSfjSfj)" 

n»l 

for  the  current  in  element  1. 

The  first  term  IMt„  is  the  current  (in  the  case  of  element  1)  that  would  be  present 
in  element  1  if  no  mutual  coupling  were  present  anywhere  in  the  system.  This  current 
is  directly  proportional  to  the  coupling  from  input  to  output  of  the  feed  network 
(represented  by  Sf,  for  element  1  and  S&  for  the  element  2),  and  to  the  magnitude  of 
generator,  G.  This  ideal  current  can  be  changed  by  the  attenuator  and  phase  shifter 
values  in  a  linear  fashion  (i.e.  1  dB  increase  in  attenuation  reduces  the  output  currents 
by  1  dB). 

1 _ T...t  is  the  contribution  to  the  current  (in  the  case  of  element  1)  from  antenna 

element  2  through  antenna  coupling.  It  is  the  current  contribution  when  there  is  mutual 
coupling  between  antenna  elements  but  no  coupling  exists  in  the  feed.  This  is  similar  to 
the  case  of  independent  generator  excitation  as  discussed  in  Sec.  2.3.1.  The  actual 
element  current  in  such  a  system  is  the  sum  of  I*,*  and  -I,nltnni00uphn,.  I.™™. is  a 
negative  current  because  it  circulates  in  the  opposite  direction  from  Iide„.  In  this  case, 
the  equations  (5.4-13)  and  (5.4-14)  show  a  linear  dependence  on  unknown  {XJ.  Thus, 
the  compensation  network  can  be  calculated  by  simply  solving  a  system  of  linear 
equations. 

For  the  third  component,  Inm  coupe.**  to  be  nonzero  the  two  output  ports  of  the  feed 
network  must  not  be  isolated  from  each  other;  i.e.,  Sf2  =  SJ,  #  0.  The  coupling  between 
these  two  ports  introduces  a  feedback  loop  in  the  system  which  makes  the  current  a 
nonlinear  function  of  compensation  variables.  This  is  seen  in  (5.4-15)  by  the  presence 
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of  the  products  X,X,.  Obviously  this  term  is  not  present  if  either  antenna  coupling, Sf^ 
or  feed  coupling, Sf^  is  zero. 

5.4.2  Comment  on  generator  excitation  dependence 

It  should  be  noted  that  the  required  attenuator  variation  with  respect  to  the 
generator  phase  as  discussed  in  the  previous  section  does  not  occur  in  this  particular 
array.  For  example,  suppose  the  two  element  array  in  Fig.  5.3-1  has  matched  feed  and 
element  ports.  The  desired  element  currents  are  those  of  (5.3-1).  The  requireed 
compensation  variables  as  with  respect  to  generator  phase  calculated  from  SANE  are 
shown  in  Fig.  5.4-2.  Clearly,  the  required  attenuators  values  remain  constant  as  the 
generator  phase  is  varied.  The  required  phase  shifters,  on  the  other  hand,  show  the 
dependence  to  the  generator  phase  similar  to  the  case  in  the  previous  section.  The 
efficiency  of  the  compensation  network  also  varies  with  respect  the  generator  phase.  This 
is  demonstrated  in  a  plot  of  the  required  generator  magnitude  versas  the  generator  phase 
shown  in  Fig.  5.4-3. 

This  phenomena,  however,  cannot  be  generalized  to  array  with  larger  number  of 
elements.  In  the  case  of  a  N-element  array,  the  generator  phase  effect  on  the  attenuator 
values  can  be  observed  even  when  all  the  elements  and  feed  ports  are  matched. 

5.4.3  Mismatched  case 

The  preceding  current  analysis  assumed  that  all  network  components  are  matched. 
When  mismatches  are  introduced,  the  incident  waves  at  each  junction  will  be  split  to  two 
components;  reflected  and  transmitted  waves.  Thus,  the  current  expression  cannot  be 
written  in  simple  summation  form.  However,  it  can  be  generalized  from  the  matched  case 
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Phase  of  Generator,  degrees 


Figure  5.4-3.  The  required  generator  magnitude  with  respect  to  the  absolute  phase 
of  the  generator  for  the  two  element  array  with  matched  element  and 
feed  ports. 
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that  even  in  the  unmatched  case  nonlinearly  will  occur  when  feed  network  coupling  is 
present. 
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VI.  Synthesis  with  a  Power  Divider  Network 


6.1  Introduction 

The  second  type  of  compensation  technique  considered  uses  a  variable  power 
divider  network.  Smith  [S-l]  has  originally  treated  this  problem  using  the  general 
nonlinear  equations  of  (2.4-14)  and  (2.4-15).  These  equations  were  solved  by  Smith 
using  computer  program  SANE-GECON,  Synthesis  of  Array  with  Network  and  Element 
Coupling  for  General  Compensation  Network  [S-IJ.  Though  the  solution  vectors  were 
found  for  different  cases  of  mutual  coupling  and  desired  element  excitation,  the  program 
converged  with  difficulty.  Convergence  problems  are  common  in  the  solution  of 
nonlinear  equations.  As  we  found  with  the  attenuator/phase  shifter  compensation 
network  in  Chapter  5,  a  general  solution  technique  cannot  be  applied  directly  to  solve 
the  nonlinear  equations.  We  must  include  the  specific  nature  and  the  physical 
constraints  of  the  compensation  network  before  solution  is  attempted. 

In  this  chapter  we  first  present  three  preparatory  steps  necessary  to  solve  for  the 
power  divider  compensation  network  using  numerical  techniques.  First,  the 
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S-parameters  of  an  N-way  power  divider  are  derived.  Then  the  nonlinear  equations  are 
simplified  for  the  power  divider  network.  Finally,  the  initial  guess  for  the  solution  vector 
which  plays  important  role  in  the  numerical  technique  is  discussed. 

The  steps  described  in  this  chapter  were  used  to  develop  the  computer  program 
SANE-PODCON,  Synthesis  of  Array  with  Network  and  Element  Coupling  for  Power 
Divider  Compensation  Network;  see  Appendix  B  for  the  details.  This  program  converges 
much  faster  than  SANE-GECON  for  various  cases  of  mutual  coupling. 


6.2  N-way  Power  Divider 

Before  we  can  proceed  with  the  synthesis  of  the  power  divider  network,  it  is 
necessary  to  look  at  the  implementation  of  an  N-way  power  divider.  We  choose  a  power 
divider  built  up  from  2-way  power  dividers.  This  is  done  because  only  2-way  variable 
power  dividers  are  available  (see  Sec.  4.4.2).  In  this  section,  the  S-matrix  of  a  4-way 
power  divider  is  derived  from  known  results  for  2-way  power  dividers.  Then,  the  result 
is  generalized  to  the  S-matrix  for  N-way  power  dividers. 

6.2.1  S-matrix  representation  of  a  4- way  power  divider 

The  4-way  power  divider  is  modeled  from  three  2-way  power  dividers  as  shown  on 
Fig.  6.2-1.  S-matrix  representation  of  a  4- way  power  divider  can  be  obtained  from  the 
S-matrices  of  2-way  power  dividers  using  the  matrix  connection  technique  discussed  in 
Section  2.1.  The  2-way  power  dividers  are  represented  by  [S®1-]  ,  [S35*]  and  [S0>c]  for 
the  left,  the  right  and  the  center  power  divider,  respectively.  The  superscript  (2)  in  the 
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notation  indicates  a  2-way  power  divider.  These  matrices  are  of  the  form  (4.4-18)  and 
they  can  be  realized  as  discussed  in  Section.  4.4.2. 

The  combined  S-matrix  for  the  4-way  power  divider  is 


[S®L]  0  0 

0  [s0^]  0 
0  0  [s^] 


(6.2-1) 


Rearranging  the  S-parameters  for  the  internal  and  external  ports  gives 


[*"] 


“ 
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s?>Ls£)L  0  0  0 

sgL  0  0  0 

sfftsgL  0  0  0 

o 

o 

o 

1? 

0  o  sf,)Rsf2)R  0 

0  Sj3R  0  0 

0  o  s$?RsgR  0 

0  sgR  0  0 

0  0  0  0  sg* 

0  o  sSfsSf 

S?i)LS?2L  ooo 

sgL  ooo 

0  o  s?>RsgR  o 

0  sgR  0  0 

0  0  0  0  sg® 

0  o  sf^s^f 

0  0  0  0  Sy0 

1 

0  0  s^s^ 

_ 

(6.2-2) 


Using  the  connection  technique  discussed  in  Section  2.1,  the  reduced  S-matrix  for  the 
4-way  power  divider  can  be  calculated  from 


[s()]-  [sPp  +  Spc(r  —  scc)  'Sep] 


(6.2-3) 
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where  the  connection  matrix  is 


r 


0  0  10 
0  0  0  1 
10  0  0 
0  10  0 


(6.2-4) 


The  port  numbers  corresponding  to  row/column  numbers  in  F  are  the  same  as  the 
definition  used  in  the  submatrix  in  (6.2-2).  For  example,  row/column  number  1 
corresponds  to  port  3  of  the  left  2-way  power  divider.  Note  from  (4.4-18)  that  a  2-way 
power  divider  has  a  matched  input  port,  so  that  SgL  *  SgR  «  Sg0  *  0  .  Evaluation  of  the 
right  hand  side  of  (6.2-3)  proceeds  as  follows: 


(r  -  scc) 


(r  ~  scc) 


-i 


0  0  1 
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1  o-s?* 
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c(2)C  «j<2)C 
321  a22 

0  1 

1  0 

0  0 

0  1 

0  0 
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Spc(r  —  scc) 
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(6.2-5) 


The  matrix  equation  (6.2-5)  represents  S-parameters  of  a  general  4-way  power  divider 
network.  This  matrix  is  used  as  a  building  block  for  the  S-matrix  of  an  N-way  power 
divider  network. 


6.2.2  Derivation  of  the  N-way  power  divider  S-matrix 

The  results  of  previous  section  can  be  extended  to  that  for  an  N-way  power  divider. 
In  this  and  subsequent  sections  the  N-way  power  divider  will  be  represented  by  [S™]. 
The  network  consists  of  two  N/2  —  M  power  dividers  and  a  2-way  power  divider  as 
shown  on  Fig.  6.2-2.  Each  M-way  power  divider  is,  in  turn,  made  of  2-way  power 
dividers  also.  We  must  have  then  N  —  1  2-way  dividers  to  realize  an  N-way  power 
divider.  The  matrix  [S^]  can  be  written  by  generalization  of  (6.2-5)  as 
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Figure  6.2-2.  Implementation  of  an  N-way  power  divider  network  from  two 
M-way  power  divider  networks  and  a  2-way  power  divider  network. 
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[«P1 


C(M)L  c(2)Cc(M)R 
al,M+laU  aM+l,M 

e(M)L  <;(2)C<.<M)R 
a2,M+JaU  aM+l,M 


c(M)L  <j(2)C«(M)R  o(M)L  ~(2)C<j<M)R 
aM,M+lal,2  aM+l.l  aM,M+laU  aM+lJ 


c(M)L  <.(2)C<,<M)R 
aM,M+laU  aM+l,M 
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o(M)R  o(2)Cq(M)L 
3I,M+1°2,1  aM+l,l 

e<M)R  «;<2)C~<M)L 
a2,M+la2,l  aM+l,l 


o(M)R  q(2)Cq(M)L 
al,M+l32,l  ^M+U 

q(M)R  c(2)Cc(M)L 
a2,M+la2,l  aM+U 


c(M)R  c(2)Cc(M)L 

31. M+la2,l  3M+l,M 

C(M)R  q(2)Cc(M)L 

32, M+la2,l  aM+l,M 


c(M)R  o(2)Cc(M)L  o(M)R  c(2>Cc(M)L 
aM,M+la2,!  aM+l,l  aM,M+la2,I  SM+U 


e(M)R  o(2)Co(M)L 
^M.M+1^2,1  aM+l,M 


c(M)R 


+  S 


(M)R  o(2)C 
l,M+la2,2 


s 


(M)R 

M+1,1 


c(M)R 


+  s 


(M)R 

l.M+1 


c(^)Cp(M)R 

a2,2  3M+1,2  •" 


c(M)R  ,  p(M)R  p(2)Cp(\l)R 
a2,l  aM+l,l 


c(M)R 

a2,2 


a2,M+l 


s 


(2)Cp(M)R 
2,2  aM+U  — 


[sft’]  - 


•<M)R 
'M, 


+  S 


(M)R  c(2)Cp(M)R 
M.M+1^2.2  3M+l.i 


S'v 


(M)R 


M.2 


+  s 


(M)R 

M.M+1 


s 


<2)Ce(M)R 
2,2  3M+U  — 


c^M)R  ,  o(M)R 

S1,M  +  al,M+l32,2  3M+1„M 

(M)R  +  S(M)R  q(2)Cq(M)R 


J2,M 


2..M+ 1^2,2  aM+l,M 


c(M)R  <;(M)R  p(2)Cc(M)R 
—  aM,M  3M,M+la24  aM+l,M 
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L 


■ 


o(M)L 

c(2)C 

al,M+l 

au 

o(M)L 

<j(2)C 

a2,M+l 

c(M)L 

c(2)C 

aM,M+l 

laU 

c(M)R 

c(2)C 

al,M+l 

a2,3 

c(M)R 

c(2)C 

a2,M+l 

a2,3 

o(M)R 

aM,M+! 

c(2)C 

la2,3 

t: 


c(N)l  _  re(M)L  -<2)0 
aVI  J  —  |_aM+l.la3,l 


e(M)L  0(2)0 
aM+l,2a3,t 


o(M)L  o(2)C 
aM+l,M  a3,l 


<j(M)R  o(2)C  q(M)R  o(2)C 
aM+l,la3^  aM+1.2a3,2 


o(M)R  c(2)C 
aM+l,Ma3,2 


] 


A  complete  S-matrix  for  an  N-way  power  divider  in  terms  of  the  S-parameters  of  2-way 
power  dividers  is  obtained  by  solving  the  above  equations  iteratively.  This  is  done  by 
letting  M  for  iteration  i  be  N  for  the  succeeding  iteration  i+  1,  and  solving  for  the  new 
matrix  using  (6.2-6).  This  process  is  continued  until  M  =  2;  in  other  words,  until  the 
entire  S-matrix  is  represented  in  terms  of  2-way  power  divider  S-parameters. 


Example:  An  8-way  power  divider 

As  an  example  of  the  use  of  (6.2-6),  let  us  solve  for  S$  ,  the  mismatch  at  port  4  of 
an  8-way  power  divider.  The  initial  network  consists  of  two  4-way  power  dividers  and  a 
2-way  power  divider  as  shown  on  Fig.  6.2-3.  First,  we  let  N  =  8  and  M  =  N/2  =  4  in 
(6.2-6).  S^l  corresponds  to  the  bottom  right  element  in  [Sf0]: 


c<«)  _  cWL  cWLc(2)Cc(4)L 
a4,4  a4,4  +  a4,5  al,l  a5,4 
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i- 


The  4-way  power  divider  on  the  left  is  constructed  from  three  2-way  power  dividers  as 
in  Section  6.2.1  (LL,  LC  and  LR  in  Fig.  6.2-3).  The  S- parameters  for  the  4- way  power 
divider  in  (6.2-7)  are  obtained  from  (6.2-5)  or  equivalently  from  (6.2-6)  for  N  =  4  and 
M-2. 


<;(4)L  <;(2)LR  <j(2)LR<j(2)LC<;(2)LR 

a4,4  “  *22  +  a2,3  *22  *32 


c(4)L  -<2)LR~(2)LC 

•*4,5  >22 


o(4)L  <.(2)LR~(2)LC 

35,4  ™  a3 2 


(6.2-8) 

(6.2-9) 

(6.2-10) 


To  illustrate  the  notation,  Sg’^  is  used  to  represent  S-parameter  for  the  right  2-way 
power  divider  which  is  contained  in  the  left  hand  side  4-way  power  divider.  Substituting 
(6.2-8)  to  (6.2-10)  into  (6.2-7), 

«<*>  _  C<2)LR  ,  o<2)LR«j<2)LC<,<2)LR  ,  ~C>Ut~G)LC~(2)CC«j<2)LR<5<2)LC  (f>  ,  .  .. 

a4.4”a2.2  +  324  *22  53J  +  “2,3  “2,3  “1.1  “3,2  “3.2  iO.Z-llj 


where  SfJ0  *=  Sg*0  which  is  apparent  from  Fig.  6.2-3.  Using  the  reciprocity  property,  the 
above  expression  simplifies  to 

S'*i  -  SSLR  +  {sSH'sfr  +  {S™LRS™LC)1S^>CC  (6.M  2) 


Finally,  the  S-parameters  in  (6.2-12)  can  be  written  in  terms  of  the  power  ratio  and 
phase  shift  at  each  power  divider  network  using  (4.4-19). 


SgLR  »(1- k<2)LR)cK<*u  ^ 


(6.2-13) 


SC2)LR  =  ^/k(2)LR  -ft/*1*) 


(6.2-14) 


*  (1  -  Ka)LC)ej(^LC-2^LC) 


(6.2-15) 
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S£)LC  _  ^yK(2)LC  ^2*°  -/V®*0) 


(6.2-16) 


sri>cc  , 

I 


(6.2-17) 


Substituting  (6.2-13)  through  (6.2-17)  into  (6.2-12)  gives 
Sg  -  (1  -  K<2>LR)eK,*M  LR  -V*®LR) 

+  _  K(2)LC^ej(^^LC  -2^LC) 

+  ^ yK( 2)LR  eK^2LR  -V**) yK(2)LC  gjWg1-0  2K(2)CCej^ 


The  procedure  similar  to  this  example  is  implemented  in  SANE-PODCON  to  solve  for 
S-matrix  of  an  N-way  power  divider  network. 


6.3  The  System  of  Nonlinear  Equations 

The  S-matrix  developed  in  the  last  section  can  be  substituted  to  the  nonlinear 
equations  of  (2.4-14)  and  (2.4-15)  to  solve  for  the  compensation  network  using  power 
divider  networks.  The  important  practical  feature  of  this  design  is  that  compensation 
network  is  contained  within  the  feed  network.  By  incorporating  compensation  at  the 
design  stage  no  new  elements  or  layouts  are  needed.  The  difference  in  the  network  with 
compensation  for  coupling  included  is  that  the  individual  2-way  power  dividers  have 
different  power  split  ratio  and  phase  shift  values  from  the  conventional  3  dB  equal  phase 
approach. 
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For  the  power  divider  implementation  the  feed  network  as  defined  in  Fig.  2.4-1 
reduces  to  a  transmission  line  that  connects  to  the  network  C  as  shown  in  Fig.  6.3-1. 
The  S-matrix  for  the  feed  network  is  then 

[sF]-[!  o]  (6-3-,) 


It  follows  that  M  ™  1  for  this  compensation  technique.  Thus,  for  power  divider 
compensation  the  network  C  in  Figs.  2.4-1  and  6.3-1  containes  the  conventional  feed 
junctions  with  compensation  imbedded.  The  network  is  excited  by  a  generator  G  with 
reflection  coefficient  S°  where  wave  induced  by  the  generator,  c°  ,  is  treated  as  an 
unknown  variable.  Using  (2.1-9)  and  (2.1-20),  the  feed/ generator  network  is  described 
by 


SF/G-SG 


(6.3-2) 


„F/G 


(6.3-3) 


There  is  one  intermediate  variable  in  the  nonlinear  equation  since  M  =  1.  Then,  (2.4-15) 
becomes 


o  -  -U,  +  SS-M.,bf  +  •"  +  sS+,.Nt>a  +  sS+i.n+,[cf'°  +  SF'°U,]  (6.3-4) 

The  S-matrix  for  the  power  divider  network  derived  in  the  previous  section  assumes  that 
the  network  has  a  matched  input  port  (SS+liN+i  =  0).  Then,  (6.3-4)  simplifies  to  a  linear 
equation  of  the  form 


U, 


’N+I,lul 


X  +  sc 


N+U^*  + 


+  S1 


N+l.N 


N 

■‘-I- 


'N+lj 


ib, 


(6.3-5) 


i-1 
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Substituting  (6.3-2),  (6.3-3)  and  (6.3-5)  into  (2.4-14)  results 


for  i  —  1 ...  N  (6.3-6) 


Thus,  N  +  M  —  N  +  1  nonlinear  equations  are  reduced  to  N  equations.  However, 
(6.3-6)  is  not  as  simple  as  the  attenuator/phase  shifter  equations  because  the 
S-parameters  of  the  compensation  network,  S£  ,  are  related  to  physical  parameters  in 
more  complicated  fashion  as  was  illustrated  in  (6.2-18)  for  an  8-element  array.  The 
nonlinear  property  in  (6.3-6)  arises  from  the  products  of  unknowns  S^N+1S^,ik.  In 
addition,  for  power  divider  compensation,  the  generator  value  G  is  treated  as  an 
unknown.  Thus,  the  term  S^+,c°  also  contributes  to  the  nonlinear  property  of  the 
equation.  The  nonlinear  equations  are  solved  numerically  by  the  program 
SANE-PODCON  using  the  Damped  Newton's  method  [D-4].  A  detailed  description  of 
SANE-PODCON  is  given  in  Appendix  B. 


6.4  Initial  Guess  to  the  Solution  Vector 

In  the  modified  Newton's  method  (which  is  the  numerical  technique  used  to  solve 
for  the  system  of  nonlinear  equations)  the  initial  guess  for  the  solution  vector  is  one  of 
the  most  important  parameters  the  user  must  supply.  Generally,  a  good  choice  of  initial 
guess  greatly  increases  the  probability  of  convergence  to  a  solution.  It  has  also  been 
proven  that  Newton's  method  converges  quadratically,  so  an  initial  guess  which  is  close 
to  the  solution  also  increases  the  convergence  speed  [D-4,J-4,J-5].  On  the  other  hand, 
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a  poor  initial  guess  could,  at  worst,  lead  to  divergence  from  the  actual  solution.  It  is 
obvious  that  the  initial  guess  cannot  be  left  to  chance. 

There  are  two  requirements  that  the  initial  guess  must  satisfy  in  addition  to  being 
relatively  close  to  the  actual  solution.  First,  the  initial  power  split  ratio  K  at  each  power 
divider  must  be  between  zero  and  one.  This  is  necessary  so  that  the  power  divider  would 
remain  a  passive  device.  The  other  requirement  is  the  unitary  property  of  a  S-matrix, 
which  is  required  for  the  network  to  be  lossless.  Thus,  the  initial  guess  must  be  selected 
such  that  the  matrix  [S0]  satisfies  these  requirements. 

A  possible  choice  for  the  initial  guess  is  to  set  the  variables  such  that  the  power 
divider  network  would  produce  the  required  element  currents  when  there  is  no  mutual 
coupling  in  the  antenna  elements.  Clearly,  this  choice  meets  the  requirements  discussed 
above.  Also,  in  the  case  where  mutual  coupling  among  antenna  elements  is  small,  the 
initial  network  will  be  close  to  the  required  compensation  network.  In  the  small 
coupling  environment  the  effect  on  the  element  currents  would  also  be  small.  Thus,  the 
required  compensation  network  should  not  be  very  different  from  the  initial  network. 

The  SANE-PODCON  program,  which  was  developed  at  Va  Tech  and  is  described 
in  Appendix  B,  is  implemented  with  a  routine  to  give  an  initial  guess  using  the  above 
procedure.  In  all  cases  tested  the  program  converged  to  a  solution.  The  program  has 
even  converged  for  the  experimental  array  described  in  Chapter  8  which  employed  tight 
element  spacings  to  intentionally  increase  mutual  coupling  between  elements.  In  most 
antenna  arrays  mutual  coupling  would  be  less  than  that  of  the  experimental  array.  This 
initial  guess  procedure,  therefore,  should  be  adequate  for  most  antenna  arrays. 
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VII.  Synthesis  with  a  Transmission  Line  Network 


7.1  Introduction 

The  third  technique  considered  uses  the  characteristic  impedance  of  transmission 
lines  to  compensated  for  mutual  coupling.  Smith  treated  this  problem  by  solving  the 
nonlinear  equations  (2.4-14)  and  (2.4-15)  using  SANE-GECON  [S-lj.  The  program 
converged  to  a  solution  for  two  element  array  cases.  However,  convergence  was  never 
achieved  for  higher  element  number. 

Initially,  we  anticipated  that  the  convergence  problem  arose  from  general  treatment 
of  the  nonlinear  equations  by  SANE-GECON  since  similar  problems  were  found  for  the 
power  divider  compensation  network.  Thus,  we  modified  the  nonlinear  equations 
specifically  for  the  transmission  line  network  (this  is  discussed  in  next  section).  Then  a 
computer  program  SANE-TRACON  (Synthesis  of  Array  with  Network  and  Element 
coupling  for  TRansmission  line  Compensation  Network)  was  developed  using  the 
modified  nonlinear  equations.  However,  we  could  not  achieve  convergence  for  the  cases 
other  than  two  element  arrays.  Although  there  is  no  mathematical  proof,  we  speculate 
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that  the  lack  of  convergence  is  due  to  a  physical  limitation  of  the  network  rather  than 
the  numerical  technique  used  to  solve  the  equations.  This  remains  an  open  question. 
Since  the  approach  used  is  simple  and  implementation  is  straight  forward,  future 
attempts  to  find  a  solution  are  warranted.  Also,  an  attempt  to  prove  in  general  a 
necessary  condition  for  a  feed  network  to  yield  a  solvable  compensation  network  would 
be  valuable. 


7.2  The  System  of  Nonlinear  Equations 


The  S-matrices  of  a  single  transmission  line  network  (discussed  in  Section  4.5)  can 
be  substituted  into  the  system  of  nonliner  equations  (2.4-14)  and  (2.4-15)  to  develop  a 
simplified  set  of  nonlinear  equations.  As  in  the  case  of  the  attenuator/phase  shifter 
network,  this  is  a  compensation  technique  using  an  added  network  between  the  element 
and  feed  network.  The  modified  nonlinear  equations,  however,  are  not  as  simple  due  to 
inherent  existence  of  reflection  coefficients  at  all  ports  of  the  compensation  network. 

The  S-matrix  of  an  N-element  array,  [SA],  and  the  S-matrix  of  N- way  feed  network, 
[SF],  are  assumed  to  be  known.  Thus,  we  are  assuming  N»  M  in  Fig.  2.4-1  and  (2.4-14) 
and  (2.4-15).  The  compensation  network  consists  of  N  transmission  lines  as  shown  in 
Fig.  7.2-1.  The  magnitude  and  phase  of  the  complex  variable  X,  represent  characteristic 
impedance  and  length  of  the  transmission  line  between  port  i  and  i+N,  respectively. 
The  S-matrix  of  this  network  is  given  by 


[SC] 


cn  m 
.  m  cn 


(7.2-1) 


where 
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for  i  *  It  is  clear  that  these  equations  are  not  as  simple  as  other  two  cases.  In 

particular,  there  are  N  intermediate  unknowns,  {U1}  and  N  unknown  variables  in  (7.2-2) 
and  (7.2-3).  Hence,  it  is  necessary  to  solve  for  2N  unknowns  in  order  to  determine 
compensation  network  values. 
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VIII.  Microstrip  Array  Antenna  Experiments 


8.1  Description  of  Experimental  Hardware 

The  microstrip  array  used  to  investigate  the  mutual  coupling  compensation  is 
composed  of  two  sections:  a  board  containing  the  array  of  radiating  elements  (the  array 
board)  and  a  feed  network  board.  The  hardware  which  was  constructed  and  tested  at 
New  Mexico  State  University  Physical  Science  Laboratory  was  operated  at  a  center 
frequency  of  2733  MHz  [J-9J.  Both  the  array  board  and  feed  network  board  were 
fabricated  on  1/16  inch  thick  teflon  fiberglass  substrate  with  relative  dielectric  constant 
er »  2.55  . 

The  array  board  consists  of  eight  quarter-wavelength  wide  patches  with  3.43  cm 
(0.312  2„)  interelement  spacings.  At  the  center  frequency  of  2733  MHz  the  wavelength 
is  =*  10.98  cm  in  air  and  2d  —  =  7.706  cm  in  the  substrate  where  the  effective 

dielectric  constant  is  £,*  =  2.03.  The  geometry  of  the  array  board  is  shown  in  Fig.  8.1-1. 
The  quarter-wavelength  patches  are  used  to  obtain  closer  than  one  half-wavelength 
element  spacing  in  order  to  increase  the  element  coupling.  This  was  done  intentionally 


VIII.  Microstrip  Array  Antenna  Experiments 


129 


Experimental  array  of  eight  quarter-wavelength  microstrip  patches. 
This  board  is  refered  to  as  the  array  board.  The  left  most  input  port 
is  No.  I. 


in  an  attempt  to  make  mutual  coupling  effects  more  readily  measurable.  The  individual 
elements  are  5.08cm  by  1.89cm  rectangular  patches  with  sixteen  shorting  pins  placed 
approximately  1.68cm  (0.25,1,,)  from  the  opposite  radiating  edges.  A  detailed  sketch  of  a 
patch  is  shown  in  Fig.  8.1-2.  The  elements  are  fed  by  50  Q  transmission  lines  which  are 
matched  to  the  patches  with  a  132  Q  quarter- wavelength  transformer  at  the  center  of 
radiating  edges. 

The  transmission  lines  attached  to  non-radiating  edges  of  the  patches  are  probe 
circuits  intended  for  measurements  of  active  element  excitation  coefficients.  The  probe 
circuits  are  made  of  158  Q  quarter- wavelength  transmission  lines  and  50  Q.  chip 
resistors.  The  input  impedance  of  the  probe  circuits  looking  in  from  the  patches  are  kept 
high  (approximately  1  KQ)  with  the  quarter-wavelength  transformers  so  that  the  probes 
would  not  load  the  patches  which  have  low  input  impedance  because  of  shorting  pins. 
The  chip  resistors  are  added  because  the  probe  circuits  have  high  input  impedance 
looking  in  from  detector  ports. 

Two  feed  network  boards  were  fabricated  to  test  different  cases  of  array  excitation: 
a  sum/difference  feed  (Fig.  8.1-3)  and  an  ordinary  endfire  feed  (Fig.  8.1-4).  A  sum 
pattern  requires  equal  magnitude  and  equal  phase  excitation,  while  a  difference  pattern 
requires  one  half  of  the  array  to  be  excited  1 80°  out  of  phase  from  the  other  half.  The 
sum  feed  and  the  difference  feed  are  integrated  and  constructed  on  the  same  substrate. 
The  combined  sum/difference  network  uses  a  quarter-wavelength  hybrid  with  two  input 
ports  attached  to  a  corporate  feed  network.  Depending  on  the  input  port  chosen,  the 
sum  pattern  or  the  difference  pattern  excitation  is  obtained. 

An  ordinary  endfire  pattern,  on  the  other  hand,  requires  equal  magnitude  and 
progressive  phase  shift  4>n  —  ±  n/Jd  where  d  is  the  inter-element  spacing.  The  endfire  feed 
is  constructed  with  seven  3  dB  power  dividers  connected  in  the  form  of  a  corporate  feed 
with  additional  transmission  lines  on  the  right  branch  of  the  power  dividers.  The 
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To  Detector 


Shorting  Pins 


50  ft  Grip  Resistor 


132  ft 


Ad/4 


50  ft  Feed  Line 


Figure  8.1-2.  Geometry  of  the  quarter-wavelength  radiating  element  with  feed  line 
and  probe  circuit  used  in  the  experimental  array. 
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transmission  lines  are  designed  to  provide  1 13°(0.312  >1,)  phase  shift  between  elements 
at  2733  MHz.  For  example,  the  difference  in  the  length  of  transmission  lines  between 
port  1  and  port  2  is  d,*  2.432  cm  as  shown  in  Fig.  8.1-4.  The  interelement  phasing 
between  two  adjacent  ports,  a  ,  is  obtained  as  follows: 


a 


£gdg 

113° 


360° 

7.706cm 


2.432cm 


(8.1-1) 


where  /?,  and  are  propagation  constant  and  wavelength  in  the  substrate,  respectively. 

The  array  board  and  the  feed  networks  are  connected  using  coaxial  cables  between 
the  feed  output  ports  and  the  antenna  element  ports.  This  design  allows  measurement 
of  S-parameters  as  well  as  placement  of  a  compensation  network  between  the  array  and 
the  feed  network  (the  compensation  network  was  discussed  in  Section  2.4).  In  addition 
to  the  above  three  microstrip  networks,  two  other  boards  were  constructed,  one 
containing  both  the  array  network  and  an  endfire  feed  and  the  other  containing  both  the 
array  and  the  sum/difference  feed.  These  boards  were  built  to  verify  that  the  coaxial 
cables  between  the  array  and  the  feed  do  not  affect  the  radiation  pattern. 


8.2  Measurements  with  the  Uncompensated  Microstrip  Arrays 

This  section  presents  measured  and  predicted  radiation  patterns  for  the 
uncompensated  experimental  array  consisting  of  the  array  board  connected  to  a  feed 
network  board  (these  boards  were  described  in  the  previous  section).  All  measurements 
(S-parameters  and  far-field  patterns)  which  were  conducted  at  PSL  were  taken  at  the 
design  frequency  2733  MHz  and  at  side  frequencies  2723  MHz  and  2743  MHz  to  study 
frequency  variation  effects.  Element  patterns  were  measured  only  at  2733  MHz  because 
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the  element  patterns  do  not  change  appreciably  with  frequency.  The  individual  element 
patterns  of  the  experimental  array  when  operating  in  array  environment  are  shown  in 
Fig.  8.2-1.  These  element  patterns  are  measured  by  feeding  one  element  while  all  others 
were  50  Q  loaded.  Figure  8.2-1  shows  that  the  element  patterns  are  significantly 
different.  This  suggests  that  an  average  element  pattern  cannot  be  used  to  calculate 
predicted  radiation  patterns  for  the  entire  array  as  in  (2.3-10).  For  the  complete  accuracy 
the  difference  in  the  element  patterns  must  be  incorporated  in  far-field  calculations.  This 
was  not  done  and  proved  to  limit  the  achievable  compensation.  Assuming  that  the 
measured  element  patterns  are  close  to  that  of  active  element  pattern  in  fully  excited 
array,  the  radiation  pattern  can  be  approximated  using  (2.3-9). 

Before  proceeding  to  actual  compensation  it  is  desireable  to  see  how  well  array 
behavior  can  be  predicted;  i.e.  we  first  study  analysis  (uncompensated)  before  synthesis 
(compensated).  Figures  8.2-2  to  8.2-4  show  the  measured  and  predicted  E-plane 
radiation  patterns  of  the  uncompensated  array  for  sum,  difference  and  endfire  excitation 
at  2733  MHz.  The  predicted  patterns  are  calculated  from  element  currents  which  are 
obtained  from  measured  array  and  feed  network  S-parameters  using  MCAP.  The 
currents  are  given  in  Table  8.2-1.  The  element  patterns  used  in  these  predicted  patterns 
are  those  of  Fig.  8.2-1.  In  addition,  the  uncompensated  radiation  patterns  shown  in 
Figs.  8.2-2(a)  to  8.2-4(a)  are  compared  to  predicted  patterns  calculated  from  ideal 
excitation  currents  which  are  shown  in  Figs.  8.2-2(c)  to  8.2-4(c).  The  ideal  patterns  are 
obtained  when  array  elements  are  excited  with  the  desired  currents.  These  patterns,  thus, 
represent  what  can  be  achieved  through  mutual  coupling  compensation.  The  excitations 
for  the  ideal  patterns  are  given  in  Table  8.2-2.  The  endfire  patterns  (Fig.  8.2-4)  have  the 
main  beams  away  from  the  endfire  direction  ($  =  90®).  This  is  due  to 
non-omnidirectional  element  pattern  in  E-plane  of  the  microstrip  array. 
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8.1-1.  liach  pattern  was  measured  with  the  remaining  seven  elements 

f  A  f\  1 _ .1  _  1  w 


la)  -30  -20  -10  OdB 


Figure  8.2-2,  Comparison  of  measured  and  calculated  radiation  patterns  for  the 
array  with  sum  pattern  feed  network,  a)  measured,  b)  calculated  and 
c)  ideal  patterns. 
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Figure  8.2-4.  Comparison  of  measured  and  calculated  radiation  patterns  to  for  the 
array  with  endfire  feed  network,  a)  measured,  b)  calculated  and  c) 
ideal  patterns. 
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Table  8.2-1.  Calculated  currents  for  the  predicted  patterns  of  Figs.  8.2-2  through 
8.2-4. 


Element 

Number 

Sum 

Array 

Difference 

Array 

Endfire 

Array 

1 

0.68664/60.74° 

0.45911/132.37° 

0.34980/221.87° 

2 

0.66823/42.03° 

0.40396/115.25° 

0.46687/81.22° 

3 

0.96257/60.61° 

0.84871/140.52“ 

0.89593/-28.680 

4 

1.00000/59.25° 

1.00000/156.20° 

0.68871/238.61“ 

5 

0.71603/35.55° 

0.80937/-29.020 

1.00000/129.15° 

6 

0.69268/30.96° 

0.6951 1/-54. 34“ 

0.78323/11.80° 

7 

0.60754/63.25° 

0.40664/-41.27® 

0.77976/267.39° 

8 

0.74779/67.59° 

0.56739/-34.52® 

0.85398/152.65° 
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Table  8.2-2.  Excitation  for  'Ideal'  patterns  for  Figs.  8.2-2  through  8.2-4. 


Figures  8.2-2  to  8.2-4  clearly  show  that  the  order  of  error  between  measured  and 
predicted  patterns  are  close  to  the  difference  between  measured  patterns  and  ideal 
patterns.  The  deviation  of  predicted  patterns  from  measured  patterns  apparently  comes 
from  errors  in  measured  S-parameters.  This  suggests  that  compensation  for  mutual 
coupling  in  the  experimental  array  may  be  difficult  to  demonstrate  in  microstrip.  The 
small  improvements  seen  in  the  ideal  radiation  patterns  from  the  uncompensated 
patterns  are  due  to  variation  in  the  element  patterns  as  shown  in  Fig.  8.2-1.  For  example, 
there  are  up  to  8  dB  differences  in  the  element  patterns  at  60°  off  broadside  which  is 
close  to  the  location  of  the  main  beam  in  the  endfire  array.  The  element  pattern 
variation  occurs  primarily  from  ground  plane  edge  effects. 

8.3  Compensation  Networks  for  the  Experimental  Arrays 

Both  attenuator/phase  shifter  and  power  divider  compensation  networks  have  been 
calculated  for  the  experimental  array  in  all  three  cases  of  excitations.  The  compensation 
network  values  were  calculated  from  measured  S-parameters  using  SANE  for 
attenuator/phase  shifter  network  and  SANE-PODCON  for  power  divider  network. 
Since  the  accuracy  of  the  S-parameters  was  shown  in  the  previous  section  to  be 
questionable,  the  compensation  networks  calculated  here  may  be  somewhat  inaccurate. 
However,  the  S-parameters  represent  typical  mutual  coupling  phenomena  in  an  array. 
Thus,  the  networks  obtained  are  representative  of  the  mutual  coupling  compensation 
technique. 
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8.3.1  Attenuator/ phase  shifter  compensation  networks 


This  section  presents  attenuator/phase  shifter  compensation  networks  for  the  sum, 
difference  and  endfire  array.  The  compensation  networks  are  calculated  from  measured 
S-parameters  of  the  sum/difference  and  the  endfire  feed  networks  and  measured 
S-parameters  of  the  microstrip  elements.  The  variables  for  this  example  are  defined  as 
follows: 

Attenuation[dB]  —  -20  log(  |  X 1 ) 

(8.3-1) 

Phase  Shift  [°]  =  [X 

The  results  are  shown  in  Table  8.3-1  which  are  the  most  efficient  compensation  networks 
when  phase  of  generator  G  is  varied  for  every  10°.  The  relationship  between  the 
generator  excitation  and  the  efficiency  of  attenuator/phase  shifter  compensation  network 
was  discussed  in  Section  5.3.  The  predicted  patterns  for  the  compensated  array  for  the 
three  cases  of  excitations  are  equivalent  to  patterns  with  the  ideal  excitation  currents. 
The  ideal  patterns  are  shown  in  Figs.  8.2-2(c),  8.2-3(c),  and  8.2-4(c)  for  the  sum, 
difference  and  endfire  array,  respectively. 

8.3.2  Power  divider  compensation  network 

Similar  to  the  attenuator/phase  shifter  compensation,  the  power  divider 
compensation  networks  are  calculated  from  measured  S-parameters  of  the  microstrip 
array  for  the  synthesis  of  a  sum,  a  difference  and  an  endfire  pattern.  The  variable  power 
divider  networks  are  defined  in  terms  of  complex  variables  as  shown  below 
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Table  8.3-1.  Attenuator/phase  shifter  compensation  networks  for  the  microstrip  array. 


case 

Sum 

Difference 

Endfire 

Element 

Attenuator 

Ph.Shift 

Attenuator 

Ph.Shift 

Attenuator 

Ph.Shift 

x, 

2.201  dB 

243.91° 

2.536  dB 

0.000  dB 

257.34° 

x, 

1.075  dB 

272.19° 

0.631  dB 

226.19° 

2.468  dB 

258.62° 

X, 

2.372  dB 

254.09° 

1.210  dB 

220.40“ 

5.206  dB 

254.51° 

2.628  dB 

6.193  dB 

184.24“ 

5.636  dB 

275.04° 

2.041  dB 

255.82° 

6.299  dB 

193.81° 

2.532  dB 

273.71° 

0.302  dB 

258.91° 

0.000  dB 

224.01° 

2.756  dB 

281.55° 

0.000  dB 

266.99“ 

0.004  dB 

218.41° 

5.401  dB 

261.24° 

2.199  dB 

239.67° 

3.482  dB 

206.21° 

6.401  dB 

286.36“ 

G 

5.55346 

40.00° 

5.23843 

0.00“ 

4.64400 

240.00° 

VIII.  Microstrip  Array  Antenna  Experiments 


145 


Fraction  of  power  to  left  branch  =  1  -  K  =  1  -  |  X I 

Fraction  of  power  to  right  branch  —  K  —  |  X I  (8.3-2) 

Phase  Shift  [°]  =  jX 

The  results  are  shown  on  Table  8.3-2.  The  compensation  networks  were  calculated  using 
quadrature  hybrid  power  dividers  discussed  in  Section  4.4  with  transmission  line  length 
tfe  =»  >1/8  (see  Figure  4.4-2).  The  calculated  patterns  for  the  compensated  arrays  are 
equivalent  to  the  ideal  patterns  which  are  shown  in  Figs.  8.2-2(c),  8.2-3(c)  and  8.2-4(c). 
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Table  8.3-2.  Results  of  power  divider  compensation  networks  for  the  microstrip  array. 


case 

Sum 

Difference 

Endfire 

K 

Ph.Shift 

K 

Ph.Shift 

K 

Ph.Shift 

X, 

0.55230 

38.17° 

0.60849 

19.01° 

0.40649 

-113.56° 

X, 

0.49141 

8.11° 

0.10034 

-11.66° 

0.29009 

-159.83° 

X, 

0.56489 

3.02° 

0.84312 

45.89° 

0.43559 

-151.01° 

0.29576 

-17.47° 

0.28341 

-14.03° 

0.01106 

-150.50° 

0.43517 

19.76° 

0.30678 

23.92“ 

0.04399 

142.75° 

0.52865 

-5.92° 

0.64275 

38.63“ 

0.14548 

145.97“ 

0.48339 

28.27° 

0.49484 

160.83° 

0.21778 

-70.76° 

G 

4.35587 

116.86° 

4.11564 

124.33° 

3.04167 

-138.73° 

i 

: 
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8.4  Synthesis  of  a  35  dB  Dolph-Chebyshev  Pattern 


As  shown  in  Section  8.2,  the  differences  between  the  uncompensated  patterns  and 
the  compensated  patterns  are  small.  A  better  way  to  demonstrate  the  effect  of  the 
mutual  coupling  compensation  on  the  microstrip  array  is  to  synthesize  a  low  side  lobe 
pattern.  A  35  dB  Dolph-Chebyshev  pattern  was  chosen  for  this  purpose.  The  required 
element  currents  are  [S-2] 

1,-0.19154  /  0.0°  1,  -  1.00000  /  0.0° 

I,  =  0.46362  /  0.0°  I*  -  0.78425  /  0.0° 

(8.4-1) 

I3  -  0.78425  /  0.0°  I7  -  0.46362  /  0.0° 

1-1.00000  /  0.0°  1,-0.19154  /  0.0° 

A  feed  network  for  the  above  current  distribution  that  does  not  include  mutual 
coupling  compensation  (blind  design  feed)  was  constructed  from  the  sum/difference  feed 
network  (using  sum  pattern  input  port)  with  six  attenuators  corresponding  to  the 
magnitude  of  required  currents  attached  to  the  output  ports.  See  Table  8.4-1.  The 
resulting  element  currents  calculated  using  measured  S-parameters  of  the  feed  network 
and  the  array  elements  are  shown  in  Table  8.4-2.  The  measured  and  predicted  patterns 
are  shown  in  Fig,  8.4-1.  The  predicted  patterns  were  again  computed  using  individual 
element  patterns  of  Fig.  8.2-1. 

To  demonstrate  that  a  feed  network  that  is  tuned  to  produce  the  desired  currents 
does  not  compensate  for  mutual  coupling  effect,  a  network  which  corrects  the  magnitude 
and  phase  error  of  the  feed  network  was  designed.  This  is  designated  as  partial  blind 
design.  The  attenuators  and  phase  shifters  required  by  partial  blind  design  are  shown  in 
Table  8.4-1.  These  parameters  are  calculated  using  SANE  with  measured  S-parameters 
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Table  8.4-1.  Attenuator/ phase  shifter  compensation  network  for  a  35  dB  Dolph 

Chebyshev  array  using  the  eight  element  microstrip  array  and  the  sum 
feed  network. 


Design 

Blind 

Partial  Blind 

Compensated 

Element 

Attenuator 

Ph.  Shift 

Attenuator 

Ph.Shift 

Attenuator 

Ph.Shift 

x, 

14.355 

dB 

0.008 

14.271 

dB 

-0.95* 

11.540 

dB 

-135.57° 

x, 

6.667 

dB 

0.00“ 

6.178 

dB 

3.053 

dB 

X, 

2.111 

dB 

0.00“ 

2.516 

dB 

2.768 

dB 

-144.95“ 

0.000 

dB 

0.00“ 

0.405 

dB 

-6.07“ 

0.261 

dB 

-142.68“ 

0.000 

dB 

0.00° 

0.000 

dB 

-2.96“ 

0.000 

dB 

-137.71* 

2.111 

dB 

0.00“ 

1.209 

dB 

-3.98“ 

1.259 

dB 

-134.74* 

6.667 

dB 

0.00* 

6.095 

dB 

-6.97“ 

1.602 

dB 

-132.37° 

14.355 

dB 

0.00“ 

13.857 

dB 

-4.00“ 

12.013 

dB 

-134.48° 

G 

2.82843 

0.00“ 

3.19917 

0.00° 

3.94536 

80.00° 
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of  the  sum  feed  network  while  setting  all  S-parameters  of  the  elements  to  zero.  The 
predicted  radiation  pattern  for  the  blind  design  array  is  shown  in  Fig.  8.4*2.  It  clearly 
shows  that  mutual  coupling  among  elements  cannot  be  neglected  in  design  of  an  array. 

An  attenuator/phase  shifter  compensation  network  was  calculated  from  the 
S*parameters  of  the  sum  feed  network  measured  at  PSL.  The  results  of  the  calculation 
using  SANE  for  the  required  compensation  network  parameters  are  summarized  in 
Table  8.4-1.  The  compensation  network  was  constructed  at  PSL  using  attenuators  and 
coaxial  cables.  Attenuators  which  closely  approximate  required  values  were  chosen  first. 
Then,  coaxial  cables  were  cut  to  achieve  the  required  phase  shifts  and  the  difference  in 
required  and  actual  attenuation  using  the  loss  introduced  by  the  cables.  The  normalized 
E-plane  radiation  pattern  of  the  compensated  array  is  presented  in  Fig.  8.4-3.  The 
currents  used  in  the  calculation  of  predicted  pattern  are  shown  in  Table  8.4-2. 

Comparison  of  measured  patterns  between  the  blind  design  array  (Fig.  8.4-1)  and 
the  compensated  array  (Fig.  8.4-3)  shows  some  improvement  in  sidelobe  level  and  main 
beam  shape.  However,  the  sidelobes  of  the  compensated  array  are  not  as  low  as  we  have 
predicted  due  to  variation  in  the  individual  element  pattern.  Nevertheless,  the  result 
offers  encouragement  that  experimental  verification  may  be  possible  for  an  array  with 
small  element  pattern  variations. 


VIII.  Microstrip  Array  Antenna  Experiments 


150 


Figure  8.4-2.  Calculated  radiation  patterns  of  a  partial  blind  design  35  dB 
Dolph-Chebyshev  array. 
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Table  8.4*2.  Calculated  element  currents  of  the  blind  designs  and  compensated  35  dB 
Dolph-Chebyshev  array  using  the  eight  element  microstrip  array  and  the 
sum  feed  network  with  G  =  1/0.0° 


Element 

Blind 

Partial  Blind 

Compensated 

1 

0.03537/77.20° 

0.03655/77.76° 

0.05359/-79.130 

2 

0.08540/61.26° 

0.06978/53.85° 

0.1 1785/-80.030 

3 

0.18667/57.07° 

0.18699/49.96° 

0.20661 /-83.62° 

4 

0.22160/56.08“ 

0.22419/53.31“ 

0.25848/-84.590 

5 

0.17864/40.42° 

0.18456/40.35° 

0.2688 1/-89.37° 

6 

0.14962/44.91° 

0.15711/32.09“ 

0.2 1 686/-86. 1 6° 

7 

0.07063/61.41“ 

0.05719/66.48° 

0.12240/-83.228 

8 

0.03916/68.61“ 

0.04427/71.86° 

0.04994/-88.440 
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IX.  Numerical  Studies  with  Wire  Antennas 


9.1  Introduction 

In  order  for  the  synthesis  techniques  developed  here  to  work,  accurate  knowledge 
of  all  coupling  effects  is  required.  This  is  an  obvious  restriction;  coupling  effects  must 
be  quantified  before  compensation  values  can  be  determined.  In  Chapter  8  experimental 
means  were  used  to  characterize  the  coupling  (both  array  and  feed  coupling).  The 
results  did  not  adequately  demonstrate  the  improvements  possible  with  the  synthesis 
methods.  Another  approach  which  offers  knowledge  of  the  antenna  array  mutual 
coupling  is  to  use  moment  methods  with  wire  antennas.  The  results  of  such  a  study  are 
presented  in  this  chapter.  The  array  is  modeled  using  a  moment  method  program 
discussed  fully  in  Section  9.2  and  the  feed  is  modeled  using  SANE  and  SANE-PODCON. 
Details  of  the  procedure  are  presented  in  Section  9.3.  Both  arrays  of  parallel  dipoles  and 
parallel  monopoles  over  a  finite  ground  plane  are  studied.  Dipole  array  synthesis  and 
monopole  array  synthesis  are  discussed  in  Section  9.4  and  Section  9.5,  respectively.  In 
both  cases  the  compensation  process  is  successfully  demonstrated. 
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9.2  Using  Moment  Methods  in  Numerical  Experiments 


The  method  of  moments  is  a  general  technique  used  to  approximately  solve 
equations  of  the  form 

F(g)-h  (9.2-1) 

where  F  is  a  known  linear  operator,  h  is  a  given  excitation  function,  and  g  is  an  unknown 
response  function.  Expansion  of  the  unknown  functions  in  a  linear  combination  of 
known  basis  function  allows  (9.2-1)  to  be  expanded  to  a  system  of  N  linear  equations 
[B-9J: 


N 

£cn  <  Wm,  F(gJ  >»<wm,h>  m  -  1 ,2 . N  (9.2-2) 

n— 1 

where  {gj  are  expansion  functions  for  the  unknown  function  g  with  coefficients  {c„}, 
{wm}  are  weighting  functions,  and  <  a,b>  represents  the  inner  product  of  a  and  b.  The 
system  of  N  equations  can  be  written  in  matrix  form  as 

[F-JECJ  -  [h  J  (9.2-3) 

where 

F,j-  <wifF(gj)> 

hj  —  <Wj,  h> 

Equation  (9.2-3)  can  be  solved  for  the  unknowns  {c„}  using  matrix  inversion: 

[«J  -  [FmJ"'[hJ  (9.2-4) 
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The  unknown  response  function  g  in  (9.2-1)  can  now  be  approximated  by 


(9.2-5) 


Using  the  method  of  moments  the  electric  field  integral  equation  of  the  form 


J 


I(s')K(s,  s')ds'  *  -  E*(s) 


(9.2-6) 


can  be  solved  for  unknown  current  l(s’)  where  kernel  K(s,  s’)  is  specified  by  the  geometry 
of  the  antennas  and  the  incident  electric  field  E'(s)  is  specified  by  the  given  excitation. 
For  example,  one  of  the  common  integral  equations  used  in  treatment  of  thin  wire 
structures  was  derived  by  Pocklington.  His  integral  equation  for  a  z-directed  wire  is  given 
as  [S-2] 


1 

j  O)t0 


d2Mz,z') 

dz2 


+  /?V(z,z') 


(9.2-7) 


The  function  ^(z,  z')  is  the  free  space  Green's  function 


Mz,  z') 


e-^R 

4ttR 


(9.2-8) 


where  R  is  distance  between  the  observation  point  z  and  the  source  location  z'.  Once 
the  induced  currents  on  the  wires  are  determined  using  the  method  of  moment  the 


radiated  electric  fields  are  computed  by 
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There  are  numerous  computer  programs  implementing  method  of  moments  for 
various  antenna  structures.  In  this  study  two  programs  are  used.  The  Electromagnetic 
Surface  Patch  code  (ESP),  which  was  developed  at  the  Ohio  State  University 
Electroscience  Laboratory,  is  a  computer  implementation  of  the  moment  method 
technique  for  solving  geometries  consisting  of  thin  wires  and  perfectly  conducting 
polygonal  plates  [N-2].  ESP  uses  piecewise  sinusoidal  functions  as  both  expansion 
function  and  weighting  function  [N-2],  The  procedure  is  known  as  the  piecewise 
sinusoidal  Galerkin's  method  which  is  computationally  more  efficient  than  other  pairs 
of  expansion  and  weighting  functions  [S-2].  The  program  approximates  the  excitation 
by  delta  gap  voltage  generators  which  can  be  located  anywhere  on  the  wires  or  at  the 
attachment  points  between  plates  and  wires  [N-2].  Then  the  excitation  function  E'(s)  is 
zero  everywhere  (since  tangential  electric  field  must  be  along  a  wire)  except  across  the 
feed  points. 

The  second  moment  method  program  used  in  this  study,  MININEC  II  written  by 
Davis  [D-10],  is  an  improved  version  of  the  Mini-Numerical  Electromagnetics  Code 
(MININEC).  MININEC  is  a  personal  computer  version  of  more  extensive  Numerical 
Electromagnetics  Code  (NEC).  MININEC  II  determines  the  current  distribution  and 
far  field  pattern  of  wire  geometries  in  free  space  and  wire  geometries  over  an  infinite 
ground  plane.  The  program  uses  pulse  function  as  both  expansion  and  weighting 
function  [D-10]. 

Before  ESP  is  used  in  numerical  experiments  to  calculate  coupling  coefficients  of 
array  antennas  it  is  necessary  to  verify  that  the  mutual  impedance  values  obtained  from 
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the  program  are  reasonable.  Figure  9.2-1  compares  the  mutual  impedance  between  two 
parallel  half  wavelength  dipoles  for  various  element  spacings.  The  dotted  curves 
represent  the  values  calculated  by  ESP  for  dipoles  with  radius  a  =  0.005^0.  The  solid 
curves  represent  the  mutual  impedance  between  two  ideal  parallel  dipoles  of  length 
L  —  20/2  which  are  obtained  from  [B-9] 


R,j 

X,a 


n 

An  sin2(/?L/2) 
V 

An  sin2(/?L/2) 


C2Ci(uo)  —  Ci(Uj)  —  Ci(u2)3 
C2Si(u0)-Si(u,)-Si(u2)2 


(9.2-10) 


where 


Uo  =  /id 

Ui  =  P(yJ  d2  +  L2  +  L) 
u2  -  /?( -  L) 

Ru  and  X,a  are  real  and  imaginary  parts  of  the  mutual  impedance  Z,„  respectively.  The 
functions  Ci  and  Si  are  cosine  and  sine  integrals.  The  differences  between  ideal  dipole 
results  and  the  values  calculated  by  ESP  for  close  element  spacings  are  due  to  the  finite 
radius  of  dipoles  accounted  for  with  the  ESP  model. 

The  programs  ESP  and  MININEC  II  are  also  used  to  predict  S-parameters  of  an 
array  of  eight  monopoles  over  a  ground  plane.  The  monopoles  have  length 
L  =  0.2375<i„  radius  a  =*  0.01  U0  and  are  spaced  d  =*  0.5>10  apart.  The  ground  plane  is  a 
5.5 52„  x  5.55^0  rectangular  conductor.  In  the  MININEC  calculation,  the  finite  ground 
plane  was  replaced  with  an  infinite  ground  plane. 

Figure  9.2-2  shows  comparison  of  the  calculated  to  measured  S-parameters  of  the 
array.  The  S-parameters  of  the  array  are  calculated  from  the  outputs  of  ESP  and 
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Mutual  Impedance,  Ohm 


1 


Figure  9.2-1.  Comparison  of  mutual  impedance  between  two  parallel  a,/2  dipoles. 

Ideal  dipoles  (solid  curves).  Calculated  using  ESP  (dotted  curves). 
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Figure  9.2-2.  Comparison  of  the  calculated  to  the  measured  S-paramctcrs  of  an 
eight  element  monopolc  array.  Measured  (solid  curves).  Calculated 
using  l:SP  (dotted  curves).  Calculated  using  M1N1N12C  II  (dashed 


MININEC  II  using  steps  1)  and  2)  of  the  numerical  experiment  procedures  described  in 
next  section.  It  is  apparent  from  the  plots  that  MININEC  II  predicts  measured  results 
more  closely  than  ESP.  ESP  is  especially  poor  at  predicting  the  return  loss,  S„,  of  the 
monopole  array.  This  result  suggest  that  the  uncompensated  patterns  of  a  monopole 
array  calculated  in  the  numerical  experiments  in  Section  9.5  may  be  worse  than  the 
actual  radiation  patterns  of  the  array. 

9.3  Numerical  Experiment  Procedures 

The  numerical  experiments  in  the  next  two  sections  were  performed  using  the 
following  steps. 

1).  The  port  to  port  mutual  admittance  of  array  antennas  are  calculated  using  ESP.  The 
admittance  parameters  are  obtained  by  exciting  one  element  and  shorting  all  other 
elements  in  an  array: 

Ymn  =  — •  where  Vj  =  0  for  j  #  n  (9.3-1) 

vn 

where  Im  is  the  appropriate  port  current  obtained  from  the  moment  method  solution. 
The  calculation  must  be  repeated  N  times  for  an  N  element  array  to  obtain  the  complete 
Y-matrix. 


2).  The  admittance  matrix  of  the  array  is  converted  to  S-matrix  by  [G-l] 

[S]  -  -  Y][Y„  +  YT'cVYcT]  (9.3-1) 
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where  [Y*]f  [>/y7]  and  [7%"]  are  diagonal  matrices  with  elements  given  by  Y#  , 
and  Jz^  -»  1/VyT  ,  respectively.  Yw  represents  the  characteristic  admittance  of  the 
transmission  line  attached  to  port  i  of  the  network. 


3)  Calculate  the  required  element  currents  using  an  ordinary  array  synthesis  technique. 

4) .  Calculate  the  S-matrix  of  the  initial  feed  network  that  does  not  include  mutual 
coupling  compensation  using  MCAP  for  the  element  currents  obtained  in  step  3  (blind 
design). 

5) .  Calculate  the  actual  element  currents  and  voltages  using  MCAP  for  the  array  and  the 
feed  network  obtained  in  the  previous  steps.  Then  calculate  blind  design  far-field  pattern 
using  ESP  with  the  voltages  as  excitation  (in  place  of  a  feed  network). 

6) .  Compute  a  compensation  network  using  SANE  or  modify  the  feed  network  using 
SANE-PODCON  to  compensate  for  mutual  coupling. 

7) .  Calculate  the  element  voltages  using  MCAP. 

8) .  Compute  the  compensated  far-field  patterns  using  ESP  with  the  voltages  of  step  7 
as  excitations  (in  place  of  a  compensated  feed  network). 


i 
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9.4  An  Eight  Element  Dipole  Array  Example 


The  eight  element  linear  array  considered  in  this  section  consists  of  z-directed 
parallel  dipoles  of  length  L  =  0.4662„  and  radius  a  =  0.0052„.  The  elements  are  spaced 
d  =  0.500/l0  apart  along  x-axis  and  are  fed  at  the  center  of  the  dipoles.  The  geometry  and 
the  orientation  of  the  array  is  shown  in  Fig.  9.4-1.  An  operating  frequency  of  300  MHz 
was  chosen  for  the  ease  of  calculation  =1  m).  A  30  dB  Dolph-Chebyshev  pattern 
with  the  main  beam  steered  to  30°  off  broadside  is  the  design  objective.  The  required 
element  currents  found  using  the  step  3  in  Section  9.3  are  {S-2J 

I,  -  0.26222/3 13a  I5  -  1.00000/-45° 

It  =»  0.51875/225°  I*  =-  0.8 1 196/- 1 35° 

(9.4-1) 

13  -  0.81 196/135°  !■>  ■»0.51875/-225° 

14  -  1 .00000/45°  I#  -  0.26222/-315° 

The  feed  network  for  the  above  current  distribution  that  does  not  include  mutual 
coupling  compensation  (blind  design)  consists  of  seven  2-way  variable  power  divider 
networks  with  the  required  phase  shift  to  steer  the  main  beam.  The  required  power 
dividers  and  phase  shifters  are  listed  in  Table  9.4-2.  The  element  currents  that  result 
from  the  blind  design  feed  network  and  the  corresponding  E-plane  and  H-plane  radiation 
patterns  are  shown  in  Table  9.4-3  and  Figure  9.4-5,  respectively. 
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9.4.1  Attenuator/phase  shifter  compensation 


In  the  attenuator/phase  shifter  compensation  network,  an  attenuator  in  series  with 
a  phase  shifter  was  placed  between  the  feed  and  each  element  (see  Fig.  9.4-2).  In  this 
compensation  technique  the  generator  phase  was  changed  at  10°  increment  to  find  the 
most  efficient  network.  (See  discussion  in  Section  5.3).  The  required  compensation 
network  parameters  and  other  results  using  SANE  are  summarized  in  Table  9.4-1.  The 
normalized  E-plane  and  H-plane  radiation  patterns  of  the  compensated  array  are  plotted 
in  Fig.  9.4-5.  The  radiation  patterns  of  the  blind  design  array  are  also  plotted  for 
comparison. 

9.4.2  Power  divider  compensation 

For  variable  power  divider  network  compensation  the  power  ratio  and  phase  shifts 
are  modified  from  those  of  the  blind  design  to  achieve  mutual  coupling  compensation. 
The  power  dividers  and  phase  shifters  are  formed  in  the  corporate  feed  network 
configuration  as  shown  in  Fig.  9.4-3.  The  required  power  divider  and  phase  shifter 
values  are  given  in  Table  9.4-2.  The  radiation  patterns  of  the  array  with  the  modified 
feed  network  along  with  those  of  blind  design  array  are  shown  in  Fig.  9.4-8. 
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Table  9.4-1.  Required  attenuators  and  phase  shifters  for  the  30  dB  Dolph-Chebyshev 
8-element  dipole  array  example. 


Attenuation 

Phase  Shift 

3.119  dB 

31.94° 

2.425  dB 

29.81° 

1.509  dB 

34.53° 

1.080  dB 

18.70° 

0.000  dB 

26.04° 

0.831  dB 

0.36° 

0.329  dB 

18.64° 

2.307  dB 

-12.84“ 

2.51264 

150.00° 
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Figure  9.4-3.  Variable  power  divider  network  for  an  eight  element  dipole  array. 


Table  9.4-2.  Required  variable  power  dividers  and  phase  shifters  for  30dB 
Dolph-Chebyshev  8-element  dipole  array  example. 


Blind  Design 

Compensated  Design 

No. 

K 

Ph.  Shift 

K 

Ph.  Shift 

1 

0.79649 

-90.00° 

0.81382 

-80.76° 

2 

0.60267 

-90.00° 

0.47938 

-100.79° 

3 

0.39733 

-90.00° 

0.22670 

-105.57° 

4 

0.20351 

-90.00° 

0.06047 

-112.71° 

5 

0.83083 

180.00° 

0.85703 

-170.25° 

6 

0.16917 

180.00° 

0.15481 

-179.89° 

7 

0.50000 

0.00° 

0.49279 

-2.37° 

G 

1.99857 

-135.00° 

2.23865 

171.65° 
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9.4.3  The  frequency  sensitivity  of  the  compensation  networks 


The  compensation  networks  derived  in  the  previous  subsections  are  of  no  practical 
value  if  they  are  very  sensitive  to  changes  in  the  frequency.  Thus,  it  is  necessary  to  test 
the  frequency  sensitivity  of  the  networks.  In  this  section  the  radiation  patterns  using  the 
compensation  networks  in  the  previous  two  sections  are  calculated  at  294  MHz  and  306 
MHz. 

In  the  attenuator/phase  shifter  network,  the  compensation  variables  are  kept 
constant  while  S-parameters  of  the  blind  design  feed  network  array  are  modified 
according  to  their  frequency  behavior.  The  array  S-matrices  are  also  recomputed  at  each 
frequency  using  ESP.  In  the  case  of  the  variable  power  divider  compensation  network, 
the  S-matrices  of  the  modified  feed  network  are  calculated  at  three  frequencies  using 
MCAP  to  include  the  frequency  variation  effects  of  the  power  dividers  discussed  in 
Section  4.4.3.  The  array  S-matrices  are  the  same  as  those  used  for  the  attenuator/phase 
shifter  case. 

Figures  9.4-4  through  9.4-6  show  the  normalized  E-plane  and  H-plane  patterns  of 
the  blind  design  array  and  the  array  with  attenuator/phase  shifter  compensation  at  294 
MHz,  300  MHz  and  306  MHz,  respectively;  these  represent  a  4%  frequency  swing. 
Figures  9.4-7  through  9.4-9  show  the  compensated  patterns  using  variable  power  divider 
networks  at  the  three  frequencies.  The  figures  show  that  though  the  compensated 
patterns  deviate  from  the  desired  pattern,  they  are  still  closer  to  the  desired  pattern  than 
the  blind  design  radiation  patterns. 

The  frequency  sensitivity  investigation  in  this  section  showed  that  the  synthesis 
techniques  are  not  so  sensitive  that  a  slight  frequency  change  eliminates  the 
improvements  achieved  at  the  design  frequency. 
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Figure  9.4-4.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  blind  design  (dotted)  and  the  attenuator/phase 
shifter  network  compensation  (solid)  at  294  MHz. 
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Figure  9.4-5.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  blind  design  (dotted)  and  the  attenuator/phase 
shifter  network  compensation  (solid)  at  300  MHz  the  design 
frequency. 
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Figure  9.4-6.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  blind  design  (dotted)  and  the  attenuator/phase 
shifter  network  compensation  (solid)  at  306  MHz. 
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Figure  9.4-7. 


H-plane 

Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  the  blind  design  (dotted)  and  the  variable  power 
divider  network,  compensation  (solid)  at  294  MHz. 
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Figure  9.4-8.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  the  blind  design  (dotted)  and  the  variable  power 
divider  network  compensation  (solid)  at  300  MHz  the  design 
frequency. 
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Figure  9.4-9.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.4-1  for  the  blind  design  (dotted)  and  the  variable  power 
divider  network  compensation  (solid)  at  306  MHz. 
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Table  9.4-3  Element  currents  used  to  produce  Figs.  9.4-4  through  9.4-9. 


294  MHz 

Element  No. 

Blind  Design 

Atten./Ph.  Shifter 

Power  Divider 

1 

0.1 2656/ 157.46° 

0.11161/189.7° 

0.10979/162.02° 

2 

0.25370/77.99° 

0.22646/109.57° 

0.26891/81.99° 

3 

0.32176/-22.37* 

0.36830/16.09° 

0.37959/- 1 2. 15° 

4 

0.5I418/-88.75® 

0.39I67/-65.89® 

0.50272/269.90“ 

5 

0.36823/168.92° 

0.38359/189.18° 

0.39586/174.04° 

6 

0.43171/107.68° 

0.30566/110.95° 

0.37092/91.76° 

7 

0.21874/-7.1 1° 

0.20139/13.16° 

0.23050/- 1.59° 

8 

0.19632/-51.160 

0.088 12/-65.900 

0.1 2030/- 76.34° 

300  MHz 

Element  No. 

Blind  Design 

Atten./Ph.  Shifter 

Power  Divider 

1 

0.13054/137.22° 

0.10436/165.00° 

0.11712/143.35° 

2 

0.23044/50.02° 

0.20646/75.00° 

0.23171/53.35° 

3 

0.3 11 10/ -45.47° 

0.323 15/- 15.00“ 

0.36269/-36.65° 

4 

0.47433/236.12° 

0.39799/255.00“ 

0.44671/233.35° 

5 

0.39349/140.01° 

0.39799/165.00° 

0.44669/143.35° 

6 

0.43402/69.29° 

0.32315/75.00° 

0.36271/53.35° 

7 

0.20806/-40.61* 

0.20646/- 1 5.00° 

0.23173/-36.65" 

8 

0.19572/260.94° 

0.10436/255.00“ 

0.11713/233.35° 

306  MHz 

Element  No. 

Blind  Design 

Atten./Ph.  Shifter 

Power  Divider 

1 

0.14284/105.88° 

0.10290/138.08° 

0.13260/111.71° 

2 

0.21113/18.36° 

0.18571/40.18° 

0.20781/21.75° 

3 

0.33308/-78.09® 

0.31038/-45.030 

0.38278/-69.470 

4 

0.44560/201.40° 

0.37339/216.76° 

0.41316/199.86° 

5 

0.40143/107.98“ 

0.41354/134.50° 

0.45570/109.74° 

6 

0.40001/33.83° 

0.30791/39.11° 

0.33384/19.68° 

7 

0.19011/-73.65* 

0.20078/-49.51° 

0.20793/-70.22® 

8 

0.17099/220.02° 

0.10378/214.10° 

0.09783/192.50° 
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9.4.4  Power  and  efficiency 


Table  9.4-4  summarizes  the  power  distribution  and  the  efficiency  of  the 
attenuator/phase  shifter  compensation  and  the  power  divider  compensation  for  the 
dipole  array  considered.  The  notation  as  defined  by  Smith  [S-2]  is  used  in  the  table. 

—  incident  power  produced  by  the  generator 
=  G2 

where  G  is  the  generator  voltage  magnitude. 

Pref  —  reflected  power  back  into  the  generator 
=  (SnG)2 

where  SN  is  the  reflection  coefficient  of  the  compensated  array  looking  into  the 
input  port. 

Pdi*  —  power  dissipated  by  the  compensation  network 

Prad  =  power  radiated 
*"  Pine  — '  Prtf  —  Pdi* 

The  efficiency  is  defined  as 

EfT.  -  x  100% 

*  inc 

The  power  distributions  are  calculated  from  the  results  of  MCAP  analysis  of  the 
compensated  arrays. 
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Table  9.4-4.  The  power  distribution  in  the  attenuator  phase  shifter  compensation 
network  and  power  divider  compensation  network  for  a  30  dB 
Dolph-Chebyshev  dipole  array  (all  powers  in  Watts). 


Atten./Ph.  Shifter 

Power  Divider 

P^c 

6.31335 

5.01155 

P«r 

0.00455 

0.00590 

P*. 

1.30290 

0.00000 

P^ 

5.00590 

5.00590 

Efficiency 

79.3  % 

99.8  % 
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In  Table  9.4-4  the  power  divider  networks  are  assumed  to  be  lossless.  Thus,  no 
power  is  lost  in  the  blind  design  feed  network  and  the  power  divider  compensated  feed 
network.  In  other  words,  Pdb  only  represents  the  power  lost  in  the  attenuators.  The  high 
efficiency  shown  in  Table  9.4-4  for  the  power  divider  compensation  results  from  this. 
When  line  losses  due  to  finite  conductivity  and  other  nonideal  effects  are  considered  in 
the  S-matrix  calculation  of  the  compensation  network,  the  power  divider  network  would 
be  somewhat  less  efficient. 

9.4.5  Computational  efficiency 

Figure  9.4-10  shows  plots  of  CPU  time  and  number  of  iterations  required  by  SANE 
and  SANE-PODCON  as  functions  of  interelement  phasing,  a,  in  desired  currents.  The 
calculations  are  based  on  the  array  of  eight  dipoles  shown  in  Fig.  9.4-1.  The  magnitude 
of  desired  currents  are  the  same  as  (9.4-1)  while  the  interelement  phasing  is  changed 
from  a  =■  0°  which  produces  a  broadside  pattern  to  a  —  1 80°  which  scans  the  main  beam 
to  endfire  direction. 

In  the  case  of  SANE  for  attenuator/phase  shifter  compensation,  a  blind  design  feed 
that  corresponds  to  each  interelement  phasing  is  used  in  the  calculations.  In  the  power 
divider  compensation  calculations  the  initial  guess  produced  by  the  program 
SANE-PODCON  is  used.  The  CPU  time  data  points  are  calculated  from  the  execution 
of  SANE  and  SANE-PODCON  on  IBM  3090  at  Va  Tech.  The  number  of  iterations  are 
calculated  for  the  error  in  the  norm  of  the  nonlinear  equations  to  be  less  than  1.0E-05. 

As  it  is  shown  in  Fig  9.4-10,  SANE  requires  larger  number  of  iterations  to  converge 
to  a  solution  than  SANE-PODCON.  This  is  generally  true  because  in  SANE  generator 
phase  is  changed  to  optimize  power  efficiency  of  auenuator/phase  shifter  compensation 
network.  In  thi«  particular  example,  generator  phase  was  changed  from  0°  to  170°  at 
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Figure  9.4-10.  CPU  time  and  number  of  iterations  required  by  SANE  and 
SANE-PODCON  to  converge  to  a  solution  as  functions  of 
interelement  phasing. 
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10°  steps.  Thus,  in  average  SANE  required  only  one  eighteenth  of  the  total  iterations  to 
converge  to  a  solution  at  each  generator  phase.  Also,  it  should  be  noted  that  though 
SANE  requires  more  iterations  than  SANE-PODCON,  the  required  CPU  time  for  this 
particular  array  are  in  the  same  order  of  magnitude.  SANE  requires  less  CPU  time  per 
iteration  than  SANE*PODCON  because  the  form  of  nonlinear  equations  is  simpler  for 
attenuator/phase  shifter  compensation. 

It  is  also  apparent  in  the  plots  that  SANE-PODCON  becomes  more  difficult  to 
converge  as  the  larger  interelement  phasing  is  required.  For  example,  the  required  CPU 
time  for  a  ■*  180°  (6.32  seconds)  is  about  5  times  longer  than  for  the  broadside  pattern 
a  =  0°  (31.84  seconds).  On  the  other  hand,  the  change  in  required  CPU  time  for  SANE 
is  small  for  different  case  of  interelement  phasing. 

The  results  presented  here  are  only  for  the  8-element  dipole  array  of  Fig.  9.4-1. 
For  an  array  with  larger  number  of  elements,  the  convergence  by  SANE-PODCON 
would  be  slower  than  by  SANE. 

9.4.6  Conclusions  from  the  dipole  array  study 

In  this  study  the  latest  wire  analysis  techniques  were  used  in  combination  with  the 
proposed  mutual  coupling  synthesis  methods.  This  was  done  for  both  attenuator/phase 
shifter  compensation  and  power  divider  network  compensation.  The  center  frequency 
radiation  patterns  for  the  30  dB  eight  element  parallel  dipole  array  are  shown  in  Figs. 
9.4-5  and  9.4-8  for  the  two  compensation  networks.  Comparison  of  the  blind  design  and 
compensated  array  patterns  in  both  compensation  network  cases  show  that  in  the 
H-plane  the  high  side  lobe  level  (-15  dB)  was  reduced  to  the  design  value  of  -30  dB.  In 
the  E-plane  side  lobe  level  of  -15  dB  was  reduced  to  -24  dB  through  compensation.  The 
desired  side  lobe  reduction  was  not  quite  achieved  in  the  E-plane  due  to  the  varying 
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shape  of  individual  element  pattern  for  dipole  antennas  in  this  elevation  plane,  whereas 
the  azimuth  plane  (H-plane)  has  omnidirectional  patterns  for  all  elements. 

9.5  An  Eight  Element  Monopole  Array  Above  a  Finite  Ground  Plane 

In  this  section  we  will  consider  synthesis  of  a  30  dB  Dolph-Chebyshev  pattern  using 
an  array  of  eight  monopoles  over  a  finite  ground  plane.  The  z-directed  monopoles  have 
length  L  =  0.233ii„,  radius  a  =*  0.005>l0  and  are  spaced  d  =  0.500AB  apart  along  x-axis. 
The  ground  plane  is  a  5.0>l0  x  5.(M0  rectangular  plate  of  a  perfect  conductor  located  on 
xy-plane  at  z  =  0.  The  geometry  of  the  array  is  shown  in  Fig.  9.5-1.  Y-parameters  for 
the  elements  are  obtained  using  ESP  and  then,  they  are  transformed  to  S-parameters 
using  the  steps  1  and  2  discussed  in  Section  9.3. 

The  desired  pattern  is  a  30  dB  Dolph-Chebyshev  pattern  with  the  main  beam 
steered  to  30°  off  broadside.  The  required  element  current  distribution  is  the  same  as  the 
dipole  array  case  in  the  previous  section  and  is  shown  in  (9.4-1).  The  same  desired 
pattern  was  chosen  for  the  comparison  with  the  dipole  results. 

In  Section  9.2,  it  was  shown  that  ESP  may  have  difficulty  in  predicting  return  loss, 
S£,  for  monopoles  over  a  ground  plane.  For  the  array  considered  in  this  section,  the 
magnitude  of  S£  calculated  using  ESP  are  approximately  -5  dB.  These  are  unrealistically 
high  return  loss  which  makes  blind  design  radiation  pattern  worse  that  what  it  really 
should  be.  Thus,  for  the  blind  design  analysis  the  diagonal  elements  of  the  array  S-matrix 
are  replaced  with  the  S-parameters  for  the  monopoles  over  an  infinite  ground  plane 
calculated  using  MININEC  II.  The  element  currents  are  calculated  from  the  blind 
design  feed  network  used  in  Section  9.4  and  are  shown  in  Table  9.5-1.  The 


IX.  Numerical  Studies  with  Wire  Antennas 


184 


I 


corresponding  E-plane  and  H-plane  radiation  patterns  are  plotted  in  Fig.  9.5-2  with 
dotted  curves. 

The  required  compensation  network  parameters  are  computed  from  the  ESP 
calculated  element  S-parameters  and  the  blind  design  feed  network  S-parameters  using 
SANE  and  SANE-PODCON.  The  resulting  parameters  are  summarized  in  Tables  9.5-1 
and  9.5-2,  respectively.  Figure  9.5-2  shows  the  E-plane  and  H-plane  radiation  patterns 
of  the  compensated  array  (solid  curves).  In  the  H-plane,  the  side  lobe  level  of  -13  dB 
was  successfully  reduced  to  the  desired  level  of  -30  dB.  In  the  E-plane,  however,  the  side 
lobe  was  only  reduced  to  -15  dB.  Unlike  the  dipole  array  case  in  Section  9.4,  the  failure 
to  synthesize  the  desired  result  comes  not  only  from  the  shape  of  element  pattern  in  this 
plane,  but  also  from  the  variation  in  individual  element  pattern  introduced  by  the  finite 
ground  plane.  This  is  similar  to  result  observed  in  the  35  dB  Dolph-Chebyshev 
microstrip  array  as  discussed  in  Section  8.4. 
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Table  9.5-1.  Required  attenuators  and  phase  shifters  for  the  30  dB  Dolph-Chebyshev 
monopole  array  example. 


Attenuation 

Phase  Shift 

0.000  dB 

48.70° 

4.257  dB 

5.39° 

0.982  dB 

41.56° 

4.319  dB 

7.68° 

1.978  dB 

22.00“ 

5.469  dB 

3.48° 

2.735  dB 

20.08° 

5.380  dB 

-0.80° 

2.55183 

130.00° 
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Table  9.5-2.  Required  variable  power  dividers  and  phase  shifters  for  30dB  Dolph 
Chebyshev  monopole  array  example. 


Blind  Design 

Compensated  Design 

K 

Ph.  Shift 

K 

Ph.  Shift 

0.79649 

-90.00° 

0.07724 

-88.48° 

0.60267 

-90.00° 

0.06040 

-92.17° 

0.39773 

-90.00° 

0.03519 

-95.22° 

0.20351 

-90.00° 

0.01446 

-101.29° 

0.83083 

180.00* 

0.85802 

-168.26° 

0.16917 

180.00° 

0.16804 

-174.55° 

0.50000 

0.00° 

0.52595 

4.15° 

1.99857 

-135.00° 

1.46877 

137.73° 
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Figure  9.5-2.  Comparison  of  E-plane  and  H-plane  radiation  patterns  of  the  array 
of  Fig.  9.5-1  for  the  blind  design  (dotted)  and  compensated  design 
(solid)  at  300  MHz. 
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Table  9.5-3.  Element  currents  used  to  produce  Figure  9.5-2. 


Element  No. 

Blind  Design 

Atten./Ph.  Shifter 

Power  Divider 

1 

0.04932/189^96° 

0.10274/185.00° 

0.17853/177.27° 

2 

0.32404/85.85° 

0.20324/95.00“ 

0.35319/87.27“ 

3 

0.14857/-9.23* 

0.31812/5.00“ 

0.5528 1/-2.730 

4 

0.78694/261.32“ 

0.39180/-85.00’ 

0.68082/267.27° 

5 

0.27949/159.52° 

0.39180/185.00“ 

0.68085/177.27° 

6 

0.91869/83.60° 

0.31812/95.00“ 

0.55281/87.27° 

7 

0.15310/-41.37* 

0.20325/5.00° 

0.353 19/-2.73° 

8 

0.451 19/-88.380 

0.I0274/-85.000 

0.17854/267.27° 
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X.  Conclusions 


This  report  reflects  four  years  of  study  into  the  effects  of  mutual  coupling  on  array 
analysis  and  synthesis.  During  this  effort  several  fundamental  array  issues  were 
addressed  as  well  as  new  array  synthesis  techniques.  Progress  was  made  in  answering 
several  fundamental  questions,  but  much  remains.  The  objective  of  this  study  was  to 
develop  and  verify  synthesis  techniques  for  array  antennas  which  include  all  coupling 
effects.  This  was  accomplished  using  an  attenuator/phase  shifter  network  and  a  variable 
power  divider  network. 

In  this  chapter  the  key  results  are  collected  and  summarized.  References  are  made 
to  the  appropriate  points  in  the  report  for  further  details.  Recommendations  for  future 
study  are  also  given  in  this  chapter. 
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10.1  Summary  of  Research  Results 

Fundamental  Issues  In  Array  Theory 

1.  Scattering  Parameter  Representation  of  Arrays  with  Coupling  Effects  Included 

The  scattering  matrix  representation  for  general  arrays  was  developed  in  Sec.  2.1 
and  array  coupling  effects  in  terms  of  scattering  parameters  were  discussed  in  Sec.  2.2, 
including  examples  of  the  coupling  effects  on  radiation  patterns  (Figs.  2.2-3  and  2.2-4). 


2.  Array  Pattern  Analysis  Including  Mutual  Coupling 

Array  pattern  analysis  techniques  in  the  presence  of  mutual  coupling  was 
developed  (Sec.  2.3).  It  was  shown  that  radiation  pattern  of  an  array  driven  by  the 
feed/generator  network  can  be  approximated  by  (2.3-9)  as  follows: 


N 

oo.i-i) 

nm\ 

I 

The  element  currents,  I„,  are  obtained  using  the  S-parameter  network  analysis  technique 

discussed  in  Section  2.1.  gw,(^,  0)  is  the  pattern  of  each  element  in  the  array  environment 

when  all  other  elements  are  match  loaded.  When  element  patterns  are  similar  in  shape  j 

and  gain,  (10.1-1)  can  be  simplified  to  (2.3-10): 


F(M)*g £lBe*-  (10.1-2) 

nm\ 

where  g*(0, <£)  is  a  typical  element  pattern  of  the  array.  These  equations  need  further 
investigation;  see  recommendation  4  in  the  next  section. 
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Computer  Aided  Microwave  Network  Analysis  To  Assist  in  Array  Design 


3.  The  MCAP  program 

MCAP,  the  Microwave  Circuit  Analysis  Program,  originally  published  in  (G-l)  was 
modified  by  PSL  and  by  Smith  during  the  effort  for  microstrip  array  feed  analysis  and 
is  described  fully  in  Sec.  4.2. 

4.  Power  Divider  Networks 

MCAP  was  used  to  study  frequency  response  of  a  two-way  variable  power  divider 
(Sec.  4.4.3).  The  results  showed  that  for  a  large  power  ratio  the  phase  difference  between 
two  output  ports  of  a  quadrature  hybrid  power  divider  deviates  rapidly  from  zero  for 
operations  off*  of  the  design  frequency  (See  plot  in  Fig.  4.4-4).  The  two-way  power 
dividers  were  then  used  as  building  block  to  develop  a  N-way  variable  power  divider 
network  (Sec.  6.2). 

General  Array  Synthesis  Including  Mutual  Coupling  Compensation 

5.  Overview  of  Pattern  Synthesis 

It  is  assumed  that  the  principle  of  pattern  multiplication  applies  and  that  the 
pattern  of  an  array  is  factorable  in  the  form  of  (2.3-10)  as  follows: 

F(e.+)-gje,M*.*)  (10.1-3) 


where 


N 

f(0,  <f>)  =  normalized  array  factor  oc  y^lnej{"  (10.1-4) 

n»l 
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where  <f#  are  defined  in  (2.3-2).  For  desired  pattern  F(0,  <£)  and  a  known  active  element 
pattern  g we  can  solve  for  the  corresponding  desired  array  factor  using  (10.1-3) 
as 


W,*) 


m<t>) 

g*(M) 


(10.1-5) 


Then  conventional  'text  book"  synthesis  techniques  that  do  not  account  for  coupling 
are  used  to  find  the  desired  currents  {I„}.  The  realization  of  these  currents  in  the  presence 
of  coupling  effects  is  the  thrust  of  the  proposed  methods  (see  following  items).  The 
important  assumption  here  is  that  the  individual  active  element  patterns  are  sufficiently 
alike  that  (10.1-5)  holds. 


6.  Summary  of  System  of  Equations 

Figure  10.1-1  summarizes  the  solution  formulation.  A  system  of  nonlinear 
equation;  which  describes  a  general  array  antenna  network  consisting  of  antenna 
elements,  a  feed  network,  a  compensation  network  and  a  generator  was  developed  in 
(2.4-14)  and  (2.4-15)  and  are 


N 

M 

M 

o - a^  +  XsjfX 

cf/°  +  £s^Ui 

for  1  <,  k  <,  N 

(10.1-6) 

i-l 

i-l 

j-l 

0  =*  — Uk_N 


N 

V 


M 

Y"1 


+ + 2,s£  N+i 

i-l  i-l 


JIG 


forN+l^k^N  +  M 


(10.1-7) 
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Knowns:  Feed /Generator  Network  [SF/G] 

Array  Excitation  [aA]  [bA] 

Compensation  Network  Form  [Sc] 

Unknowns:  Reflected  Waves  from  Compensation  Network  {UJ 

Compensation  Network  Parameters  {X;}  where  Sc  =  F(X) 
Induced  Waves  from  Generator  {cf/G} 

Figure  10.1-1.  General  compensation  approach  for  a  N-element  array  with  a 
M-way  feed  network. 
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where  {a£ }  and  {bf }  are  obtained  from  (2.4-7)  and  (2.4-8)  using  desired  element  currents 
{If}  and  known  scattering  parameters  of  the  antenna  elements  [SA].  The  feed/generator 
network  has  combined  S-parameters  [SF/°]  and  the  generator  induces  waves  [cF/G]  at  the 
feed/generator  output  ports.  The  unknown  variables  {X,}  which  represent  the  physical 
specification  for  the  compensation  network  are  imbedded  in  the  scattering  parameters 
[S0].  The  variables  {U,}  are  intermediate  unknowns  necessary  to  expand  the  system  of 
equations  in  closed  form.  The  equations  are  nonlinear  because  of  the  terms  S^N*iU, 
which  are  products  of  the  unknown  variables.  In  the  case  where  the  generator  value  is 
treated  as  a  unknown,  the  products  S^cf'0  also  contribute  to  the  nonlinear  property 
of  the  equations. 

7.  Mutual  Coupling  Effect  in  Symmetric  Excitations 

Smith  [S-l,Sec  8.4]  showed  that  an  ideal  N-way  power  divider  could  be  used  to 
produce  a  sum  pattern  which  would  be  insensitive  to  the  presence  of  mutual  coupling. 
The  statement  can  be  extended  to  any  symmetrical  excitation.  That  is,  if  the  feed 
network  produces  an  array  excitation  which  is  symmetric  (amplitude  and  phase 
symmetric  about  center  of  the  array)  mutual  coupling  effects  (which  can  be 
unsymmetric)  do  not  alter  the  symmetry  of  the  resulting  currents.  Note  that  if  a  linear 
phase  taper  is  imposed  the  excitation  is  unsymmetric  and  mutual  coupling  efTects  alter 
the  current  distribution.  This  property,  however,  is  only  applicable  to  ideal  power 
dividers  which  have  S-parameters  as  derived  by  Smith  [S- 1, (4.4- 1 9)].  The  variable  power 
divider  network  derived  in  this  effort  (item  4)  has  S-parameters  different  from  the  N-way 
power  divider  derived  by  Smith.  Thus,  the  property  does  not  apply  to  the  variable  power 
divider  network. 

A  sum  pattern  feed  network  that  has  Smith's  form  of  S-parameters  can  be  realized 
using  50Q  coaxial  T-junctions  connected  in  a  corporate  feed  configuration.  Such  feed 
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network  would  be  least  effected  by  mutual  coupling.  However,  the  network  would  have 
high  VSWR  at  the  feed  output  ports  resulting  in  low  power  efficiency. 


8.  The  Decomposition  of  Element  Currents 

The  currents  in  each  element  of  an  array  can  be  decomposed  into  three  terms: 


I 


^ ideal  ^ antenna  coupling  ^full  coupling 


(10.1-8) 


This  was  derived  in  Sec.  5.4  for  an  attenuator/phase  shifter  compensation  two  element 
array  of  matched  antennas  but  should  hold  with  mismatches  present  as  well.  In  (10.1-8) 
IMwi  is  the  current  that  would  be  present  if  no  mutual  coupling  were  present  anywhere 

in  the  system.  I... _ is  the  contribution  to  the  current  in  one  element  induced  by 

the  other  elements  through  antenna  coupling.  This  current  exists  when  there  is  mutual 
coupling  between  antenna  elements  but  no  coupling  through  the  feed  network.  I*,,  _,.Tllnt 
is  present  when  there  is  a  coupling  between  output  ports  of  the  feed  network.  The 
coupling  through  the  feed  network  introduces  a  feedback  loop  in  the  system  which 
makes  the  currents  a  nonlinear  function  of  compensation  variables;  this  is  discussed  in 
item  9  below. 


9.  The  Nonlinear  Property  of  Attenuator/  Phase  Shifter  Compensated  Arrays 

The  nonlinear  relation  between  attenuator/phase  shifter  compensation  network 
parameters  and  element  currents  was  demonstrated  for  a  two  element  array  of  matched 
antennas  in  (5.4-10)  and  (5.4-11).  It  was  shown  that  the  compensation  network  can  be 
obtained  by  solving  a  system  of  linear  equations  when  only  antenna  coupling  is 
considered.  When  feed  network  coupling  is  present  the  relationship  between  the  element 
currents  and  the  compensation  network  parameters  becomes  a  nonlinear  function. 
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10.  The  Generator  Level  Setting  Effect 

In  solving  the  system  of  nonlinear  equations  for  attenuator/phase  shifter 
compensation,  it  was  shown  that  the  constraint  in  (5.2-6)  must  be  imposed  on  the 
unknown  variables  to  ensure  that  a  physically  realizable  network  results.  Smith  [S-l] 
showed  that  the  magnitude  of  generator  can  be  adjusted  to  meet  this  constraint.  The 
process,  however,  is  a  nonlinear  function  because  of  the  relation  between  element 
currents  and  generator  values.  It  was  shown  that  the  constraint  can  be  used  to  eliminate 
an  attenuator  from  the  attenuator/phase  shifter  compensation  network. 

Specific  Solutions 

11.  Attenuator j Phase  Shifter  Compensation  Network 

A  system  of  equations  for  the  synthesis  of  attenuator/phase  shifter  compensation 
networks  that  was  derived  by  Smith  [S-l]  was  summarized  in  Section  5.2.  In  addition, 
our  investigation  showed  that  multiple  solutions  exist  and  efficiency  can  be  improved 
by  changing  the  absolute  phase  of  the  generator;  see  Section  5.3.  The  computer 
program  SANE  [S-l]  was  modified  to  include  this  effect  in  the  solution  for  the  system 
of  equations.  A  description  of  changes  made  in  SANE  is  discussed  in  Appendix  A. 

12.  Power  Divider  Compensation  Network 

A  general  system  of  nonlinear  equations  for  the  synthesis  of  compensation 
networks  (item  6)  was  modified  for  power  divider  networks  so  that  they  can  be  solved 
numerically;  see  (6.3-6).  Also,  the  initial  guess  to  the  solution  of  the  nonlinear  equations 
which  plays  an  important  role  in  numerical  technique  was  discussed  (Sec.  6.4).  The 
computer  program  SANE-PODCON  was  written  to  solve  the  nonlinear  equations  using 
a  modified  Newton's  method;  see  Appendix  B.  The  code  converges  much  more 
efficiently  than  SANE-GECON  written  by  Smith  [S-l]. 
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13.  Transmission  Line  Compensation  Network 

Mutual  coupling  compensation  with  transmission  lines  was  attempted  (Chap.  7). 
However,  except  for  the  case  of  two  element  arrays,  convergence  to  a  solution  for  the 
nonlinear  equations  was  never  achieved.  This  failure  seems  to  come  from  the  physical 
nature  and  the  specific  structure  of  scattering  parameters  of  the  network  rather  than 
from  the  use  of  inadequate  numerical  techniques. 

Verification 

14.  Analysis  Experiments  with  a  Microstrip  Array 

The  S-parameter  analysis  technique  developed  in  Section  2. 1  and  the  array  pattern 
analysis  technique  developed  in  Section  2.3  were  applied  to  an  eight  element  microstrip 
array  in  Section  8.2  for  sum,  difference  and  endfire  excitations.  Comparison  of  measured 
and  calculated  patterns  (Figs.  8.2-2  through  8.2-4)  suggests  that  the  techniques  are 
basically  valid.  However,  close  match  of  measured  and  calculated  patterns  were  never 
obtained  due  to  uncertainty  in  accuracy  of  measured  S-parameters  and  variation  in 
active  element  patterns  of  the  microstrip  array. 

15.  Synthesis  Experiment  with  a  Microstrip  Array 

Attenuator/phase  shifter  networks  and  power  divider  networks  were  designed  in 
Section  8.3  to  compensated  for  mutual  coupling  in  the  eight  element  microstrip  array  for 
sum,  difference  and  endfire  excitation.  In  Section  8.4  an  attenuator/phase  shifter 
compensation  network  was  designed  to  synthesize  a  35  dB  Dolph-Chebyshev  pattern 
using  the  microstrip  elements  and  sum  feed  network.  The  compensation  network  was 
built  and  tested  at  PSL.  The  measured  patterns  (Figs.  8.4-1  and  8.4-3)  showed  reduction 
in  the  sidelobe  level  from  -17  dB  to  -21  dB  and  improvement  in  main  beam  shape.  The 
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desired  sidelobe  level,  however,  was  not  achieved  because  of  the  variation  in  active 
element  patterns. 

16.  Wire  Antenna  Computer  Studies 

A  30  dB  Dolph-Chebyshev  pattern  with  main  beam  steered  to  30  degrees  off 
broadside  was  synthesized  using  an  array  of  eight  dipoles  and  an  array  of  eight 
monopoles  over  a  finite  ground  plane.  The  dipole  and  monopoles  arrays  were  modeled 
using  the  ESP  moment  method  computer  code.  The  required  attenuator/ phase  shifter 
network  and  power  divider  network  parameters  are  given  in  Tables  9.4-1  and  9.4-2  for 
the  dipole  array  and  in  Tables  9.5-2  and  9.5-3  for  the  monopole  array.  The  calculated 
E-plane  and  H-plane  patterns  of  the  blind  design  and  and  compensated  dipole  array  are 
shown  in  Fig.  9.4-5  for  attenuator/phase  shifter  compensation  and  in  Fig.  9.4-8  for 
power  divider  compensation.  The  calculated  patterns  of  the  monopole  array  are  shown 
in  Fig.  9.5-2.  In  both  dipole  and  monopole  arrays,  complete  compensation  was  achieved 
in  H-plane.  However,  the  sidelobes  were  not  reduced  to  the  desired  level  in  E-plane  due 
to  varying  shape  of  individual  element  pattern  in  the  elevation  plane.  (See  discussion  in 
Sec  9.4.6). 

Comparison  of  Compensation  Methods 

17.  Accuracy  and  Convergence 

Attenuator/phase  shifter  and  power  divider  compensation  network  parameters  were 
successfully  calculated  for  eight  element  arrays  in  Chapters  8  and  9.  However,  solutions 
were  obtained  only  for  two  element  array  cases  using  transmission  line  characteristic 
impedance  network  technique  (See  Chapter  7). 

Smith  found  that  convergence  of  the  power  divider  network  equations  are  very 
sensitive  to  the  initial  guess  to  the  solution.  The  convergence  problem  was  solved  by 
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using  the  initial  network  that  would  be  necessary  without  considering  mutual  coupling 
(See  Section  6.4). 

18.  Frequency  Sensitivity 

Comparison  of  calculated  patterns  for  the  dipole  arrays  at  the  center  frequency  of 
300  MHz  and  side  frequencies  of  294  MHz  and  306  MHz  showed  that  the  compensation 
networks  are  not  extremely  sensitive  to  small  variations  in  the  operating  frequency.  See 
Figs.  9.4-4  through  9.4-9. 

19.  Computer  Time 

The  required  CPU  time  and  number  of  iterations  for  the  eight  element  dipole  array 
were  compared  for  different  cases  of  interelement  phasing  in  desired  currents  in  Section 
9.4.5.  It  was  shown  that  SANE  for  attenuator/phase  shifter  compensation  requires  more 
iterations  than  SANE-PODCON  for  power  divider  compensation  because  SANE 
optimizes  the  compensation  network  for  power  efficiency  by  varying  the  phase  of 
generator.  The  required  CPU  time,  however,  were  similar  for  both  SANE  and 
SANE-PODCON  for  the  particular  array  considered.  It  was  also  shown  that 
convergence  to  a  solution  by  SANE-PODCON  becomes  slow  as  the  interelement 
spacing  is  increased  (See  Fig.  9.4- 10). 

20.  Efficiency  of  Synthesis  Network 

Comparison  of  power  efficiency  of  attenuator/phase  shifter  network  and  power 
divider  network  is  presented  in  Table  9.4-5.  The  power  divider  network  was  shown  to 
have  very  high  efficiency. 
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10.2  Recommendations  for  Future  Work 


1.  Experimental  Verification 

The  next  major  step  necessary  in  this  study  is  experimental  verification  of  the 
technique.  If  the  present  scattering  parameter  models  for  the  array  can  be  realized  then 
other  types  of  compensation  network  can  be  developed.  For  example,  the  technique  may 
be  extended  to  allow  series  feeding  of  the  elements.  Otherwise  more  rigorous  model  for 
the  array,  feed  and  compensation  networks  must  be  developed. 

The  present  experiments  as  described  in  Chapter  8  use  measured  scattering 
parameters  to  compensated  for  mutual  coupling.  Such  measurements  are  very  difficult 
for  large  element  array.  The  Moment  methods  can  be  used  to  approximate  for  the 
mutual  coupling  among  antenna  elements  as  shown  in  Chapter  9.  This  needs  to  be 
verified  also  through  experiments. 

2.  Array  Factoring 

During  the  course  of  this  study  it  was  discovered  that  fundamental  questions  in 
basic  array  pattern  formulation  remain.  For  example,  a  thorough  study  on  (10.1-1), 
(10.1-2)  and  the  developments  in  Section  2.3  could  be  examined  with  wire  antenna 
moment  method  computer  experiments. 

3.  General  Study  of  Sufficient  Conditions  for  Compensation 

In  this  study  we  found  two  network  architectures  which  permit  compensation  of 
mutual  coupling  effects.  A  worthwhile  study  would  be  to  find  the  necessary  conditions 
for  compensation  in  general. 
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4.  Inclusion  of  Individual  Element  Patterns 

The  synthesis  methods  developed  here  work  well  only  if  array  factorization  as  in 
(10.1-2)  holds.  A  challenging  area  of  study  would  be  to  develop  synthesis  techniques 
with  different  individual  element  pattern  as  in  (10.1-1). 

5.  Termination  of  Synthesis  Codes 

The  termination  condition  of  the  iteration  in  the  computer  programs  SANE  and 
SANE-PODCON  need  to  be  modified.  The  present  programs  use  error  in  the  element 
currents  as  the  condition.  This  should  be  replaced  with  the  error  criterion  in  radiation 
pattern. 

6.  Extension  to  Active  Devices 

It  may  be  possible  to  apply  this  technique  to  use  active  component  in  the 
compensation  network.  This  will  lead  to  synthesis  technique  where  antenna  elements, 
the  feed  network  and  compensation  network  as  well  as  any  other  networks  such  as 
amplifiers  are  implemented  on  the  same  substrate. 

7.  Shaped  Beam  Synthesis 

The  results  presented  in  this  report  are  for  narrow  beam  (or  difference)  patterns. 
Synthesis  of  shaped  main  beams  is  also  very  important.  Shaped  beam  patterns  were 
successfully  synthesized  using  both  attenuator/phase  shifter  and  power  divider 
compensation  networks  for  a  16  element  dipole  array  calculated  using  ESP.  Experiments 
must  be  conducted  using  array  with  larger  number  of  elements. 
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XII.  Appendix  A:  Program  SANE,  Version  2.0 


This  chapter  presents  a  brief  description  of  the  modifications  to  and  a  user's  guide 
for  the  program  SANE  (Synthesis  of  Array  with  Network  and  Element  coupling) 
originally  written  by  Smith  [S-l].  SANE  is  designed  tc  solve  the  system  of  nonlinear 
equations  in  (5.2-5)  for  attenuator/phase  shifter  compensation  network  parameters  given 
the  desired  element  currents  and  the  S-parameters  of  the  feed  network  and  the  antenna 
elements.  SANE  can  also  be  used  instead  of  MCAP  to  analyze  the  effects  of  mutual 
coupling  on  element  currents  in  an  uncompensated  array  using  techniques  discussed  in 
Section  2.1.  A  review  of  the  theory  used  in  SANE  is  presented  in  Chapter  5. 


12.1  Description  of  The  Code 

Two  modifications  were  implemented  to  SANE  (forming  Version  2.0)  in  this  effort. 
The  first  modification  was  the  addition  of  a  subroutine  to  handle  feed  line  effects  as 
discussed  in  Section  2.5.  The  subroutine  PREPRC  calculates  the  S-matrix  of  antenna 
elements  given  the  S-matrix  of  a  feed  line  and  the  S-matrix  of  the  elements  that  include 
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feed  lines.  This  is  a  FORTRAN  implementation  of  equations  (2.5-19)  through  (2.5-22). 
Using  the  S-matrix  of  elements  without  feed  lines,  the  desired  currents  at  the  elements 
can  be  transformed  to  the  desired  currents  at  the  input  to  feed  lines  by  a  process  similar 
to  (2.5-3). 

The  second  change  to  SANE  was  implementation  of  generator  phase  variation  as 
discussed  in  Section  5.3.  This  is  coded  as  an  additional  DO  loop  in  the  main  program  to 
change  initial  phase  value  of  the  generator.  This  allows  users  to  find  an  attenuator/phase 
shifter  network  that  is  power  efficient.  This  DO  loop  can  be  replaced  in  future  with  an 
optimization  routine  to  find  the  compensation  network  with  the  maximum  power 
efficiency. 


12.2  User's  Guide  to  SANE,  Version  2.0 

The  modifications  to  SANE  were  implemented  in  such  way  that  the  program  can 
be  operated  with  the  same  set  of  input  variables  as  the  original  SANE.  However,  the 
input  file  format  has  been  changed  slightly  to  accept  additional  information  necessary 
in  SANE  2.0.  An  user  must  perform  two  tasks  before  SANE  can  be  executed.  The  first 
task  is  to  modify  the  parameter  statements  in  the  program.  The  parameter  statements 
are  of  the  form  PARAMETER(NE  =  n)  where  n  is  the  number  of  antenna  elements  in 
the  array  being  considered.  There  are  ten  parameter  statements  in  SANE,  Version  2.0. 
The  parameter  statements  can  be  changed  using  a  text  editor. 

The  second  task  is  to  prepare  an  input  data  file.  The  data  file  is  set  up  as  a  series 
of  alphanumeric  character  lines  which  are  read  through  I/O  unit  5.  The  definition  and 
format  of  the  variables  in  the  data  file  are  listed  in  Table  12.2-1.  There  are  two  blocks 
of  variables  required  in  SANE.  The  first  block  of  numbers  controls  how  SANE  runs. 


XII.  Appendix  A:  Program  SANE,  Version  2.0 


213 


Table  12.2-1.  Input  file  format  for  SANE. 


Variable 

Format 

Description 

FIND 

A3 

Find  Compensation  Network 

PNTITS 

A3 

Print  iterations 

ERROR 

E10.3 

Maximum  error  in  the  nonlinear  equations 

AERROR 

E10.3 

Maximum  attenuator  deviation,  c 

MAXIT 

13 

Maximum  iterations 

MAXTRY 

13 

Maximum  generator  tries 

CHGPHS 

A3 

Change  generator  phase 

NPHASE 

13 

Number  of  generator  phase  changes  t 

DPHASE 

F10.5 

Increment  in  generator  phase  t 

GAM  MAG 

2FI0.5 

Reflection  coefficient  of  generator 

G 

2F10.5 

Initial  generator  value 

SF(1,1) 

2F10.5 

Feed  network  S-parameter  Sfa 

SF(N+  1,N  +  1) 

2F10.5 

Feed  network  S-parameter  S£+1>N+1 

SA(1,1) 

2F10.5 

Element  S-parameter  S£, 

SA(N,N) 

2F10.5 

Element  S-parameter  Sft  N 

REFRNC 

A3 

Feed  lines  exist 

SFL(1,1) 

2F10.5 

Feed  line  S-parameter  Sft  * 

SFL(2,2) 

2F10.5 

Feed  line  S-parameter  S^  * 

COORD 

A5 

Polar  or  complex  format  for  the  desired  currents 

KD 

2F10.5 

Desired  current  at  element  l 

I(N) 

2F10.5 

Desired  current  at  element  N 

X(l) 

2F10.5 

Compensation  variable  X,  X 

X(N) 

2F10.5 

Compensation  variable  XN  $ 

t  not  necessary  if  CHGPHS  *  'NO' 
X  not  necessary  if  FIND  =*  'NO' 

*  not  necessary  if  REFRNC  *  'NO' 
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The  first  two  lines  are  flags  FIND  and  PNTITS.  FIND  is  'YES'  if  a  compensation 
network  is  to  be  found  and  'NO'  otherwise.  PNTITS  is  'YES'  if  a  summary  of  the  results 
at  each  iteration  is  to  be  printed.  The  third  line  contains  the  absolute  error  allowed  in 
norm  of  the  nonliner  equations  (5.2-5).  The  allowable  deviation  of  the  magnitude  of  the 
maximum  attenuator  value  from  unity,  e,  as  defined  in  (5.2-6)  is  specified  in  the  fourth 
line;  a  tipical  value  for  e  is  0.01  which  corresponds  to  the  minimum  attenuator  of 
approximately  0.1  dB.  The  fifth  and  sixth  lines  contain  the  maximum  number  of 
iterations  for  Newton's  method  and  determining  the  necessary  generator  magnitude, 
respectively.  The  seventh  line  is  a  flag  CHGPHS  which  is  set  to  'YES'  if  the  phase  of 
initial  generator  is  to  be  changed  to  find  the  most  power  efficient  network  as  discused 
in  Section  5.3.  The  following  two  lines  are  required  if  CHGPHS  is  YES'.  The  first  of 
these,  NPHASE,  contains  the  number  of  generator  phase  changes  required,  and  the 
second,  DPHASE,  contains  the  angle  in  degrees  where  generator  phase,  [G ,  will  be 
incremented;  typically  NPHASE  is  set  to  18  and  DPHASE  is  set  to  10.0. 

In  the  second  block  the  array  characteristic  is  specified.  The  first  two  lines  in  the 
second  block  are  the  complex  reflection  coefficient  of  the  generator  (SG)  and  the 
magnitude  of  the  generator  excitation  ( I  c*3 1 ).  The  following  (n  +  1)J  lines  contain  the 
S-parameters  of  the  feed  network,  S£.  The  next  nJ  lines  contain  the  S-parameters  of  the 
elements,  Sj)  .  Following  these  is  a  flag  REFRNC  which  is  'YES'  if  feed  line  effects  as 
discussed  in  Section  2.5  are  to  be  considered.  In  this  case,  S-parameters  of  a  feed  line 
are  specified  in  the  four  lines  following  REFRNC.  A  flag  COORD  is  specified  in  the 
next  line.  COORD  is  set  to  'POLAR'  if  the  required  element  currents  are  specified  in  the 
magnitude-phase  format.  COORD  is  'COMPLEX'  if  the  required  currents  are  specified 
in  the  real-imaginary  format.  The  required  element  currents  are  specified  in  the  n  lines 
following  the  flag  COORD.  IF  the  flag  FIND  is  'NO'  then  the  next  n  lines  contain  the 
value  of  the  compensation  network. 


XII.  Appendix  A:  Program  SANE,  Version  2.0 


215 


As  an  example  consider  the  eight  element  microstrip  array  with  the  sum  feed 
network  as  described  in  Chapter  8.  Figure  12.2-1  shows  a  part  of  the  input  data  file. 
The  corresponding  output  file  which  is  written  on  I/O  unit  6  is  shown  in  Figure  12.2-2. 
Note  that  these  values  corresponds  to  those  in  Table.  8.3-1. 
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Figure  12.2-2.  The  resulting  output  from  data  in  Fig.  12.2-1  (see  Table  8.3-1). 


2  SSSffSff 

WlllWIUUflUltfUl 

oepooooo 


fooooc 

ffSifl 

a  oooooc 

f  OOOOQ< 


£  oooooooo 

_  MJUIWUJUJUIUIUI 

j  Illliiji 

•  oooooooo 

\ . 

j  ^NM«IA4N« 

I  SSiHHi 

|  uuuuuuuu 


500000000 
♦♦♦♦♦♦♦♦ 
UlUWWWWWItf 

•3 

a  HNNN««OOU> 
MNHON<ftAN 

c  OOOOOOOO 

i  I  I  I  I  I  I  I  I 

i  oooooxno 

•  OOOOOOOO 

umiuuiiuuwiii 

9  Krow^ooio 

•  ?••••••? 


A  ISSSSSSS 

•  Cii.Ct.LWL 

I 

_  ouuuuuuu 


i  *s 

I  *2 

I  u 


•4fNKNf« 

iclssass! 

assssssa 

35iSin  18525 

N I MNMMNN 


NMONOOf 

JWH^p'OO'i 

NjNNNroN 


SNANIflNHH 

HHAAxOAAN 

NMNN^NAO 

^aK)N^f9>N 


111*1111 


fNNN^AAN 

^«K-*NC«rt 

hmammwnn 

*KtoiA*oino 

MONHHHOM 


I  I  I  I  I  1  I 


nnm#ia«n« 


XXXXXXXX 


XII.  Appendix  A:  Progrun  SANE,  Version  2.0 


2t9 


Figure  12.2-2.  (Continued) 
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Figure  12.2-2.  (Continued) 


J 

XIII.  Appendix  B:  Program  SANE-PODCON 


This  chapter  presents  a  brief  description  and  user's  guide  to  the  program 
SANE-PODCON  (Synthesis  of  Array  with  Network  and  Element  coupling  for  POwer 
Divider  Compensation  Network).  SANE-PODCON  solves  the  system  of  nonlinear 
equations  (6.3-6)  for  power  divider  compensation  network  parameters  given  the  desired 
element  currents  and  the  S-parameters  of  the  antenna  elements.  A  discussion  on  the 
theory  used  to  develop  SANE-PODCON  is  presented  in  Chapter  6. 

13.1  Description  of  the  Code 

The  program  SANE-PODCON  consists  of  five  primary  blocks:  I)  initialization, 
2)  set  up,  3)  solve,  4)  verify  and  5)  output,  as  shown  in  Fig.  13.1-1.  The  first  block 
initializes  a  scaler  variable  and  two  array  variables.  It  consists  of  three  subroutines: 
GETEPS,  CALCII,  and  CALCMX.  The  subroutine  GETEPS  calculates  machine 
epsilon  of  the  computer  running  SANE-PODCON.  The  machine  epsilon  is  a  real 
number  which  a  computer  can  distinguish  the  difference  between  1  and  1  +  e.  It  is  used 


Xlll.  Appendix  B:  Program  SANE-PODCON 


221 


J 


Figure  13.1-1.  Block  diagram  of  SANE-PODCON. 


1.  Initialize 

2.  Set  Up 

5.  I/O 

3.  Solve 

4.  Verify 

5.  I/O 
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in  calculation  of  Jacobian  of  the  nonlinear  equations.  The  subroutines  CALCII  and 
CALCMX  create  indicies  used  to  represent  N-way  power  dividers  which  are  similar  to 
index  system  as  defined  in  Section  6.2.2. 

The  second  block  reads  data  from  a  input  file  and  establishes  the  nonlinear 
equations.  The  block  contains  four  subroutines:  INPUT,  PREPRC,  GENAB,  and 
GENINX.  The  subroutine  INPUT  reads  data  from  an  input  file  through  I/O  unit  5. 
PREPRC  is  called  when  S-parameters  of  the  elements  contain  the  feed  lines.  The  process 
is  discussed  in  Section  2.5.  The  routine  GENAB  calculates  the  incident  and  reflected 
wave  variables  at  the  elements,  {a*}  and  {b*}  from  the  desired  element  currents  using 
(2.4-7)  and  (2.4-8).  The  initial  guess  to  the  solution  from  desired  element  currents  are 
generated  in  GENINX.  This  is  discussed  in  Section  6.4. 

The  third  block  in  SANE-PODCON  solves  the  nonlinear  equations  set  up  by  the 
previous  block.  The  block  is  consists  of  NEWTON,  ALPHA,  CAIJAC,  CALCSC  and 
GETSCI.  The  subroutine  NEWTON  is  the  driver  routine  that  approximates  the 
solutions  to  the  nonlinear  equations  using  Damped  Newton's  method  [D-4].  The 
damping  factor  a  is  calculated  in  subroutine  ALPHA.  The  subroutine  CALJAC 
calculates  the  Jacobian  of  the  nonlinear  equations  at  a  particular  domain  represented 
by  {X„}.  It  uses  a  finite  difference  method  to  approximate  for  the  Jacobian  {D-4j.  The 
routine  calls  CALCSC  and  GETSCI  which  calculates  the  S-matrix  of  the  compensation 
network  using  equation  (6.2-6).  In  particular,  CALCSC  calculates  the  entire  S-matrix 
and  GETSCI  calculates  i-th  column  of  the  S-matrix.  This  was  done  to  increase  the 
computational  efficiency  since  the  compensation  network  S-matrix  has  many  redundant 
variables. 

The  forth  block  calculates  the  element  currents  with  the  compensation  network 
incorporated  into  the  array  network.  This  provides  verification  that  solutions  are  correct 
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and  are  not  affected  by  round  off  error.  The  element  currents  calculations  are  performed 
by  GETCUR. 


The  final  block  contains  two  routines:  ECHO  and  DISPLY.  The  subroutine  ECHO 
outputs  the  input  data  and  other  parameters  calculated  prior  to  solving  the  nonlinear 
equations.  DISPLY  on  the  other  hand  writes  the  solution  of  the  nonlinear  equations. 
Both  ECHO  and  DISPLY  writes  output  to  I/O  unit  6.  The  subroutine  DISPLY  also 
writes  the  resulting  S-matrix  of  compensated  feed  network  through  I/O  unit  7  in  a 
format  compatible  with  MCAP. 


13.2  User's  Guide  to  SANE-PODCON 


An  user  must  perform  two  tasks  before  SANE-PODCON  can  be  executed.  The 
first  task  is  to  modify  the  parameter  statements  in  the  program.  The  parameter 
statements  are  of  the  form  PARAMETER(NE=n)  where  n  is  the  number  of  antenna 
elements  in  the  array  being  considered.  There  are  eighteen  parameter  statements  in 
SANE-PODCON.  The  parameter  statements  can  changed  using  any  text  editor. 

The  second  task  is  to  prepare  an  input  data  file.  The  data  file  is  set  up  as  a  series 
of  alphanumeric  character  lines  which  are  read  through  I/O  unit  5.  The  definition  and 
format  of  the  variables  in  data  file  are  listed  in  Table  13.2-1.  The  first  line  in  the  data 
file  is  flag  PRTITR,  which  is  set  to  'YES'  if  a  summary  of  the  results  at  each  iteration 
is  to  be  printed.  The  second  line  contains  the  maximum  number  of  iterations  for 
Newton's  method.  The  third  line  contains  the  absolute  error  allowed  in  norm  of  the 
nonliner  equations  (6.3-6);  a  typical  value  of  ERROR  is  1.0E-05.  The  fourth  line  is  a 
flag  GENINT  which  is  set  to  'YES'  if  the  initial  guess  to  the  solution  vector  is  to  be 
generated  by  the  program.  If  GENINT  is  'NO'  user  must  supply  the  initial  guess  in  the 
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Table  12.2-1.  Input  file  format  for  SANE-PODCON. 


Variable 

Format 

Description 

PRTITR 

A3 

Print  iterations 

MAXITR 

14 

Maximum  iterations 

ERROR 

F8.5 

Maximum  error 

GENINT 

A3 

Generate  initial  guess 

X(l) 

2F10.5 

Initial  guess  for  power  divider  network  1  f 

X(N*1) 

2F10.5 

Initial  guess  for  power  divider  network  N-l  t 

G 

2F10.5 

Initial  generator  value 

GAM  MAG 

2F10.5 

Reflection  coefficient  of  generator 

REFRNC 

A3 

Feed  lines  exist 

SFL(1,1) 

2F10.5 

Feed  line  S-parameter  Sft  t 

SFL(2,2) 

2F10.5 

Feed  line  S-parameter  Sgr  t 

SA(1,1) 

2F10.5 

Element  S-parameter  S£, 

SA(N,N) 

2F10.5 

Element  S-parameter  Sft^ 

COORD 

AS 

Polar  or  complex  format  for  the  desired  currents 

1(1) 

2F10.5 

Desired  current  at  element  1 

KN) 

2F10.5 

Desired  current  at  element  N 

t  not  necessary  if  GENINT  *  'YES' 
$  not  necessary  if  REFRNC  *  'NO' 
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f 


following  n-1  lines.  The  next  two  lines  are  the  reflection  coefficient  of  the  generator 
(Sc)  and  the  initial  magnitude  of  generator  excitation  (|cG|).  The  following  nJ  lines 
contain  the  S-parameters  of  the  elements,  S$.  A  flag  REFRNC  follows  next.  REFRNC 
is  'YES'  if  feed  line  effects  as  discussed  in  Section  2.5  are  to  be  considered.  In  this  case, 
the  four  lines  following  REFRNC  are  the  S-parameters  of  a  feed  line.  A  flag  COORD 
defines  the  format  of  the  desired  currents  specified  next.  COORD  is  'POLAR'  if  the 
required  element  currents  are  specified  in  magnitude-phase  format.  COORD  is 
'COMPLEX'  if  the  required  currents  are  specified  in  real-imaginary  format.  The  desired 
element  currents  are  specified  in  n  lines  following  the  flag  COORD. 

As  an  example  consider  the  eight  element  microstrip  array  for  the  sum  pattern 
current  distribution  as  described  in  Chapter  8.  Figure  13.2-1  shows  a  part  of  the  input 
data  file.  The  corresponding  output  file  which  is  written  on  I/O  unit  6  is  shown  in  Figure 
13.2-2.  Note  that  these  values  corresponds  to  those  in  Table.  8.3-2. 
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NO 

0150 

00.00010 


print  iterations 
nsx  iterations 
nax  error 


YES 

0.00000 

0.00000 

generate  ir 
Generator 

0.00000 

0.00000 

6a— a  Q 

YES 

feed  lines 

0.45757 

-0.05515 

fl  1,1 

-0.04479 

0.88415 

fl  2,1 

-0.04479 

0.88415 

fl  1,2 

0.44725 

0.10151 

fl  2,2 

-0.07940 

0.12100 

array  1,1 

0.20500 

-0.24200 

array  2,1 

-0.15800 

0.01110 

array  5,1 

0.07190 

0.08280 

array  4,1 

0.00784 

-0.08980 

array  5,1 

-0.04110 

0.05250 

array  4,1 

0.04200 

0.02940 

array  7,1 

0.00045 

-0.05710 

array  8,1 

0.00045 

-0.05710 

array 

1,8 

0.05750 

0.05240 

array 

2,8 

-0.05490 

0.01750 

array 

5,8 

0.01790 

-0.04480 

array 

4,8 

0.08210 

0.04890 

array 

5,8 

-0.14000 

0.02750 

array 

4,8 

0.14700 

-0.24700 

array 

7,8 

-0.04270 

0.14400 

array 

8,8 

AAR 

current  f< 

1.00000 

0.00000 

XI 

1.00000 

0.00000 

12 

1.00000 

0.00000 

15 

1.00000 

0.00000 

14 

1.00000 

0.00000 

15 

1.00000 

0.00000 

14 

1.00000 

0.00000 

17 

1.00000 

0.00000 

18 

Figure  13.2-1.  SANE-PODCON  input  file  for  the  eight  element  microstrip  array 
with  sum  pattern  excitation  (see  Fig.  8.1-1). 
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Figure  13.2-2.  (Continued) 


13-3  Program  Listing 


om»«m»»m»»»m*mmm*m*mmmm»mmmmmmsm*mmm 

c* 

c*  SANE  -  POOCON 

c* 

c*  Synthesis  of  Arrays  with  Network  and  El ament  coupling 

c*  for  PQwer  Oivider  Compensation  Networks 

cm 

c*  Programmer  :  Koichiro  Takaaixawa 

cm  Revision  :  August  4,  1968 

cm 
cm 

cmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm 
cm  1.3.6  6/04/86  aetcur  added 

cm  1.3.S  4/22/86  fixed  error  in  ealese  and  gatsci 
cm  1.3.4  10/25/87  fixed  error  in  geninx 
cm  1.3.5  10/14/67  fixed  error  in  geninx 
cm  1.3.2  10/11/87  fixed  bugs  in  ealese 

cm  1.3.1  10/10/87  sibrautine  getsci  added  to  increase  the  efficiency 
cm  1.3.0  10/09/87  subroutine  ealese  added  to  replace  function  sc 
cm  1.2.2  8/08/87  subroutine  prepro  added  to  taka  care  of  x -mission 

cm  lines  on  the  array  board 

cm  1.2.1  7/01/87  geninx  added  (creates  initial  guess  solution  vector! 

Cm  l.l.o  5/15/87  SANE -POOCON  written  from  SANE-GECON 
Cm 

Cmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm 


program  sane  oodcon 
parameter  INE*SI 

complexmlA  x(NE  ) ,  iacobf  NE  ,NE  ),f(NE ) 

complexml6  stli  2,2) ,sa2(  NE  ,NE  )  ,eurr2(  NE  ) ,sa( NE  ,NE  )  .curmtl  NE  ) 
real  error 
integer  maxitr 

character»3  yes,prtitr,genint , ref  me 
characterm26  vers 

common/inknwn/x ,  jacob 
common/ vr* ion/vers 

common/opt ion/error >  maxitr ,prt i tr ,gsnint 
common/pa  tch/st  1 ,  sa2  ,curr2 ,  ref  me 
common/antena/sa  ,curmt 
cowmon/tima/itimel,itime2 


call  time on 
e 

c***"vtps  i ofi  nunber 

O  12345678901234567890123456 

vers  ■  '1.3.6  August  4,  1988' 

yes»'YES' 

Q 

c — —initialization 
c 

call  maceps 
call  calcii 
call  calcmx 
c 

c— read  input  variables 

call  input 
o 

c— — get  antenna  patch  S-metrix 
c 

if  ( ref  me.  aq.  yes)  than 
call  prspre 
end  if 
e 

c - set  up  the  nonlinear  equations 

c 

call  genab 
c 

c- - generate  initial  guess 

e 

if  (gsnint.eq.yes )  than 
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000  OOO  OOO  OOO  OOO  OOO  00000  000  000  ooo 


call  ganinx 
and  if 

- print  the  input  variable*  and  tha  nonlinear  equations 

call  echo 

call  tiaeckt iti— 1  l 
—-solve  the  equetian* 
call  newton 

- print  the  result 

call  disply 

— -vari fy  tha  eoapansation  network 
call  getcur 


subroutine  input 

read  input  data  thru  i/o  unit  S 


parameter  (NE*S) 

coaplex*16  x( ne ) ,  jacabl  ne >ne  I  .sal  ne >n* )  ,curmt ( no ) 
cok1*x»U  g.geem g.saZt  ne  »ne  >  >currZI  ne )  »stl(  Z  >  Z  )  >sn  >cn 
real  error 

real**  d2r ,  wag,  ang, snag. ephasa, dang 
integer  naxitr 

character* 3  yas.prtitr  .ganint  .ref  me 
character*5  polar .coop! .coord 
so—aon/unknwn/x , jacob 
coaaaon/genrtr/g  ,gawag 
co— on/an  tena/sa  ,curmt 
coaaaon/netwrk/sn.cn 

cow on/opt  ion/arror  >— xi  tr  .prtitr ,  ganint 
co— on/pa  tch/s  tl  >s*Z  ,currZ  »ref  me 

dZr  a  3 . 141592654  /  160.0 
yes* 'YES' 
polar  *  'POLAR' 
coapl  a  'RECTA' 


- read  in  the  options 


read!5»10)  prtitr,— xitr, error 
10  for— tla3/i4/f6.5> 

- initial  guess 

readl5>30)  ganint 
SO  for— tl  *3 1 

if  I ganint. ne. yes )  then 
raadl 5.40  )  <x(  i  ),i*l.NE-l> 
for— tlZflO.5) 
and  if 

irator 


read!  5.401  g 
raadl 5.40 )  gaanag 
snaqaweaq 
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OOO  OOO  000  ooo  oio  oooolo  ooo  OOO  000 


cn«(  1.0,0.01 
x(NE )  «  g 

— antenna  food  lino* 

rood!  5  *50)  rofmc 
50  foroaf ( aS ) 

if  I  ref  me.  aq.  yes)  than 

road  (5,40 1  I (till i, j  I, j«l,2 ),i«l>2 I 
and  if 

- — antorma  • -matrix 

roadt  5 >40 1  ( (sat i,j  l»i*l,NE  ),i«l,NE  ) 

- das i rod  currant 

read!  5,60)  coord 
40  format!  a5) 

if  I  coord. aq. polar )  than 
do  100  i«l,NE 

raadl 5,40  )  nag ,ang 

curmtl  i  )  *  dcnplxl  aag«dcoil  ang*d2r )  .aagodainf  ang*d2r ) ) 
100  continue 

also 

road 15,40)  (curmtl  i  )»i»l,N£  I 
and  if 

- calculto  tha  initial  guaos  for  tha  Sonora tor 

if  IcdabsIxINE  I  ).la.l.0d-05)  than 
ewag  a  edabs I  curmtl  1 1 1 
ephaaa  «  dblol  phase!  curmtl  11)1*  d2r 
dang  ■  3. 141592654/4. 0» leva 1 -ephaaa 
x( no  )>dcaplxl caag*dcoa( dang )  ,caag*dcos( dang ) ) 
andif 
return 
and 


jfl 


I 


subroutine  prop re 

calculate  Smetrix  of  alaaants  without  faad  lines 
saa  Section  2.5 


paranetar  (NE*0) 

coaplax*14  sal  no, no  ) , sap!  no, no  ( , curmtl  no  I 
coaplax*16  curr2(na),sa2(ne,na  )>stl(2,2),ono 
eharactar*3  ref  me 
coamon/antona/sa  ,curmt 
conmon/pa  tch/s  tl  >sa2  ,curr2 ,  rofmc 

aneall. 0,0.0) 


—•taka  a  copy  of  Antom  S' 

do  100  i>l,na 
do  110  J*1 >ne 

sa2(  iti)>sa(  i»j) 
110  continue 
100  continue 

- iterate  from  i»l  to  HE 

do  200  i»l »na 
— calcualta  S*(i,i) 


itrix 


I 


sap(iti)  a  Isa2(  i.i  )-stl(  1,1 ) )/ 

<  I stll  2,2  )•(  sa2(  i  ,i  )-stl(  1,1 )  l«stl(  1,2  )«stl(  2,1 ) ) 
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000  0^000  OgO  000  000  ooo 


— -calculate  S’(i,j)  and  S'l  j»i)  for  j»l  to  NE  oxcopt  j«i 


210 


do  210  j*l»no 
if  ( j.no.i  1  than 
Mp(i>}  I<u2(iii  )•( 
tap!  jit  >«m2I  j , i  )*l 
and  if 
oontinua 


til  2 >2 )«*apl i»i  )  )/*tl( 1,2 ) 
til  2,2 )*sapl i,i ) )/atll 2,1 ) 


— ealoulata  S'ljtk) 

do  220  j*l,na 
do  221  k*l,na 

if  1 1  j.rm.  i  J.and.  (k.na.  i  ) )  than 
sapi j,k  >**a2(  j  ,k  1— stl(  2,2  )**ap(  j,i  )**ap(  i,k  1/ 
t  tona-stlt 2,2)*sapl  i,i ) I 

andif 

221  continua 

- kaap  a  copy  of  S' 

220  oontinua 

do  230  ]*l,nt 
do  ZSl  k*l,na 

sa2t  j,k)**ap<  j  ,k  I 
231  continua 

230  oontinua 


200 


itii 


ratum 

and 


subroutina  ganinx 

ganarata  initial  guaaa  to  tha  solution 


paraaatar  INE«8) 

eos^i*x»W  xINE),  jacoblNE  ),salNE,NE  hcurmtINE),  cdzaro 
raal*8  caagl NE  )  .cpnasal NE )  >paagl NE  )  ,pphasal NE  )  ,zaro , ona ,dang 
intagar  laval.ct 
co— on  /unkrwn/x ,  jacob 
co— on  /antana/sa,eurmt 

zaro*0.0 
ona* 1.0 

cdzaro  a  10.0,0.0) 

—ealculata  tha  aagnituda  and  tha  phasa  of  raquirad  currant 

do  100  i«l,NE 

caagl  i  )*cdabKleurmtl  i ) ) 

cphasal  i  )*dbla<  phasa ( curmtf  i ) )  )*3 . 141592654/180 . 0 
100  oontinua 
e 

c - -ealculata  tha  nuaber  of  lavala  of  peuar  dividars  nacassary 

o 

laval  a  inti  log!  float! NE  )  )/log(  2.0)+0.5) 
et»o 

do  200  i«l, laval 
jand*NE/2«»i 
do  210  jal,jand 
o 

o - ooasxjta  tha  aayiituda  and  phasa  of  initial  guass 

o 

kcoas>l*(  j-1  )«2**i*l 
kcoap2»j*2«ai-2**l  i-l)+l 
if  1 caagl kcoap2 ) . na . zaro )  than 
paagl  j+ct  )*ona/(  one+l  caagl  kcoapl  )/c—  g(  kcoap2 )  >**2  ) 
c—gl  kco—1  lacaagl  kcoapl  )/dsqrtl  ona  p— gl  j*ct ) ) 
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ooo  oto  oooto  H  to  ooooto  000 


•1m 


p— g(  jfet  I*  zero 
endif 

pphasel  jrct  )*cpha«a<  kcoap2  l-epha*e<  Kconpl ) 

210  continue 

el<eUi«nd 
200  continue 

- convert  the  initial  guasa  in  the  rectangular  coordinate  fora 

do  600  i-l,ne-l 

xl  i  )*dcaplx(  paagl  i  )«dco* I pphasel  i  )  )rp<*agl  i  )*dsinf  pphasel  ill) 
600  continue 
return 
and 


subroutine  eecap* 
coapute  machine  apailon 


raal  acheps»  sqrep* 
coaaon/aac par/  echaps»sqraps 
acheps  a  1.0 

10  acheps  a  achapa  /  2.0 

if  (nchepsal.O.gt.1.0)  goto  100 
achapa  *  achapa  •  2.0 
sqreps  a  sqrtl  acheps) 
return 
and 


subroutine  ganab 

generate  a'a  and  b'a  of  antenna  array 
using  atpa.  (2.4-7)  and  12.4-01 


pa  rams  tar  (HE  *6 1 

coaplax*16  al  NE  )  >b(NE  )  >aal  NE  >NE  )  ,curmt(  ME  )  ,g»garaaag 
complax*16  tanpvl  NE  I  >  tenpml  NE  »N£  I 

co*pl*x*16  sa21NE,NE  I  ,curr2(NE  I  ,stll  2,2  I  ,ap(  NE  )  >bp(NE  I 
character*!  ref  me,  yes 

conaon/wavat/a >b 
comaon/genr  tr/g  >  gamsag 
common/antena/sa  icurmt 
co«non/patch/stl,sa2,curr2,refmc 

yes* 'YES' 

luaaation  of  ( sa*currant  I 

do  200  iai,NE 

taapvl  i  )*(0. 0,0.0) 
do  100  j*l,NE 

if  I  ref  me.  aq.  yes)  than 

taapvl  i  )»ta*g»v(  i  )«sa2(  i,j  l*curmtl  j  ) 
also 

ta^rvl  i  I* tanpvl  i  )*sal  i>  j  )*curmt(  j  ) 
andif 
100  continue 
200  continue 

c - 11-SAl 

c 

do  400  i*l,NE 
do  300  3*1 ,NE 

if  (rafme.oq.yaa)  than 


XIII.  Appendix  B:  Program  SANE-PODCON 


235 


{ 


30 0 

400 

e 

c - 

c 

e 

c-— 

o 


it  i ,  j  )*-sa2t  i»  j  > 
i,j)»-sal i>}) 

i(  i>i  )♦(  1.0  tO.O  ) 


500 


t, 

•1m 

t(_ 

and  if 
continue 
tempmt 
continue 


tempm**(  -1  )*t« 

call  matsol  1  ME » tempv , tempi ) 

•^current  »b 

do  500  i*l,NE 
bl  i  )“ tempv ( i  ) 
at  i  )«curr«t(  i  )*bt  i  I 
if  (  refmc.eq.ye* )  than 
apt  i  )*•(  i  ) 
bpf  i  )*bt  i  ) 

at  i  )*( apt  i  )-stl(  2,2  )*bpl  i  )  )/stll  2,1 ) 
bl  i  l*stl(  1,1  )*a(  i  )  tat  It  1,2  )*bp(  i  ) 
curr2t  i  )*curmtt  i ) 
curmtl  i  )«•<  i  )-bt  i  ) 
endif 
continua 
return 


a 

c* 

c 

e 

o 

e 

c* 

c 


(dtroutine  acKo 
echo  tha  input  file 


parameter  INE“0) 

conpl«x*16  xt  NE ) ,  jacobt  ME  ,NE  )  ,at  ME  )  ,bt  NE  ) 
complax*16  g,gaamiag,sal  NE  >NE  ),curmttME  ) 
complex*16  sa2l NE  ,NE  I  ,curr2l  NE  )  ,*tl( 2,2) 
complex*  16  sctNE+l,NE+l I 
raal  error 
integer  maxitr 

character* 5  prtitr ,ganint,refmc,yas 
charactar-*26  vers 


common/unknwn/x , jacob 

common/ no twrk/ sn.cn 

common/waves/a >b 

cossson/genrtr /g  >  gammeg 

common/ an tana/sa , cur m t 

comaon/opt ion/error , maxitr, prtitr ,genint 

common/macpar/mcheps > sqrsps 

cotnmon/vrs  ion/vers 

co«Mon/pateh/stl,M2,curr2>rafme 

yaa»'YES' 
write!  6,10)  vers 

10  format! •  ******  SANE-PODCON  *****  Version  :  • ,A26,//) 
writa(6,20)  prtitr , error, maxitr 
20  format! '  Iteration  dump  •  ' ,A3/ 

A  '  Error  for  tha  norm  of  system  *  * »E12.5/ 

t  '  Maximum  number  of  iterations  *  ’,14) 

write)  6,30)  mcheps 

30  format! '  Machine  Epsilon  ■  ' ,E12.5) 
e 

e— - initial  guess 
e 

if  tganint.aq.yas )  than 
write! 6,40  I 

40  formatl '  Initial  guess  ware  calculated. ’  ) 

else 

write) 6,451 

45  formatt 1  Initial  guess  ware  inputted. '  ) 
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and  if 

write!  6  >50) 

50  FORMAT! 1 -  The  initial  guaas  for  the  solution  vector1  ) 

write!  6,60)  I  i>xl  i  )  ,cdabst  x(  i  )  ), phase!  xl  i  )  ),i*l  ,NE  ) 

60  format! 3X, 'XI ' ,12. ' I  a  I  * ,E12.5> * > * ,E12 .5, *  )  ■  * ,F8.5, '/' ,F8.3,/ 

*  ,'*',4ax,' _ _•> 

write!  6 ,70)  NE  lOiwig  >cdabs!  gamwag  1  .phase!  gawag ) 

70  format! // , 15X , ' -  Network  Parameter  - '  / 

6  *  The  antama  network  containas  ’  >13 , '  elements1/ 

<  '  reflection  coefficient  of  the  generator  *  (',£12.5,’,', 

«  E12.5,  ’  )  »  ' f fB.5> '/' >f8.3/> '♦'  .BOx.  * _ '  ,// 

l  ’  Tha  S-aatnx  of  tha  antenna  array'  ) 

do  85  i»l,NE 

write! 6 >80)  I  sal i > j  )  > i=l >NE  ) 

80  format! 3x>8(  'I  ' ,f5.3> V >f5.3> * > 

85  continue 


)) 


call  calcsclsc,x) 
write!  6  >200) 

200  formatl//>'  The  S-natrix  of  initial  feed  network') 
do  210  i=l.na*l 

write! 6  >220  )  I  Sc! i > i  )> )*l>ne«l ) 

220  format!  3x>9!  '!  *  >f5.3 , '  ,  '  ,  f5.3>  *  ) '  )  ) 

210  continue 

write!  7,230)  1 1  sc!  i>jl>i,i>i3l,NE9l),i*l,NE+l) 

230  format! 2f 10.5, ‘  Power  Divider  Feed  *,il,'>'>il) 

if  I  refmc.eq.yes )  then 
write! 6, 90 1 

90  format!//'  The  desired  currents  at  the  patch  edges  '  ) 

write! 6d00)  I  i  .curmtl  i  ),cdabsl eurmtl  i  ) ) ,phese(  curmtl  i  ) ), 
t  i*l,ne) 

100  format!'  II ' >i2. '  )  «  I ' ,el2.S, ' , • >el2.5, •  )  »  • >f8.5, '/' >f8.3,/ 

A  ,  '  ♦ '  ,48x  > ' _ ') 

write! 6,1151 

110  format!//'  The  desired  currents  at  the  reference  plane  '  I 

write! 6 >100 )  i i>curr2l i  ) >cdabsl currZl i  )  )>phase( curr2! i  )  )>i>l,NE  ) 

else 

write! 6 >120 ) 

120  format! //,'  The  desired  currents') 

write!  6 >100 )  I  i>curmt(  i  )>cdabst curmtl  i  )  )> 

•  phase!  curmtl  i  > )  >1*1  >NE  ) 
endif 
return 
end 
e 

c 

subroutine  calcf<f>x) 

c  calculate  the  nonlinear  equations  16.3-6) 


parameter!  NE  =8 ) 

eomplex*16  f(NE  )>x!NE  )>felrat>a! NE  )>b!NE  )>sn>cn>zl ,sc(NE+l  ,NE  +  1 ) 


common/netwrk/sn  >cn 
common/waves/a  >b 
c 

c — — get  scattering  matrix  of  coasaensation  network 
o 

cell  calcsclsc>x) 


c - -calculate  equation  16.3-6) 

c 

do  100  i*l,ne 
zl  *  0.0 
falmt  *  -a! i ) 
do  110  j*l,ne 

falmt  s  feint  *  sc! i > i  )*bl j  ) 
zl  ■  zl  -  sc!  ne*l > j  l»bi j  ) 

110  continue 

fl  i  )  »  faint  ♦  sc!  i  >ne+l )»( cn*x(  ne  )+sn*zl ) 
100  continue 
return 
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e» 

o 

e 

e 

o 

c* 

e 


function  falaat!  i  ,x  I 

calculate  i-th  nonlinear  aquation 


paraaatar  (NE*8> 

coa*>lax«16  x(NE  )  ,a(NE  )  ,b(NE  )  >sn.cn>zl 
coeplox«16  foleat >sc( NE+1 >NE*1 ) 
coaaon/natwrk/«n»cn 
cossnon/wavos/a >b 


*1  »  0.0 

call  getscili»sc>x) 
falaat  *  -a I i ) 

do  100  j»l,NE 

falaat  »  falaat  ♦  *c<  i  »  i  )«bl  j  ) 
zl  »  zl  -  sc(NE+l>  j  )*bl  j  ) 

100  continue 

falaat  »  falaat  a  ac(  i  >NE+1 )»( cn*x(  ME  >+sn*zl ) 
return 
and 
e 

c*»e«»»saeaea»ea*»»a*««*»«»a*»aa 

o 

subroutine  calcscl sc ,X ) 

o 

o  calculate  S-aatrix  of  paster  divider  casmensation  network 

o 

c 

paraaatar  I NE*8  I 
conMon/iii/ii  ^nx 

coaplox*l6  xl  NE ) »sc( NE+1 ,NE+1 ) ,kl NE*2 )  ,p( NE*2 ) .one, zero 

coaplax«16  sctaap,sctnp,plc,pld,piovr2 

integer  i >1 la  .level 

integar  il , jl.iif NE ,NE  1  ,«xl NE ,NE  ) 

real*8  batalcidzero 

one«< 1.0,0.01 
zero*!  0.0,0.01 
dzoro  *  0.0 

be talc  a  3.141592654  /  4.0 
plc=  cdoxpl  dcnplxl dzero , -batalc ) ) 
pld*  cdexpl dcnplxl dzaro  ,-2 . 0*betalc ) ) 
piovr2  *  10. 0>1. 01 
c 

c — — - compute  poster  ratio  and  phase  shifts 

°  do  100  1*1 .NE-1 
lo  »  ( 1-1  )*2tl 
la  *  1»2 

kl  la )*cdabsl x( 1 ) ) 
if  (kl  la  l.ne.zaro  >  than 
pi  le)*xl  1  )/kl  le) 

alsa 

pi  le  )sone 

i 

kl  lo  )*cna-k(  le  I 
p(  lo )>one 
100  continue 

laval*intl  log!  float! NE  )  l/logl  2.0 1+0.5) 
c 

e — — input  port  S-paraater 

sclNE+l,NE+l)=zero 

c 

c - SCvi  and  SCv  in  (6.2-6) 

do  200  i*l,NE 
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OOO  000  000 


210 


200 


•ctenpsone 
do  210  1*1.1 

il*iil  i.l)  ...  , 

»cta«p**cta»i>*cd*qrtt  k(  il )  l*p<  il  )*plc*piovr2 
continue 

tel  i  >NE*1  )*sctamp 
sc(NE+l,i  l=sctemp 
continue 


- SCI i.i ) 


do  300  i*l,NE 
tetMo’MPO 
do  310  1*1, level 
il*ii< i.l ) 

ll*il*eod<  il.2  )i*od(  il+1,21 
*cimp*kt 11 )*pl il  )*p( il  )*pld 
11*1-1 

do. 320. 11*11,1.-1 

sctm>**ctmp*k( il  )*p< il  )*p( il  )*plc*plc*( -1.0, 0.0  ) 

320  continue 

scte*e>**cta*p+*ct*p 
310  continue 

ac( i,i  )*scteep 
300  continue 

- SCii  and  SCiii  in  1 6.2-6 ) 

ibgn*N2/2*l 

iend=NE 

do  400  i*ibgn,iend 
do  410  j*f,NE/2 
*ctaep*one 
setup *one 
do  420  1*1, level 
il*iil i,i ) 
jl*iitj,l) 

»ct«ap**ct«;»*k(  il  )«k(  jl ) 

•ctep*setnp*p( il l«pl  jl )*plc*plc*( -1.0, 0.0  I 
420  continue 

•cl i , j l*cdsqrt( seteap  )*sctmp 
•cl  },t l*scl i,  j  I 
410  continue 
400  continue 

- SCi  in  (6.2-41 

iend*NE/2-l 
do  500  i*l,iend 
jbgn*i*l 
iend=NE/2 

do  510  j*jbgn,jend 
lbgn=mxl i , j I 
il*iil i  ,lbcri ) 

jl*ii<  j .lbgn  I  .... 

sctmp*-cdsqrt(  k(  il  )»kl  jl  J  )*p(  il  )*p(  jl  )*pld 
wend* lbgn- 1 
do  520  n*l,aand 
il*ii(  i,e) 

«ctep*»ctnp«cdsqrt( k(  il  )«k(  jl  I  )*pl  il  )»p(  jl  )*ple»ple*<  -1.0,0 1 
520  continue 

scteep*sct«p 
lbgn*lbgntl 
do  530  1* lbgn, level 
il*ii( i.l  I 

ll*il+*iod( i  1,2) -modi il+1,2) 

«ctwp*k(  11  )«p(  il  l«p(  il  )*pld 
nend=l-l 
do  540  e*l,mend 
il*ii( i»e t 

«ctepi*cb*P<Kxisqrt(k(  il  l«k(  jl )  )«pl  il  )*p(  jl  )«plc*ple*(  -1.0  J 
540  continue 
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k. 


« 


l: 


I 


sc  temp*sc  iMp<sc  tep 
5S0  continue 

•cl  i»  j)*sctamp 
scl j>i  )* sc temp 
510  continue 
500  continue 
c 

c — — SCiv  in  16.2-6) 
iend*NE-l 

do  600  i-NE/2+liiend 
jbgn*i+l 
jend=NE 

do  610  jsjbgn.jend 

lbgn*mx( i-NE/Z> j-NE/2 1 
il=ii<  i  .ltxyi ) 

jl»ii< i»lbgn )  .... 

sctmps-cdsqrt!  k(  il  )*K(  jl )  )*pt  ll  >*pl  jl  )*pld 
raendslbgn-1 
do  620  m«l  ,mand 
il=iil  i >n ) 

sctmps*ctrap*cdsqrt(k(  il  )»k(  jl )  )*p<  il  )*pl  jl  )*plc*ple*<  -1.01 
620  continue 

sctempssctmp 

lbgn»lbgn*l 

do  630  Islb^i, level 
il*ii( i»l ) 

ll*il+«edl il >2)  modi il*l,2) 
sctapskl  11  )*p(  il  )*p<  il  )*pld 
mand*l-l 
do  640  m»l(mand 
il*iit i>e ) 

jl*ii(j,m)  ... 

sctmp*sctmp»cd»qrt( k!  il  )«kl  jl )  )*p<  il  )*pl  jl  )*plc*plc*l  -1.0  ) 
640  continue 

sctempssctemp'tsctmp 
630  continue 

scl  i  ,  j  )=sctesp 
scl j.l  )»sctemp 
610  continue 
600  continue 
return 


e 

c< 

c 

c 

c 

c 

c< 

c 


subroutine  getscil i,sc,X) 

fl*t  i  -th  row  of  compensation  network  S -matrix 


parameter  (NE-81 
common/ i i i/i i ,nx 

complex*16  x!NE  ), scl  NE*1,NE*1>, kl  NE*2  )  ,p<  NE*2  >, one, zero 
cotnpiox*16  setemp,  setup, pic  >pld,piovr2 

real«8  betalc>dzero 
integer  i , j ,l,m,lo,le, level 
integer  il,jl>ii(N£,NE),ax(NE,NE) 


one=( 1.0, 0.0 ) 
zeros!  0.0  >0.0) 
dzmra  *  0.0 

&iovr2  »  (0.0,1.01 
•talc  «  3.141592654  /  4.0 
pic*  cdexpldcmplxl dzero>-betalc ) ) 
pld*  cdexp(dcmplx(dzero,-2.0»betalc) ) 
c 

c — —compute  power  ratio  and  phase  shifts 

''  do  100  1*1  >NE-1 

lo  s  (1-1)»2*1 
le  a  1*2 

k(  le  )scdabs(  x<  1 1 ) 


I 
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L 


00  0  OOO  OOO  000  000 


if  ( kt  la). na. zero)  than 
plla)*xll)/k(la) 

•1m 

pi  la  >*ona 
and  if 

k( lo J*ona-kl la ) 
pi  lo)»ona 
100  continua 

laval*int(  log!  float!  NE  )  )/log(  2.0  )*0.5) 
— SCI  i  ,NE  )  and  SCI  ME,  i) 


do  200  jrl >na 
scta«p=ona 
do  210  131 ,laval 
jl*iil 3,1 ) 

set*Mp3sctaap«cdsqrtlkt  jl )  )*p(  jl  )«ple*piovr2 
210  continua 

scl 3  ,na+l )3sctamp 
scl na*l , j  )asctamp 
200  continue 


- SC(NE«1,ME*1) 

if  I i. gt.no)  than 
SC|NE+1,NE+1 )=Zwro 

•Isa 

- SCI i,i ) 

sctMp'zaro 
do  300  l3l,laval 
il*iili>ll 

ll3il«asod(  il,2)-mod(  il+1,2) 
sctmpakl 11  )«pl il  )*pl il  )«pld 
11*1-1  . 

do  310  11*11,1,-1 
il3iil i,ll I 

sctm>*sctmp*kl  il  )«p(  il  )<*p(  il  )*plc*pic*l  -1.0 ) 
310  confirm* 

setetss « sctmp 
300  continua 

scl  i,i  I  =sctatnp 

— — SCii  and  SCiii  in  16.2-6) 


410 

400 


if  I  x.ga.na/2+1 ) 
jbgn*l 
jand=NE/2 

•Isa 

jagn*NE/2*l 

jcnd=NE 

•ndif 

do  400  j=jbgn,jand 
sctemp=ona 
sctmp =ooa 
do  410  l*l,laval 
il«iii i,l ) 
jl«iil  j,l) 

setaapasctaapakl il  )«kl jl ) 
setnp»sct«a»*pl  il  )*pl  jl  )*ple»plc*l  -1.0) 
continua 

scl  i ,  j  )*cdsqrtl  sctaap  )*sctm> 
continua 


- SCi  in  16.2-4) 


if  li.la.NE/2)  than 
do  500  j*l,NE/2 
if  I i.na. j )  than 
lbgn=«tx(  i ,  j  ) 
il*iil i,lbgn) 
jl3iit j,lbgn) 

setup* -cdsqr tl kl  il  )*kl  jl )  )«pl  il  )»p<  jl  )*pld 
aand*lbgn-l 


I 


| 
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JO  O  O  O  JO  000 


do  510  ■«!>— nd 
il*iit  i.a) 
jl»iitj,a) 

*ctnp**cfap* 

*  cdsqrtt  K(  il  )*kl  jl  I  >*p<  il  )«p<  jl  )*plc*plc*(  >1.0) 

510  continue 

SCtSMP=SCtsp 
lbgn*lbgn+l 
do  520  lalbgnf level 
il*iit if  1 ) 

ll*ilniiodf  il  ,2  I -nodi  il  +  1,2 ) 
cctnp=k(  11  )«p(  il  )»pt  il  )*pld 
mond=l-l 
do  530  m*l>mond 
il*ii( i >n ) 
jl*ii(  jtil 
sctap^sctnp* 

<  cdsqrtf  kl  il  )«kl  jl )  )«p(  il  )*pt  jl  )*plc*plc*l  -1.0) 

530  continue 

sctanp*sctsap+sctmp 
520  coniine* 

set i» j  )3sctsnp 
set  jti  I* sctsap 
sndif 

500  continue 

•Iso 


— SCiv  in  16.2-0) 

do  600  j»NE/2+lfNE 
if  t  j.na.i  )  than 

lbgn>mx( i-NE/2, j-NE/2 ) 
il*ii( itibgn) 
jl*iil  j»lbgnl 

sctmp*-cdsqrttk(  il  )*kf  jl )  )«p(  il  )«pt  3I  )*pld 
nand*ibgn-l 
do  610  a»lfaand 
ilsiilita) 
jl»iit  j,a) 
sctap*setnp* 

t  cdsqrtt  kt  il  )«kt  jl )  )*p(  il  )«pt  jl  )«plc*plc*l  -1.0) 

ilO  confirms 

setaapssefap 
lbgn*lbgn*l 
do  620  lalbgnflaval 
il*iit ifl ) 

ll*il«nodt  il  ,2  )-mod(  il*l,2  ) 
scfmp=kl  11  )«pl  il  )*pt  il  )#pld 
«and*l-l 
do  630  ■■IfMand 
il=»ii< i f« ) 
jl3iif  jfO) 
sctmp=sctnip* 

t  cdsqrtt  kt  il  )*kl  jl )  )*p(  il  >»p(  jl  )*plc*plc*l  -1.0 ) 

630  confirms 

sctefi*>=sctsfflp+sc'tmp 
620  confirms 

set  if  j  l^scti 
set  jf i  )=scf« 
sndif 

600  continue 
sndif 
•ndif 
rsfum 


4 


subroutine  ealeii 

calculafa  index  sysfssi  us  ad  in  power  dividar  network 


c 


P* 


remoter  ( NE*8  I 

1 
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co— on/iii/ii  .ax 

intigtr  iilrwintlimlMiiwIiltvtl 


level*intl  log(  float! NE  )  )/log(  2.0  1*0.5 ) 
i  court*!  NE-1  )*2 
do  100  l*l«v«lil>*l 
do  110  i*l.NE 

&i€  i >1  )*icount-(  na-i  )/!  2** 1 1-1 ) ) 

1X0  coniinLM 

icount*icount-NE/l 2—1 1-1 1) 

100  continue 
pa  turn 
and 
e 

e 

subroutine  calc— 
c 

c  calculate  index  syst—  usad  in  powar  divider  network 
a 

para— ter  (NE*A ) 

integer  iilNE.NE).— INE.hS  l.level 
co— on/iii/ii.— 


level* inti  log!  float! NE  )  )/log(  2. OHO. 5) 


130 

120 

110 

100 


c 

c* 

c 


c 

c 

c 

c* 

o 


do  100  kalevel.2.-l 
jbgn*2*»k/4*l 
iend*2»«k/2 
xstep*2**( k-1 I 
do  110  i*l.na/2»istap 
lbgi*i 

lend*i*<  jend-  jbjpi  > 
do  120  l*lbgn.land 
do  130  i*ibgn.iend 
—i l»j>i-l)*k-i 
—I j*i-l»l)*k-l 
continue 
continue 
continue 
continue 
return 


subroutine  cal iac 
para— ter  INC  *8! 

calculate  Jacobian  of  the  nonlinear  equations 


coeplex*16  xl Nt )  .jacob! NE .NE  ) 

coaplex*16  flNH l.rtenpxINE  ).itenpx(NE ) 

real«8  restep. im tap. rate— .iateap. re jac.iajac .zero. typx 

real  acheps.sqreps 

c  o—o n/cnkrwn/x  >  jactA* 

c  o— on/— epa  r/— chap .  sqreps 


o - save  current  solution 

c 

zaro*0.0 
typx*0.01 
call  ealcflf.x) 
do  1000  i*l,NE 
rtampx!  i  )*x(  i  ) 
itenpx! i  )*x( i ) 

1000  continua 
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c - 

o 


e 

c — 
o 


e 
o— — 
o 


2100 

2000 


m  finite  difference  in  real  part 
do  2000  j«l,NE 

if  I  daosl  draall  xl  j  ))  hgt.typx)  than 
res  tap*draal  I xt j  )  )**qr#p* 
retaap*dreal< xl  j ) )♦ res tap 
restap*reteap-dreell xl j ) ) 
alsa 

restapadsi jpd  typx,  draall x(  j  >ll«sqreps 
andif 

rteapxl  j  )*xl  j  )*deq>lxl  raatap ) 

*u«a  finita  dif farartca  in  imaginary  part 

if  (dabsldiaegtxl  j  I  ll.gt.typx)  than 
imtapsdixagt  x(  J  )  l*sqreps 
iateap*diaagt  xl  J  )  Him  tap 
iastep»iataap  diaagl  xl  j  H 

ala* 

intap>dtiml  typx ,diaag(  xl  j ) )  )*sqraps 
andif 

itaopxl  jl«xl  j  )*deaplxl2aro,ia*tap> 

-estiaete  Jacobian 
do  2100  i-l>NE 

re  jac«<  draall  falaatl  i>rtaapx  I— f f  i  )  )  l/rastep 
iajac>( draall  falaatl  i.itempx  )-fl  i  )  I  )/iastap 
jacobl  l»  j  ladaaplxl  rejact-iajac ) 
continue 
rtaspxl  j  )*x(  i  ) 
itaapxl  jl*xl  j) 
eontinua 
return 


a 

c« 


subroutine  newton 
o 

e  aodifinad  newton's  aathod  driver 
o 

caaaaaaaaaaaaaaaaaaaaaaaaaa 

a 

paraeeter  (NE4I 

com  lax*  16  xl  NE  )  >  jacobl  NE  >NE  )  >  dal  taxi  NE  ) , f  I  NE  ) 

real  error (  norar  phase 

integer  aexitr,  iter 

character*!  yes  >  prtitr»  ganint 

logical  alphck 


cowaon/unkr wn/x  >  jacob 

coaaon/opt  ion/arror ,  aax  itr  »prt  i  tr  >  ganint 
coaaon/dpinfo/raslop .alphofc >deltax , f 
consort/ result/iret , i tar 
e 

e— -initialise 
o 

iret*0 

itar»l 

Q 

c - is  error  lass  than  specified 

c 

100  call  caleflf>x) 
checkv*dtore<  f  >na ) 
if Icheckv.lt.error )  goto  W 
if! iter. ne.l. and. prtxtr.eq. 'YES' )  than 
witsItilOl  iter.  chackv >  alpha 

10  foraatt 1  iteration) ' >i4»3x> 'non*  *  '.el2.5>3x. 'alpha  *  ',el2.5l 
andif 
e 

c— -gat  Jacobian 

c 

call  caljae 
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c— --got  next  (tap 

do  200  ial.no 

del taxi i  )»-f ( i ) 
200  continue 

— -aolvo 


coll  aatsollne.deltax,  jacob) 

- gat  dipping  factor 

call  damped!  alpha  I 

- compute  next  ostioaticn  to  tha  solution 

do  300  i*l.ne 

x! i l*xl i  )4alpha*deltax( i  I 

if  I  I  i.na.NE  hand.  (dcabslxl  i  )  l.gt.1.0  ) )  irat*2 
300  continue 


- show  tha  result  of  this  iteration 

write( 10.201  iter,  checkv.  alpha 

20  format! *  iteration: ' >i4.3x. ‘norm  •  ' >al2.5.3x. ‘alpha  a  '.al2.5) 
xr i tat  10.21)  xtnal 

21  fomotl  2el2.5) 

if  (irat.aq.2l  than 
write! 10.30  I 

30  fonaet!  1  abs(  x(  i  I )  >  1 .  Program  tanainatad  '  ) 
goto  999 
endlf 

—is  iteration  lass  than  aax 


itar  ■  itar  ♦  1 
if  I itar. It. maxi tr)  goto  100 
i rat»l 
999  return 
and 


subroutine  daapadt alpha ) 
compute  damping  factor  alpha 


o 


e 


parameter  I  NE*8 ) 

real  lands,  u.  1.  minalp.  norm,  laftsd.  rigtsd.  If.  insqrt 
logical  alphok 

complax*16  xl  NE ) , jacobt NE  »NE ) .xtampl  NE ) , dal taxi NE ) , f f  NE ) , f tempi  ME ) 
complex*16  tamp] , inelop 

comnon/inkrMn/x .  jacob 

common/ct>  info/ res  lop.  alphok  .dal  tax,  f 

- — initialize 


slphaal.O 
landa«l . 0a-02 
u*0.5 
1»0.1 
foalflf) 
alphoka.trua. 
it«0 

inslopal 0.0,0.01 


do  15  ial.NE 

tamp  j»l  0.0, 0.0 1 
do  10  3*1 ,NE 
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e 

c* 


10 


taapj*taopi+canjg(  ft  j  )  )• jacobl j,i ) 
continue 

inalopajnalcp i teap  j«daltaxl i ) 

15  cent inua 

reslop  *  real!  inslop ) 

20  do  30  i«l.NE 

xteapi i  )*x( i  )+alpha*daltax( i ) 

30  continue 

call  ealef I ftaap,xtaap ) 
fpr*v»lf ( f  teap ) 

if  I  fprsv.lt. fo*alpha*landa*res lop)  ra turn 
if  (alpha. It. 1.0e-05)  than 
alphofc*.  falsa, 
ratum 
andif 

altaapa-raslop/l 2*1 fprev~fo-reslop ) ) 

alprev*alpha 

f2prav«fprav 

if  (altaap.lt.0.1*alpha)  than 
alphe*0 .  l*alpha 

alsa 

alphaaaltaop 

andif 

it*it+l 

if  I  it. It. 10)  goto  20 
alphoka. falsa, 
ratum 


c 

c 

o 

c* 

e 


subroutine  gatour 

coaaxita  the  coapansatad  currants  for  verification 


parameter  (NE*8I 


coaplax*16  x(NE  )>  jacobl  NE  >NE  )  >sal  NE  ,NE  )  ,eurmt(  NE  I 
com>lax*16  q>gawaag»sa2l NE ,NE  ) ,curr2(NE  ) >stl(  2,2  >,*c(NE*l,NE*l ) 
complex*16  cntINE),  snttNE  >NE  ),te«nps(  NE*2,NE*2  ) , tempcl  NE*2 ) 
complex* 16  calcur(NE)*  ap(NE),  bp(NE) 
character*!  ref  me >  yes 


c  oaacn  AsdUTHn/x  >  jaceb 
conaon/ganrtr/g  ijiang 
coaaon/antene/'se  icurmt 
coaaon/petch/stl  >  sa2  »curr2  >  ref  me 


yes  ■  'YES’ 
c 

c — — initialize  taape  and  taaps  to  zero 
c 

do  100  i=l,NE*Z 

tempcl i)  *  1 0.0 >0.0  I 
do  110  j*l >NE*Z 

teepsfiaj)  »  (0.0,0.0) 

110  continue 
100  continue 
e 

e— — get  sn  and  cn  of  the  coapensation  network 
e 

call  calcsc(so>x) 
do  200  i*l,NE 

ont(i)  ■  xINE)  *  scli>NE*l) 
do  210  j*l >NE 

sntli,j)  *  set  i>i)*gs*amq»scli,NE*l)»sclNE*l,j) 

210  continue 
200  continue 
e 

c - set  up  the  aetrix  equation  describing  the  entire  network 

do  300  i«l,NE 

t amici  i  I  >  cntl  i  ) 
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do  310  j«l,NE 

teams!  i^NE.J^tC)  «  -wli.jl 
taapfliij)  *  -snt!  i,  j ) 

310  continue 

imp-sA  i,NE  +  i  )  *  (1.0,0.01 
tempsINE^i.i  )  *  (1.0,0.01 
300  continue 
ne2  *  NE*2 

coll  matsol!  no2 ,  tempo >  tawgs ) 
do  400  i»l,NE 

calcurl  i  I  *  iwoel  i«-NE  )-tampcl  i  I 
400  continue 
write!  4,20) 

20  format!//'  The  calculated  currants  at  the  refaranca  plane  ') 
do  500  i»l,NE 

wri tel 6,10 )  i, calcurl i ),deabsl calcurl i  1 1 , phase! calcurl i  ) ) 

10  foreetl 3x, 'll  ' ,i2» '  I  *  I ' , F10.5, ' , ' >F10.5, ' )  * 
t  F10. 5, 1  /  * ,F10.51 

500  continue 

if  I  ref  me.  eq.  yes )  than 
do  600  i«l,N£ 

bp  I  i  )  a  I  teepcl  i  1-stll  1,1  )»tampc!  NEti )  l/stll  1,2 ) 
apli)  *  stll  2>1  )*tewpc(  NE+i  )+stl(  2,2  )*bpi  i  I 
calcurl  i )  a  apli  J-bp!  i  ) 

600  continue 

write! 6,30 1 

30  fora— tl//'  The  calculated  currents  at  the  patch  edgas  ') 

do  700  i«l,NE 

write! 6,10 )  i, calcurl i ),dcabs I calcurl i  ) ), phase! calcurl i ) ) 
700  continue 
endif 
return 
and 
e 

c**«»»**»»»*****»*»»»*»*«»***»*»**i 

c 

subroutine  disply 
c 

c  output  the  results 
c 

o 

parameter  INE><) 

complax*16  xINE),  jacoblNE.NE  1  ,sc(  NE  +  1  ,NE*1 ) 
common/ unknwn/ x , jacob 
common/ resul t/irat , i ter 
co— aon/tiee/itima l,itima2 


write!  6,4) 

4  format! ’1  - RESULT - •///) 

if  liret.eq.O)  then 

write!  4, 5)  iter 

5  format! *  Program  was  successful. 

else 

writel6,6l  iret 

6  format! '  program  was  unsuccessful, 
endif 

call  timeckl itimeZ) 
rtime  «  f loat( i timaZ-itimel 1  /  100.0 
write! 6,7)  rtime 

7  formatl/'CPU  time  ■  ',f6.Z,'  seconds  ' 
write!  6,10) 

10  format!  '-The  retired  power  divider  specification  '/) 
write! 6,20 )  ( i »xl i  )  .cdabslxl i  ) ), phase! xl i  ) )  ,i=l,ne-l ) 

20  format!5x, 'XI ' ,i2, * )  »  I ' ,el2.5, ' , ' ,el2.5, '  )  «  ' ,f®.5, V ,f8.3,/> 
t  '♦'.SOX.'  ' I 

write! 4,40)  xl ME  I  .cdabsl xl NE  ) ) .phase! xl NE  > ) 

40  format!/'  Generator  »  ' ,f8.5, ’ , ' ,f8.5, ‘  *  ‘ ,f8.5, '/' ,f8.3,/. 


Nunber  of  iteration  *  ’,i4) 
iret  *  '.ill 


) 


»41X» 


) 


call  calcsclsc.x) 
writal 6,200 ) 

200  format!//,'  The  S-matrix  of  modified  feed  network'  ) 
do  210  isl.r— ♦! 

write! 6,220)  ISC! i »j  ) ,  j«l ,ne*l ) 
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220 

210 

230 


240 


fonaetl  3xi  9t  '(  '  »f5. 3 » '  >  *  ,f5.3 . '  )  •  )) 
continue 

writ*! 3(230)  ( < sc! i.j  ),i t j,  j»l,NE*l ) .i«l,NC«l ) 

FORMAT! ZtlO.S, •  Power  Divider  Feed 

cell  timeck! itime2 ) 

r tints  *  float! itime2 )  /  100.0 

writs!  6(240  )  r tints 

format!///'  Total  CPU  time  a  ‘ ,f 6.2»*  seconds  ' ) 
return 


e 

e* 

c 

e 

e 

c 

am 

a 


real  function  lfif) 
norm  squared 


ter  !NE*3) 


10 


complex*!*  f I NE ) 
lfsO.O 

do  10  i«l(NE 

lfalftf!  i  )*conjgl  f I  i ) ) 
continue 
If *0. Saif 
return 


o 

c* 

o 

o 

o 

c 

c* 

o 


real  fine t  ion  drtorm  ( vector  >n ) 
daUble  percision  norm 


camplex*16  vector!  n) 
roal*8  tamp 


t  strip  *0 . 0 

do  100  i>l,n 

tamp^cdabsl  vector!  i  i  )**2+ts«o 
100  continue 

dnorm=dsqrtl temp ) 
return 
end 
c 

e 

subroutine  matsol  (n(btpass) 
c 

c - matrix  aquation  solvar 

c 

c***««**«»****«*****»*******»*»»»i 

c 

parameter  (NE*8) 

complax«16  pass!4*NE*NE  )  ,a(  2*NE  ,2*NE  )  ,tenp,bl  2*NE  ) 
real  di>mex 

integer  ipivotl 2«ne-l ) >n 
o 

o - transfers  the  information  in  the  vector  pass  to  the  matrix  a 

o 

do  100  iai,n 
do  110  j>l(n 

a!  j(i  lapassl  I  i-1  )*n+i  ) 

110  continue 
100  continue 
e 

c - factor  the  matrix  into  an  upper  triangular  matrix  using 

e  implicit  pivoting, 
o 

do  200  k*i,n-i 
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c — — pmrfom  th*  implicit  pivoting 
o 

nxifl ,  o 

do  210  irow*kfn 

di*dcabsl  al  irow.k ) I 
do  211  i=k*l,n 

if  Idi.lt.dcabslat  iroMfj)))  di^dcabslal  irow>3 ) I 
211  continuo 

if  tdnbstal  irow.k  ll/di.gt  .Max)  than 
ipivot  lkl*irow 
eax*dcabs(al irow>k >)/di 
andif 

210  continue 
e 

c — — «*ap  the  two  rows 
o 

do  220  isk»n 
iam>alkiil 
alk,i  l=a(  ipivotlk  l.i  ) 
al  ipivotlk  I >i  )>taap 
220  continue 
e 

c - zaro  tha  k  coluan 

o 

do  230  i*k*l »n 

tcap*-al  i  >k  )/al  k  *k  I 
•liikl'taw 
do  231  j*k«l>n 

al  i  >  3  >»al  i  *  3  )*t«ap*at  k  >  3 1 
231  continue 
230  continue 


200  continue 
e 

c — —establish  tha  corresponding  b  vector 

C  do  300  k*l»n-l 
c 

c- - pivot  tha  b  vector 

c 

teinp*bl k  I 
blk  )*bl  ipivotlk  I ) 
bl  ipivotlk  I  l*tenp 
c 

c — -perform  tha  row  Multiplications  on  b 
o 

do  310  i*k*lfn 

bl  3  l=bl k  l*a<  3  >k  )*bl  3  I 
310  continua 
300  continua 
e 

c— backsolva  for  tha  solution  vector 
c 

btn)=b(n>/a(n.n) 
do  400  iis2>n 
i=n+l-ii 
tempabl i  I 
do  410  i=i*l>n 

tenpataap-bl 3 )*a( 1,3) 

410  continua 

bl  i  l>tesp/a(  i»i  I 
400  continua 
return 
and 
a 

caa»*»«»»»»**»**s*«******* 

c 

function  phase! x I 
c 

o  calculate  phase  of  complex  variable  x  in  dagi 
c 

c**»*****»***»*«**»*«»**»* 

c 

comp  lax* 16  x.zaro 
res 1*8  rZd 
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rid* 180 . 0/3 . 141592654 
za ro«0.0 


e 

c* 

e 

o 

o 

c 

c* 

e 


if  ( fdraall  x  l.eq.zaro  ).and.  (diwaglx l.eq.zaro  I )  than 
phaaa*0 . 0 

•Is* 

phase*sngl(  datan2(  dinayl  x )  rdreall  x )  )*r  2d 1 
end  if 
pa turn 
and 


raal  function  dcabslx) 

doubla  parcision  complex  absolute 


complex* 16  x 

deabs*sngl(  cdabsl  x ) 1 

return 

and 
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